Final Report on NASA contract NAS8- 36169 


Determination of Atmospheric Moisture Structure and Infrared Cooling 
Rates from High Resolution MAMS Radiance Data 

for the period of 

7 November 1986 to 18 September 1991 

submitted by 
W. Paul Menzel 
Christopher C. Moeller 
William L. Smith 

Cooperative Institute for Meteorological Satellite Studies (CIMSS) 

University of Wisconsin 
1225 West Dayton Street 
Madison, Wisconsin 53706 


December 1991 



Table of Contents 

I . INTRODUCTION 1 

II. RESEARCH ACCOMPLISHMENTS 2 

A. COHMEX Data Collection 2 

B. HIS Research Results 2 

1. HIS/FASCODE Comparisons 2 

2. HIS Retrievals 3 

C. MAMS Research Results 4 

1. MAMS Instrument Improvements and Data Comparisons 4 

2. COHMEX Data Archive 5 

3. MAMS/VAS Retrieval Products 5 

4. MAMS/VAS Fluxes 7 

5. MAMS/VAS Cloud Clearing Activities 9 

III. LIST OF PUBLICATIONS RESULTING FROM NAS8-36169 10 

IV. REFERENCES 13 


APPENDIX A - NASA SUPPORTED PUBLICATIONS 


15 





I. INTRODUCTION 

With the advancement of remote sensing from space as a tool for 
monitoring the earth's atmosphere and surface, new and better 
information is constantly becoming available for investigating 
meteorological processes and their spatial and temporal scales of 
preference in the earth- atmosphere system. Much information has been 
gained in understanding synoptic and mesoscale variability; the 
Multispectral Atmospheric Mapping Sensor (MAMS) is contributing to this 
understanding by looking at very small mesoscale (meso 7 : 2-20 km) 
horizontal variability in the atmosphere and at the earth's surface. 

The High- resolution Interferometer Sounder (HIS) is contributing by 
producing high vertical resolution profiles of atmospheric temperature 
and moisture. These two research instruments, flown on NASA U2 and ER2 
high altitude aircraft, are helping to identify advantageous design and 
deployment characteristics for future operational satellite programs. 

This program has applied MAMS high resolution data to the problem 
of monitoring atmospheric quantities of moisture and radiative flux at 
small spatial scales. MAMS, with 100m horizontal resolution in its four 
infrared channels, was developed to study small scale atmospheric 
moisture and surface thermal variability, especially as related to the 
development of clouds, precipitation, and severe storms. HIS data has 
been used to develop a high spectral resolution infrared sounding 
retrieval algorithm for producing high resolution vertical profiles of 
atmospheric temperature and moisture. The results of this program are 
that high spectral resolution soundings significantly improve the 
vertical resolution of atmospheric prpfiles, and that high spatial 


resolution data shows small horizontal scale variability is clearly 
present in the atmosphere and on the surface. 


II. RESEARCH ACCOMPLISHMENTS 

A. COHMEX Data collection 

MAMS and HIS were flown on NASA U2 and ER2 aircraft during the 
Cooperative Huntsville Meteorological Experiment (COHMEX) (Williams et 
al. 1987). A listing of the flights is shown in Table 1. Funding from 
the contract was used to support science teams for the HIS and MAMS, as 
well as instrument field teams for the HIS during COHMEX. Science 
support teams contributed to the experiment through flight go/no-go 
decision-making based on the objectives of the field experiment, and by 
specifying flight patterns and locations for MAMS and HIS data 
collection as well as coordinating data collection activities among the 
many data sources available during COHMEX. The data collected during 
COHMEX was the primary source of data for scientific studies under this 
contract . 

B. HIS Research Results 

1. HIS/FASCODE Comparisons 

Comparisons of HIS spectra to FASC0D2 (Clough et al. 1986) spectra 
have been under intense study. In general, it is concluded that HIS and 
FASCODE agree remarkably well, a tribute to the current state of 
spectral line files and line-by-line code. However, several 
discrepancies have been identified and are reported in Re vercomb et al. 
1990. A copy of this paper is included in Appendix A. 


2 . 


HIS Retrievals 


A significant effort has been made to develop efficient retrieval 
algorithms for producing atmospheric profiles of temperature and 
moisture using HIS spectra. Atmospheric profiles for COHMEX days of 
June 15 and June 19, 1986 have been produced. These days were chosen 
because good quality HIS data of contrasting atmospheric conditions 
(boundary layer cumulus development on June 15; clear conditions on June 
19) was available. Statistics of HIS atmospheric profiles, along with 
comparisons to simulated atmospheric profiles (from HIS spectra) for the 
HIRS and GOES I/M instruments are given in Huang et al. (submitted). In 
this work, it is shown that the HIS atmospheric profiles are 
approximately five times better in vertical resolution and two times 
better in accuracy than HIRS and GOES I/M atmospheric profiles. The 
retrieval methodology is also transferable to the retrieval of 
atmospheric trace gas concentrations . Future EOS infrared sounder data 
has also been simulated and studied using HIS spectra (Huang et al. 

1991). Efforts to improve the retrieval methodology are ongoing. 



Table 1. Listing of MAMS and HIS flights during COHMEX. 


Date 

location 

MAMS? 

HIS? 

Comments 

6-15-86 

COHMEX 

Y 

Y 

both on U2 

6-17-86 

Local ocean 


Y 


6-18-86 

COHMEX 

Y 



6-19-86 

COHMEX 

Y 

Y 

repeat coverage 

6-20-86 

Chesapeake 


Y 

Tstorms 

6-23-86 

COHMEX 


Y 

Tstorms 

6-26-86 

COHMEX 

Y 

Y 

repeat coverage 

6-29-86 

Atlantic 


Y 

Ocean buoys 

7-1-86 



Y 

test hop 

7-3-86 

COHMEX 

Y 

Y 

both on U2 

7-4-86 

Chesapeake 


Y 

TMS mission 

7-05-86 

Atlantic 

Y 

Y 

both on U2 

7-08-86 

COHMEX 

Y 

Y 

both on ER2 

7-11-86 

COHMEX 

Y 

Y 

both on ER2 

7-12-86 

COHMEX 

Y 

Y 

both on ER2 

7-14-86 

COHMEX 


Y 

night flight 

7-15-86 

Lake Oneida 


Y 


7-17-86 

COHMEX 


Y 

aborted 

7-20-86 

COHMEX 

Y 

Y 

both on ER2 

7-21-86 

East coast 

Y 

Y 

ocean buoys 


C. MAMS Research Results 

1, MAMS Instrument Improvements and Data Comparisons 

Improvements to the radiometric performance have been made during 
the life of this contract. These include the use of a full axe blade 
mirror, 10 bit digitization, and improved external blackbody sampling. 
While hardware improvements were not funded by the contract, tests of 
MAMS data quality following hardware and procedural improvements were 
conducted. Science efforts under this contract benefitted through 
improved data quality. Comparisons to VAS , AVHRR, and HIS radiometric 
data were made. These improvements and the impact they have had on the 
data quality, as well as instrument intercomparisons are discussed in 
NASA Tech Memo 100352, which was produced in part under this contract. 



2 . 


COHMEX Data Archive 


Data from six selected MAMS COHMEX flights were archived at the 
Cooperative Institute for Meteorological Satellite Studies (CIMSS) at 
the University of Wisconsin. The archive procedure consisted of 
breaking the data up into separate straight line flight tracks, fixing 
identifiable data errors in housekeeping (calibration, time, etc.) and 
video (missing lines) data, and attaching earth location information to 
the data ("navigating" the data). This saved considerable data 
processing during subsequent use of the data. Table 2 lists the MAMS 
archived data. 


Table 2. MAMS COHMEX archived tape inventory. 


Tape 

Date 

Location 

Comments 

MC136 

6/15/86 

COHMEX 

popcorn cumulus 

MC138 

6/18/86 

COHMEX 

high pressure overhead 

MC140 

6/19/86 

COHMEX 

developing cumulus 

MC147 

6/26/86 

COHMEX 

mature Tstorm 

MC155 

7/03/86 

COHMEX 

synthetic navigation 

MC159 

7/05/86 

ATLANTIC 

many missing lines 


3. MAMS/VAS Retrieval Products 

Mesoscale atmospheric variability was investigated by making 100 
meter resolution atmospheric temperature and moisture profiles using 
MAMS and geostationary satellite borne VISSR Atmospheric Sounder (VAS) 
data. MAMS and VAS data collected during COHMEX (June 19, 1986) were 
combined in a retrieval algorithm that allowed the high resolution MAMS 
data to be substituted for similar VAS spectral channels. VAS CO 2 
sounding channels not contained on MAMS provided most of the vertical 
resolution in the profiles, while the MAMS data, with its 100 meter 




resolution, provided the horizontal resolution. MAMS/VAS derived image 
products of skin temperature, precipitable water (atmospheric moisture), 
and lifted index (atmospheric stability) were made from the retrievals. 
The imagery showed several interesting features. In particular, a 
tongue of low level moist air pushing into the COHMEX region from the 
south was captured in the precipitable water product. This feature was 
also observed in the COHMEX radiosonde network, however not with the 
definition provided by the 100 meter resolution MAMS/VAS moisture 
product. Interestingly, an area of low level cumulus formed in the 
moist tongue region within a few hours of the imagery. The MAMS/VAS 
lifted index product showed a bias towards unstable values when compared 
to the COHMEX radiosondes, as has been seen previously with VAS lifted 
index products. The skin temperature product showed a definite heat 
island in the urban land-use region of Nashville, TN. 

Meso f) (20-200 km) and meso 7 (2-20 km) scale gradients in the 
MAMS/VAS moisture, stability, and surface skin temperature products were 
evaluated using a structure function analysis. Among the significant 
findings of this analysis, the meso 7 scale moisture and skin 
temperature gradients clearly exceeded the strength of the Meso scale 
gradients. This is important information for investigations of 
mesoscale meteorological events. Meso 7 scale stability gradients were 
marginally stronger than the Meso scale gradients. However, this may 
only have occurred because the VAS data (with it # s much coarser 8 km 
horizontal resolution) influences the stability product more than the 
other products. Also, for all the products, detectable gradient 
information was found at resolutions down to 300 meters, the smallest 
resolution at which the gradients were evaluated. These results, while 


from only a single case study, indicate that significant atmospheric 
variability exists at very small scales. The importance of this 
variability to small scale meteorological phenomena has not been 
thoroughly assessed. Further investigations into mesoscale atmospheric 
gradients are warranted by these findings. 

A comparison of the MAMS/VAS moisture gradients was made to 
gradients of VAS-only atmospheric moisture data. It was found that the 
VAS-only gradients were significantly weaker, especially for meso 7 and 
smaller meso /? scales (<50 km) . This indicates that significant 
atmospheric moisture gradient information is being missed in operational 
practice by the VAS instrument and suggests a need for an improved 
resolution operational sounding instrument in geostationary orbit. A 
paper has been submitted to the Int. Journal of Remote Sensing on this 
topic; a copy is included in the appendix of this report. 

4. MAMS/VAS Fluxes 

Atmospheric radiative fluxes were studied by producing 100 meter 
resolution MAMS/VAS net fluxes at 1000 mb, 700 mb, and 300 mb. 

Regression coefficients from previously developed relationships between 
VAS data and simulated radiation fluxes were applied, using VAS data 
with MAMS spectral channels substituted for similar VAS spectral 
channels as predictors. A repeat overpass of MAMS data collected over 
COHMEX on June 19, 1986 was used. Results showed that 1000 mb and 700 
mb net fluxes were highly related demonstrating a strong influence of 
surface temperature on the net flux at these levels. The 300 mb net 
flux was smoother than the 700 mb or 1000 mb net fluxes, demonstrating 
the effect of atmospheric absorption. Cooling rates from the net flux 


data showed very little variation in the 1000 - 700 mb layer, reflecting 
the high correlation between the 1000 mb and 700 mb net fluxes. Cooling 
rates for the 700 - 300 mb layer were more variant. Also, midday net 
fluxes at 1000 mb _and 700 mb decreased with time in the eastern third of 
the flight track. This agreed with the location of a moist tongue being 
advected into the COHMEX region (described above) and also with the 
development of the low level cumulus field. Net fluxes at 1000 mb and 
700 mb increased with time in the western two- thirds of the flight 
track. Simulations showed that a decrease of net flux at a level can 
result from increased moisture below the level; an increase in net flux 
can be caused by decreasing moisture below the level. 

A shortcoming of the radiation net flux effort was its inability 
to produce reasonable fluxes over small low level cumulus clouds. This 
occurred because the 8 km resolution of VAS channel 5, which is 
susceptible to low level cloudiness, was insufficient to resolve 
individual cumulus clouds. Other VAS channels affected by the low level 
clouds were replaceable by MAMS data. This points out the importance of 
high resolution data for quantitative analysis in and around cloudy 
fields of view. 

The net flux products (and cooling rates) in this effort combine 
information available separately in satellite temperature and moisture 
products. The net flux product however represents a different way of 
packaging this information which may be useful to radiative budget 
studies of the earth-atmosphere system and modelling efforts which 
attempt to assimilate radiative data. Using the radiometric data 
directly in a regression algorithm*, instead of producing atmospheric 


temperature and moisture profiles from the radiometric data and then 
computing fluxes, saved considerable computation time. 

5. MAMS/VAS Cloud Clearing Activities 

MAMS COHMEX data was used to investigate the effect of low level 
clouds on VAS retrievals. MAMS data from June 15, 1986 was used. On 
this day the COHMEX region was covered by an extensive low level 
scattered to broken cumulus cloud field. MAMS 11/im and 12/im "cloud 
cleared" radiances were used to replace cloud contaminated VAS 11/im and 
12/im radiances. MAMS data was "cloud cleared" by applying brightness 
temperature thresholds to the MAMS data, which can see between the 
clouds with its high resolution. The non-cloudy MAMS fields of view 
were then averaged to simulate the VAS field of view. MAMS averaged 
non-cloudy field of view temperatures were 1 - 3 K warmer than the VAS 
cloud contaminated data. Cloud contaminated retrievals were compared to 
"cloud cleared" retrievals. The impact of cloud clearing the 11/im and 
12/im data was most evident in the moisture profile. Temperature profile 
modification was generally negligible. The moisture profile showed a 
general 5% increase in moisture in the lowest 300 - 500 mb (about 1 - 
2mm of total column precipitable water). While the response of the 
moisture profiles to the cloud clearing is encouraging, a lack of 
vertical resolving power in the lower atmosphere in the retrievals 
limited the potential impact. Temperature and moisture profile impact 
would likely be increased if the VAS 13.3/im data (VAS channel 5) were 
cloud cleared by a similar MAMS channel (MAMS currently does not have a 
13.3/im channel). This channel is an important contributor to the low 
level temperature profile but is sensitive to low level cloud 


contamination. The impact of unresolved cloud on sounding data from 
instruments carrying more lower atmospheric sounding channels than VAS 
may be significantly greater, increasing the need to use high spatial 
resolution data for atmospheric sounding. 
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Radiometric calibration of IR Fourier transform spectrometers: 
solution to a problem with the High-Resolution Interferometer 
Sounder 


Henry E. Revercomb, H. Buijs, Hugh B. Howell, D. D. LaPorte, William L. Smith, and L. A. Sromovsky 


A calibrated Fourier transform spectrometer, known as the High-Resolution Interferometer Sounder (HIS), 
has been flown on the NASA U-2 research aircraft to measure the infrared emission spectrum of the earth* 
The primary use-atmospheric temperature and humidity sounding — requires high” radiometric "precision 
and accuracy « of the order of 0.1 and I°C, respectively). To meet these requirements, the HIS instrument 
performs inflight radiometric calibration, using observations of hot and coid blackbody reference sources as 
the basis for two-point calibrations at each wavenumber. Initially, laboratory tests revealed a calibration 
problem with brightness temperature errors as large as 15°C between 600 and 900 cm' 1 . The symptom of the 
problem, which occurred in one of the three spectral bands of HIS, was a source-dependent phase response. 
Minor changes to the calibration equations completely eliminated the anomalous errors. The new analysis 
properlv accounts for the situation in which the phase response for radiance from the instrument itself differs 
from that for radiance from an external source. The mechanism responsible for the dual phase response of 
the HIS instrument is identified as emission from the interferometer beam splitter. 


I. Introduction 

The capability of measuring absolute radiance with 
an infrared Fourier transform spectrometer is impor- 
tant to the High-Resolution Interferometer Sounder 
(HIS) program, 1 - 2 which is applying interferometry to 
measure the upwelling emission spectrum of the earth 
for retrieving the temperature, humidity, and other 
parameters of the atmosphere and surface. The pri- 
mary objective of the HIS program is to improve signif- 
icantly the vertical resolution of temperature and hu- 
midity profiles determined from satellite platforms by 
increasing the spectral resolution of observed ra- 
diances. The interferometer is a natural choice for the 
job. Interferometers have proved their merit in space 
with successful measurements of the emission spectra 
of planetary atmospheres, not only of the earth, ^ but 
also of Venus, 3 Mars, 6 Jupiter, Saturn and Uranus. 7 
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The multiplex and throughput advantages of the in- 
terferometer 8 make it possible to make radiometrically 
precise observations at a much higher spectral resolu- 
tion than that of current filter radiometers (tens of 
cm“ ! ). To improve the vertical resolution in the tro- 
posphere and lower stratosphere by at least a factor of 
2 over that of current temperature profiling radiome- 
ters (3-8 km), a resolving power (X/AA) of M000 is 
needed in the spectral range from 3.7 to 17 gin. (The 
15-/itm CO 2 absorption band is the primary band for 
temperature sounding, the combined 4.3-^m CCU and 
4,5-^m NoO bands are used to enhance the tempera- 
ture resolution in the lower troposphere, and the 6.7- 
Mm H 2 O absorption band gives humidity soundings.) 
High radiometric precision is required because radio- 
metric noise and time-dependent and wavelength-de- 
pendent calibration errors are magnified in the inver- 
sion process to derive atmospheric parameters. To 
obtain temperature profiles with rms errors of <1°C 
from high-resolution measurements requires noise- 
equivalent-temperature errors and calibration repro- 
ducibilities of the order of 0.1°C and absolute errors of 
less than about 1°C. 

Several important steps have been taken toward 
developing an improved temperature and humidity 
sounding satellite instrument. The feasibility of 
building an instrument to achieve the radiometric per- 
formance needed from a geosynchronous satellite plat- 
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ig. . unctionaJ schematic of HIS optics. The primary’ mirror, collimating mirror, and focusing mirrors are shown as lenses to simplify the 
rawing. ane re ecting surfaces are shown as dashed lines, and the two paths of the interferometer are functionally represented. The 
complete instrument is -2.7 m (9 ft) long and fits into the 45.7-cm (18-in.) diam wing pods of the NASA U-2 aircraft. 


form has been demonstrated. 9 An aircraft model HIS 
has been developed and flown on dozens of flights of 
the NASA U-2 research aircraft as a direct demonstra- 
tion of the scientific value of such an instrument 10 and 
that the radiometric noise and calibration require- 
ments are achievable. 11 The instrument was flown 
with many other atmospheric sensing instruments in 
two major NASA field programs, the combined Hunts- 
ville Meteorological Experiment (COHMEX) for 
studying severe storms and the First ISCCP Regional 
Experiment (FIRE) for studying the effect of cirrus 
clouds on climate. The unique ability of this instru- 
ment to measure accurately the emission spectrum 
from a flexible high-altitude platform with a large 
complement of other instrumentation should make it 
an important resource for many types of experiment 
for several years. Currently, plans for modifying the 
newest NOAA geosynchronous sounding instrument 
to incorporate an interferometer are being investigat- 
ed to realize improvements in the satellite sounding 
data for weather forecasting before the next century. 

This paper describes information about the radio- 
metric calibration of the HIS aircraft instrument, 
which may be useful to other IR interferometer appli- 
cations. Early laboratory testing of the aircraft in- 
strument revealed that wavelength-dependent cali- 
bration errors of up to 15° C occurred between 600 and 
900 cm 1 when blackbodies at 300 and 77 K were used 
to determine the radiance of a 280 K blackbody from 
the magnitude spectra of each source (Sec. II). These 
large errors, which were originally believed to require 
hardware changes to the instrument, have been com- 
pletely eliminated by the processing technique de 
scribed in Sec. III. The instrument characteristic re- 
sponsible for the initial calibration errors is also 
identified there. The implications of the new calibra- 
tion analysis technique for interferometer design and 
performance are summarized in the final section. 


n. HIS Radiometric Calibration and the Problem 


A. Calibration Approach and Laboratory Test Results 

The basic approach for determining absolute ra- 
diances from the HIS nadir-viewing interferometer is 
the same as that used for filter radiometers and has 
been used successfully for other interferometric appli- 
cations. 3-7 ' 12 The detectors and electronics are de- 
signed to yield an output which is linear in the incident 
radiance for all wavenumbers in the optical passband 
of the instrument, and two blackbody reference 
sources are viewed to determine the slope and offset 
which defi ne the linear instrument response at each 
wavenumber. 

In the HIS U-2 instrument, calibration observations 
of the two onboard reference blackbodies are made 
every 2 min. There are four double-sided optical-path 
scans of each reference source for every twelve scans of 
the earth. As shown in Fig. 1, which summarizes the 
optical configuration, the blackbodies are viewed by 
rotating the telescope field of view (FOV) from below 
the aircraft to inside a blackbody aperture using a 45° 
plane mirror. There are no uncalibrated optical sur- 
faces, since the earth is viewed through an open aper- 
ture in the pod, which provides an aerodynamic shell. 

The small size of the optical beam at the blackbody 
positions makes the design of accurate radiation stan- 
dards relatively easy. The reference blackbodies are 
thermoelectrically controlled blackened copper cavi- 
ties. The insulated copper walls of the blackbody 
cavities give good temperature uniformity, and be- 
cause of the cavity effect, the normal emissivity is very 
close to one (Table I). The temperatures are sensed 
with accurately calibrated platinum resistance ther- 
mometers (PRTs) embedded in the base of each cavity. 
(During testing, a second PRT in the side of the cavity 
was used to verify adequate temperature uniformity.) 


Tabl« 1. Characteristics of ths HIS Aircraft Instrumsnt 

Spectral range (cm' 1 ) 0 


Band I 

590-1070 

Band II 

1040-1930 

Band III 

2070-2750 

Field of view diameter imrad) 


Telescope 

100 

Interferometer 

30 

Blackbody reference sources 


Emissivity 

>0.998 

Aperture diameter (cm> 

1.5 

Temperature stability (Ki 

±0.1 

Temperatures <K) 

240, 300 

Autoaligned interferometer: Modified BOMEM BBDA2.I 

Beam splitter 


Substrate 

KC1 

Coatings P/ 4 X at 3.3 Mm) 

Ge + Sb^Si 

Maximum delay (double-sided current configuration (cm) 

Band I (hardware limit is ±2.0) 

±1.8 

Bands II and III (limited by data system) 

4*1.2, -0.8 

Michelson mirror optical scan rate (cm/s) 

0.6- 1.0 

Aperture stop (at interferometer exit window) 


Diameter (cm) 

4.1 

Central obscuration area fraction 

0.17 

Area (cm*) 

10.8 

Area-solid angle product (cm- sr) 

0.0076 

Detectors 


Type 

Ar-doped Si 

Diameter (cm) 

0.16 

Temperature (K) 

6 

Nominal instrument temperature tK) 

260 

1 The ranges shown are design ranges. The current bandpass 
filters were chosen from available stock filters and will be changed as 

new filters are acquired. 


One important additional requirement when apply- 
ing a two-point calibration with blackbody references 

to an interferometer, as opposed to an 

instrument 

measuring spectra directly, is that the instrument re- 
sponsivity should be independent of optical delay (or 

that any delay dependences should be 

accurately 

known). Avoiding sources of delay-dependent re- 
sponse was a major objective in designing the HIS 
instrument. To accomplish this, care was taken in the 
optical stop design and alignment to prevent the effec- 
tive aperture stop size from changing with the motion 


of the Michelson mirror. The best location for the 
aperture stop, which is focused on the detectors, was 
found to be at the exit window of the interferometer 
module (see Fig. 1). Furthermore, the FOV of the 
interferometer is restricted to 30 mrad to limit self- 
apodization. 

Now we turn to the mathematical expressions for the 
calibration. First, a formalism which leads to the use 
of magnitude spectra in the expression for calibrated 
radiance is presented to show where this commonly 
used approach can create a problem, as occurred with 
the early (1985) HIS calibrations. Assuming linearity 
as expressed above, the output interferogram i^can be 


expressed in terms of the incident spectral radiance L v 
as follows, using a continuous representation: 

F(x) * 2 I exp[i0fi')) exp(t'2jr*'x)di', (1) 

where the uncalibrated magnitude spectrum (C r s C -, ) 
is given by 

C ( . * 1/1 = r XL, + L'X (2) 

and where* = optical path difference (delay), 
v = wavenumber, 

<p(v) = phase response of instrument [4>{v) = 
(Pi-*)], 

r v = responsivity of instrument, 

L ° = offset from instrument emission, re- 
ferred to input, 

~ = complex Fourier transform. 

The phase characterizes the combined optical and 
electrical dispersion of the instrument and here is as- 
sumed to be the same for scene and background emis- 
sions. Although this assumption is commonly made, 
it will later be shown to be invalid for the HIS instru- 
ment. Equation (2) expresses the linear relationship 
between the uncalibrated spectrum and spectral radi- 
ance. The two unknowns to be determined from the 
two calibration observations are the responsivity and 
the offset radiance. The offset radiance defined here 
is the radiance, which, if introduced at the input of the 
instrument, would give the same contribution as the 
actual emission from various parts of the optical train. 
Equation (2) written for both the hot and cold black- 
body views can be solved to yield 

r^iC he -Cj/[B v (T h )-B,(T c )), (3) 

L° = CJr t . - 5,(7^) = C Jr v - (4) 

where B v is the Planck blackbody radiance, and sub- 
scripts h and c label the quantities associated with the 
hot and cold blackbody. [Note that for simplicity the 
blackbodies are assumed to have unit emittance here. 
To account for actual emittances t, the Planck ra- 
diances should be replaced with eB + (1 - t)B(T a ) t 
where T a is the ambient temperature.] Solving Eq. (2) 
for the source radiance and substituting from Eqs. (3) 
and (4) yield the basic calibration relationships: 

L v = CJr , - C; (5) - 

L t . * [<C, - CJ/<C*„ - CJ][B v (T h ) - B V (T C ) ] + B,(T e ). (6) 

The ground calibration tests to be reported here 
consisted of measuring the radiance from a blackbody 
at ^280 K using calibration blackbodies at 300 and 77 
K. The uncalibrated spectra are shown for two of the 
three HIS spectral bands in Figs. 2 and 3 (see Table I 
for the nominal spectral coverage of each ba nd). 
These bands were chosen because band I displays the 
problem to be discussed here, and band II does not. 
<Band III also does not.) Both the magnitude and 
phase are determined directly from a complex Fourier 
transformation of the measured two-sided interfero- 
gram. As will be explained further in Sec. Ill, the 
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600 700 800 900 1000 1100 

WAVENUMBER Ccm' T ) 

Fig. 2. Uncalibrated magnitude (a) and phase (b) spectra of biack- 
bodv sources for HIS spectral band I. The curves are labeled by the 
temperature of the source. The phase labeled 300K-77K Differ- 
ence is the phase of the difference of the 300 and the 77 K spectra. 
The magnitude spectra are shaped by a Gaussian numerical filter 
and display CO? absorption features and channeling as discussed in 
the text. Note the substantial differences among the phase responses. 

linear phase contribution arising from an ambiguity in 
the sample offset from the zero path difference (ZPD) 

has also been removed. 

The magnitude spectra have various features which 
need explanation. Because testing was conducted in 
air t the instrument transmittance is significantly af- 

— fected by the 1-2 m of air in the path from the blackbo- 
dies to the detectors. (Of the total optical path be- 
tween the blackbodies and the detector/Dewar 
assembly, only the interferometer itself is enclosed and 
backfilled with dry nitrogen to protect the beam split- 
ter; at flight altitudes of 20 km, atmospheric absorp- 
tion is insignificant.) Therefore, the magnitude spec- 

trumfor band I (600-1 100 cm -1 ) shows CO? absorption 

between 600 and 750 cm -1 , and band II (1100-1800 
cm" 1 ) shows water vapor absorption beyond 1300 
cm" 1 . The absorption for the 668-cm" 1 COo line and 

— for several H 2 O lines is so strong in air that the signal is 
almost zero, and the phase is poorly defined. The 
general Gaussian shape of the magnitude spectra is 
caused by the numerical filtering which is performed in 
the instrument digital electronics. (A hardware con- 


1100 1200 1300 1400 1500 1600 1700 1800 

WAVENUMBER (cm' 1 ) 

Fig. 3. Uncalibrated magnitude (a) and phase (b) spectra of black- 
body sources for HIS spectral band II. The curves are labeled as in 
Fig. 2. Here also the magnitude spectra are shaped by a numerical 
filter and display channeling as discussed in the text. The very deep 
and numerous lines in the magnitude spectra are due to H^O absorp- 
tion. In contrast to band I, the phase spectra are very similar and 
quite linear. 

volution is performed for signal-to-noise preserving 
sample volume reduction by factors of 14, 8, and 8 in 
the three spectral bands.) 

The sinusoidal components superimposed on the 
magnitude spectra are channeled spectra caused by 
the parallel surfaces of the arsenic-doped silicon detec- 
tors. (The band III detector has an antireflective 
coating, and the spectra do not display channeling.) 
The .channeled spectra are very stable because the 
detectors, which are operated in a liquid He Dewar, 
have a very small coefficient of expansion and experi- 
ence only small temperature variations. As will be 
shown in the next section, channeling does not affect 
calibrated spectra. 

The phase spectra for bands I and II differ markedly. 
For band II the phases are nearly linear, the behavior 
expected with an ideal beam splitter having zero dis- 
persion and with an electrical response having a pure 
time delay. Band I phases, on the other hand, show 
significant deviations from linearity and also vary as 
the source radiance is varied. The nonlinearity is 
optical, not electrical, in origin because the sign of the 






deviations from linearity depends on delay scan direc- 
tion. The phase spectra shown are for one scan direc- 
tion. The corresponding phases for the other direc- 
tion are approximately a mirror image about z rad. 

B. Problem 

It is probably not surprising, in light of the anoma- 
lous band I phases, that direct application of the nor- 
mal calibration procedure (based on magnitude spec- 
tra) to the laboratory test data does not work 
uniformly well. The band I spectrum derived in this 
way is shown in Fig. 4. The spectrum is presented as a 
brightness temperature to make any errors stand out 
as a deviation from the measured blackbody tempera- 
ture of 280.2 K. Errors of tens of degrees occur, and 
the similarity of the wavelength dependence of the 
errors to the nonlinearities of the phase is apparent. 
Figure 5 shows that the calibrated spectra for the HIS 
bands with nearly linear phase spectra are extremely 
accurate when calibrated in the same way. The region 
of large variance from 280.2 K beyond 1400 cm' 1 for 
band II and centered near 2350 cm' 1 for band III is 
caused by water vapor and CO 2 in the optical path and 
does not occur when the instrument is operated at 
flight altitude. Apparently, the normal calibration 
approach cannot handle measurements with phase 
anomalies like those of band I. 

III. Solution 

A. New Calibration Analysis 

The band I phase spectra shown in Fig. 2(b) indicate 
that the phase becomes more nearly linear as the 
source radiance becomes larger. This behavior sug- 
gests that the total measured interferogram and the 
corresponding complex spectrum have two compo- 
nents; a normal phase component for radiance from 
the source and an anomalous phase component for 
radiance from the instrument itself. Of course, the 
radiance from the instrument can originate from many 
individual components with different phase character- 
istics, but they can be lumped together into one term 
with one phase. 

The new calibration analysis is a generalization 
which properly accounts for a dual phase response, if it 
is present. It is really a minor modification of the 
procedure presented in Sec. II. The difference is that 
the complete complex spectra obtained from Fourier 
transforming the measured interferograms are used 
for calibration, not just the magnitudes. Let the com- 
plex uncalibrated spectrum be represented by 

C ¥ = Py - r \L , + L* 1 exp[i< 2 > (*')]! expfiXv)], (7) 

where 0 0 is the difference from the normal phase of the 
anomalous phase associated with the combined radi- 
ance from the many emitting components of t he in - 
strument. Then it is clear, under the same instrument 
stability assumptions on which the basic calibration 
approach depends, that the anomalous phase contri- 
bution can be eliminated along with the instrument 
radiance offset by differencing complex spectra from 



Fig. 4. Brightness temperature spectra for band I. The spectrum 
with the large deviations from the measured 280.2 K temperature of 
t he blackbody source is derived from the original calibration analysis 
using magnitude spectra. The spectrum which accurately repro- 
duces a constant brightness temperature is from the modified analy- 
sis presented here. The large variance between 600 and 650 cm " 1 is 
caused by the low instrument transmittance in that region. 
Changes have since been made to improve the throughput in this 
region. 
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WAVENUMBER (cm' 1 ) 

Fig. 5. Brightness temperature spectra for band II fa) and band III 
<b). The calibration for these bands is good and essentially identical 
tor the standard and modified analyses. The larger variance beyond 
1400 cm’ 1 for band II and centered at 2350 cm* 1 for band III are due 
to H 2 O and CO 2 in the path of the interferometer. 


different sources. The difference spectra are identical 
to the difference spectra which would result if there 
were no anomalous phase contribution. The equa- 
tions for the difference spectra from Eq. (7) are 

C - C = rjlr, - B V {T C )\ exp(i0(W|, (8) 

C, p - C, * r\B v {T h ) - B t .(T c )] exp[*0<*)]. ( 9 ) 

The new expression for the responsivity, which follows 
immediately from Eq. (9) by taking the magnitude of 
both sides, is 

r , * \C hw - C cw \l[B ¥ {T h ) - B,(T f )]. ( 10 ) 

Compared to Eq. (3), the difference of magnitude spec- 
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tra is replaced by the magnitude of difference spectra. 
The offset which follows directly by substituting the 
responsivity into Eq. (7) is 

L?exp[iX*')] * C ht . exp[-i>M]/r 4 . - B£T^. (11) 

Note that while the instrument responsivity remains a 
real function in this formalism, the offset is now com- 
plex. Finally, the basic calibration expression which 
follows by taking the ratio of Eq. (8) to (9) is 

L, * R e[(C„ - Cj/(C h , - C f ,)][fi (7\) - B t (T c )] + B r (T v l (12) 

For ideal spectra with no noise, this expression for the 
calibrated radiance would be real, since the phases of 
the ratioed difference spectra are the same. This can- 
cellation of the phases avoids the square root of two 
noise amplification, which can be associated with tak- 
ing the magnitude of spectra with nonzero phase. Be- 
cause the phase of the ratio of difference spectra is zero 
to within the noise, the calibrated spectrum can equal- 
ly well be defined in terms of the real part of the ratio 
(as shown) or in terms of the magnitude of the ratio. 

This technique solves the problem. Application to 
the measurements presented in the last section yields 
an accurately calibrated spectrum for band I. The 
new brightness temperature spectrum, compared to 
the old in Fig. 4, shows no sign of the large dip centered 
at 740 era" 1 . Also, the phase associated with radiance 
from the source (the phase of difference spectra) is 
found to be reasonably linear as shown in Fig. 2. The 
impact on the band II and III spectra for a 280 K 
blackbody is negligible. However, the small phase 
dependence on source for band II would probably 
cause a detectable error for colder sources. 

It should be pointed out that ambiguities in the 
phase from one spectrum to the next must be eliminat- 
ed before applying this technique. For the HIS instru- 
ment. which does not use a white light source to main- 
tain a fixed offset between ZPD and the delay of a 
numerically filtered point, a discrete ambiguity occurs. 
This translational ambiguity is linear in wavenumber 
and takes the form 

Ot - rmk/MR , (13) 

where k is an integer, m is the spectral sample number, 
M is the number of points in the spectrum, and R is the 
numerical filter sample reduction factor. The ambi- 
guity is removed by determining the k for each spec- 
trum, which nearly eliminates the phase slope. Ambi- 
guities of 7r rad, which occur in instruments with 
ambient temperature detectors, are not usually a con- 
sideration for the HIS application which uses LHe 
cooled detectors. 

Also, note that if a white light source is used to 
control numerical filtering so that the same delays are 
measured on subsequent scans, the differencing to 
eliminate instrument background can be performed on 
the interferograms. Then anomalous phases would 
have no effect on spectra. Differencing interfero- 
grams is also the way to handle a potential dual phase 
problem, if one-sided OPD scanning is used. 


B. Source of the Anomalous Phase Spectrum 

It is of interest to explore the cause of the anomalous 
band I phase to assess the applicability of the new 
analysis technique to other FTIR applications. Also, 
it may be possible to relax constraints on interferome- 
ter optical design, since some apparent hardware prob- 
lems can be eliminated by analysis. 

Phase nonlinearities can of course occur from dis- 
persion in beam splitter/compensator substrates, if 
their thicknesses are not well matched. This effect 
cannot explain the HIS data. The linearity of the 
phases for bands II and III suggests that the beam 
splitter/compensator matching is quite good. More 
important, dispersion from this origin does not give 
phase spectra which depend on the source. 

To explain the HIS data, a mechanism is needed for 
which radiance from the instrument can yield a phase 
spectrum which is different from that for an external 
source. Apparently, the plane of wavefront division in 
the beam splitter depends on the source of the radia- 
tion. We have identified two mechanisms by which a 
source dependence might occur: (1) dependence of 

the beam splitter coating properties on the angle of 
incidence and (2) emission from the beam splitter coat- 
ings. 

The first mechanism, angle dependence, has been 
ruled out as an important factor for the HIS configura- 
tion. A test was performed to measure off-axis radia- 
tion, It consisted of viewing a LNo source with the 
instrument in its normal configuration and with the 
field stop removed. The spectrum with the field stop 
removed is dominated by radiance from the Dewar 
window entering the interferometer exit window, not 
radiance from the field stop as in the normal configura- 
tion. The difference in the spectra for these two con- 
figurations is the spectrum of radiation from the field 
stop. The phase of the difference is essentially identi- 
cal to that of radiance from an external source (Fig. 2), 
indicating no significant dependence on the angle of 
incidence. 

The second possible mechanism, beam splitter coat- 
ing emission, appears to be a likely candidate for ex- 
plaining the anomalous phase. Beam splitter emis- 
sion could lead to the effective wavefront division 
occurring at the point of emission with coherent radia- 
tion being emitted into both legs of the interferometer. 
Emission at a different depth in the beam splitter than 
normal wavefront division would create an anomalous 
phase. 

There is evidence that the HIS beam splitter has 
absorption in the wavenumber region where phases are 
anomalous. Although the absorption of the beam 
splitter used for the measurements presented here has 
not been measured directly, transmittance measure- 
ments for other beam splitters of the same construc- 
tion are available. Figure 6 shows the transmittance 
of one such beam splitter compared to a sinusoidal fit 
used to extrapolate the expected transmittance in the 
absence of absorption into the region of the narrow 
absorption feature centered at —740 cm -1 . The mini- 
mum reflectance at 3000 cm" 1 is consistent with the 
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Fig. 6. Transmittance of a beam splitter of the same construction 
as the HIS instrument beam splitter. The absorption feature cen- 
tered at ~740 cm' 1 is responsible for the duai phase response of the 
instrument. The sinusoidal fit to the transmittance is used to 
approximate the reflectance for beam splitter emittance and effi- 
ciency calculations. 

thickness of the beam splitter coatings, 0.25 wave- 
length at 3.3 Mm of Ge and the same of Sb 2 S 3 . While 
germanium has a weak absorption feature centered 
near 740 cm -1 , the absorptance of up to 0.1 or 0.2 
implied by these data is surprising, because the thick- 
ness of the coatings is not large enough for substantial 
bulk absorption. Maybe the indicated absorption is 
caused by a surface effect or contamination. 

It is possible to obtain information on the beam 
splitter emittance and efficiency directly from mea- 
surements of the HIS instrument itself, if a model of 
the beam splitter is assumed. The two simple models 
illustrated in Fig. 7 have been used for this purpose. 
Rays illustrating the reflection process for an external 
source and the emission process are shown for both. 
Model 1 represents the beam splitter/compensator as a 
single plane with nonzero reflectance p and emittance 
f. It is too simple to explain a phase anomaly, because 
there is a single surface for wavefront division due to 
both reflection and emission. Model 2 uses two plane 
surfaces, one which can reflect but not emit and one 
which can emit but not reflect, A beam splitter repre- 
sentable by this model would create an anomalous 
phase. 

Emittances and efficiencies are de r -rmined from the 
ratio of the uncalibrated magnituue spectrum for 
emission from the beam splitter (referred to interfer- 
ometer input) to that of a blackbody at the same tem- 
perature. The blackbody spectrum, free of emission 
contributions, is obtained in the normal manner by 
differencing the complex spectrum of an ambient tem- 
perature blackbody from the spectrum of a liquid ni- 
trogen blackbody. The beam splitter spectrum is ob- 
tained by subtracting a scaied blackbody spectrum 
from the liquid nitrogen spectrum with the scale factor 
chosen to give approximately zero from 950 to 1100 
cm' 1 where absorption is expected to be small. The 
ratio defined in this way is equal to the ratio of the 
expressions for the emitted and reflected output given 





MODEL 1: ONE PLANE SURFACE 



MODEL 2: TWO PLANE SURFACES ^ 
Fig. 7. Simplified beam splitter models for emittance and efficien- 
cy calculations. The expressions in terms of the emittance t and 
reflectance p represent the amplitudes of the beam at various loca- 
tions. The rays on the left represent the emission process, and those 
on the right represent the passage of an external beam. 


in Fig. 7 (e.g., ratio = e/p for model 1). Using the 
reflectivity from the sinusoidal fit of Fig. 6 [fit = (1 — 
p)], the expression for the ratio can be solved for the 
emittance. The efficiency of the beam splitter is then 
given by the product of the complete beam splitter 
reflectivity and transmittance divided by the ideal 
output amplitude of 0.25 for a reflectance of 0.5 and no 
absorptance. 

The emittances and efficiencies deduced from inter- 
ferometer measurements using these models are 
shown in Fig. 8. The emittances estimates are in 
reasonable agreement with those obtained from the 
direct transmittance data of Fig. 6, giving further sup- 
port for the conclusion that there is beam splitter 
emission acting approximately as modeled. Notice 


3216 APPLIED OPTICS / Voi. 27. No. 15 / 1 August 1988 





Fig. 8. Beam splitter emittance and efficiency estimates labeled by 
the model assumed. Model 2' is the same as model 2 with the 
emitting and reflecting surfaces reversed in order. The dashed 
emittance curves are from the data of Fig. 6, and the solid curves are 
from measurements with the HIS instrument. 


that due to the multiple passes through the absorbing 
layer, the efficiency can get as low as 50% at 740 cm“h 

Further verification of the influence of beam splitter 
absorption on the HIS instrument is provided by com- 
parison of the measured instrument responsivity to 
that calculated based on the transmittance of each 
optical element and on the detector responsivities, as 
shown in Fig. 9. The Gaussian numerical filter func- 
tion has been divided out of the measured responsivity 
(note the effects of COo absorption and channeling 
discussed earlier). The responsivity calculated as- 
suming the beam splitter to be uniformly efficient is 
noticeably too large between 650 and 900 cm* 1 where 
the band I phase anomaly occurs. The beam splitter 
efficiencies from Fig. 8 reduce the responsivity over the 
correct wavenumber region, with the model 2 calcula- 
tion giving very good agreement with the measured 
responsivity. 

Both the deduced emittances and the efficiencies 
give a strong case for the emission explanation of 
anomalous phases. However, a piece of puzzle that 
does not fit very well is the size of the phase anomalies. 
The optical thickness (x = 5<j>/27rv) corresponding to 
the observed phase anomalies (5<p) at 740 cm -1 is 3.9 
Atm. Model 2 suggests that, if the absorption is in the 
coatings, the additional optical path difference should 
be <1.9 /im, the optical thickness of the combined 
beam splitter coatings divided by the cosine of the 30° 
angle of incidence. This peculiarity and the mecha- 
nism for the absorption are remaining mysteries. 

IV. Implications and Summary 

It is not known whether the problem identified here 
is common or rare in other FTIR applications. The 
errors resulting from this problem could possibly be 
small enough to go unnoticed but might be significant 
when evaluating detailed performance. Comparing 
the phase spectra for a cold and a warm source is a 
simple test for diagnosing the problem. 



Fig. 9. Comparison of instrument responsivity calculated from 
optical component transmittances and detector responsivity to th 
end-to-end responsivity from calibration measurements. The beam 
splitter efficiencies used in the calculated responsivities are those 
deduced from measurements with the HIS interferometer itself. 
The calculations using model 2 seem to account properly for actual 
beam splitter efficiencies. 


The new calibration analysis technique has been 
discussed here in the context of the HIS application, 
involving absolute radiance measurements. It may 
also be important for applications in which the trans- 
mittance or reflectance of a sample is measured. As 
pointed out by Tanner and McCall, 13 emission from 
the samples or the surroundings can create significant 
errors for these types of measurement. However, the 
radiance contribution from the background can be 
eliminated, if the spectra for two different reference 
sources are differenced. By determining the differ- 
ence spectra using complete complex spectra as de- 
scribed in the last section, any anomalous phase contri- 
butions of radiance from the interferometer will also 
be eliminated. 

There are some aspects of interferometer design that 
mav be influenced by the analysis approach presented 
here and by knowledge that beam splitter emission can 
modify the phase response. First, the constraints on 
beam splitter coating materials can probably be re- 
laxed to include materials with some absorption. Sec- 
ond. interferometer designs may be able to rely on a 
better capability to handle errors associated with the 
phase. It is conceivable that having a well-behaved 
phase is not so important when a proper calibration is 
performed with double-sided scanning. 

In summary, the phase response of one spectral band 
of the HIS FTIR instrument appears to vary with the 
source radiance, becoming more linear the higher the 
radiance contribution from the source. This depen- 
dence creates substantial radiance errors when magni- 
tude spectra are used to perform a two-point calibra- 
tion using blackbody reference sources. The errors 
are eliminated, yielding good calibration results, when 
a minor modification to this technique is applied. The 
new technique explicitly accounts for the possibility 
that a dual phase response can occur. That is, the 
phase response for radiance from the instrument can 



differ from that for radiance from a source. The origin 
of the different phase response for radiance from the 
HIS instrument is identified as emission from the 
beam splitter. 
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team at the Space Science and Engineering Center, 
BOMEM, Inc., the Santa Barbara Research Center,, 
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Murcray for many helpful discussions. The HIS pro- 
gram is jointly funded by NOAA contract NA-84- 
DGC-00095 and NASA contract NAS5-27608. 

References 

1. W. L. Smith, H. E. Revercomb, H. B. Howell, and H. \1. Woolf, 
“Recent Advances in Satellite Remote Sounding/’ Internation- 
al Radiation Symposium 84: Current Problems in At mo- 
spheric Radiation, G. Fiocco. Ed. (A. Deepak, Hampton, VA, 
1984), p. 388. 

2. W. L. Smith, H. E. Revercomb, H. B. Howell, and H. M. Woolf. 
“HIS — A Satellite Instrument to Observe Temperature and 
Moisture Profiles with High Vertical Resolution,” in Fifth Con- 
ference on Atmospheric Radiation (American Meteorological 
Society, Boston, 1983). 

3. R. A* Hanel B. Schlachman, F. D. Clark, C. H. Prokesh, J. B. 
Taylor, W. M. Wilson, and L. Chaney, ‘The Nimbus III Michel- 
son Interferometer,” Appl. Opt. 9, 1767 (1970). 


4. R. A. Hanei, B. Schlachman, D. Rodgers, and D. Vanous, “Nim- 
bus 4 Micheison Interferometer,” Appl. Opt. 10, 1376 (1971). 

5. D. Oertei et al. , “Infrared Spectrometry of Venus from Venera- 
15 and Venera-16,” Adv. Space Res. 5, 25 (1985). 

6. R. A. Hanel et al., “Mariner 9 Micheison Interferometer,” Appl. 
Opt, 11,2625(1972). 

7. R. A. Hanel et al., “Infrared Spectrometer for Voyager,” Appl. 
Opt. 19, 1391 (1980). 

8. J. W. Brault, “Fourier Transform Spectroscopy,” High Resolu- 
tion Astronomy, Proceedings, Fifteenth Advanced Course in 
Astronomy and Astrophysics, Saas-Fee, M. Huber, A. Benz, and 
M. Mayor, Eds. (1985). 

9. “A Design Feasibility Study for the High-Resolution Interfer- 
ometer Sounder (HIS),” Santa Barbara Center Final Report for 
contract UAA 871R55 5 (10 July 1981, updated 19 July 1982, 
updated 15 Feb. 1983). 

10. W. L. Smith, H. E. Revercomb, H. M. Woolf, H. B. Howell D. D. 
LaPorte, and K. Rageyama, “Improved Geostationary Satellite 
Soundings for the Mesoscale Weather Analysis/Forecast Opera- 
tions,” in Proceedings , Symposium on Mesoscale Analysis and 
Forecasting , Vancouver, Canada, 17-19 Aug. 1987, ESA SP-282 
(1987). 

11. H. E. Revercomb, D. D. LaPorte, W. L. Smith, H. Buijs, D. G. 
Murcrav, F. J. Murcray, and L. A. Sromovsky, “High-Altitude 
Aircraft Measurements of Upwelling IR Radiance: Prelude to 
FTIR from Geosynchronous Satellite,” Mikrochim. Acta in 
press, 000 (Springer-Veriag, Wien, 1987). 

12. D. D. LaPorte and R. Howitt, “Ambient Temperature Absolute 
Radiometry using Fourier Transform Spectrometers,” Proc. 
Soc. Photo-Opt. Instrum. Eng. 364 (1982). 

13. D. B. Tanner and R. P. McCall “Source of a Problem with 
Fourier Transform Spectroscopy," Appl Opt. 23, 2363 (1984). 


APPLIED OPTICS / Voi. 27. No. 15 / 1 August 1988 


3218 



SPECTROSCOPIC INFERENCES FROM HIS MEASUREMENTS 
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1225 West Dayton Street 
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H. M. Woolf, and H. B. Howell 
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1, _ INTRODUCTION 

Radiometrically accurate observations of the earth's emission spectrum 
from 3.8 to 16.6 microns have been made using the High- resolution 
Interferometer Sounder (HIS) to look downward from the NASA U2/ER2 aircraft or 
upward from the ground. These observations have been used to demonstrate the 
substantially improved vertical resolution of temperature and water vapor 
soundings derived from high resolution spectra (resolving power from 1800 to 
3800), as compared to soundings from the low resolution filter radiometer 
observations used in current satellite sounders ^ ^ 

The HIS observations have also demonstrated that Fourier Transform - 
Infrared (FTIR) instruments are especially well suited to absolute emission 
measurements of broad spectral bands at high resolution .^’ ^ ° A fundamental 
advantage of FTIR instruments for accurate calibration is wavelength 
integrity, the same property which has made FTIR the standard for very high 
resolution absorption measurements. 

While the HIS was primarily designed for atmospheric temperature and 
water vapor profiling, its broad spectral coverage and high accuracy (absolute 
brightness temperature errors of less than 0.5®C and a reproducibility 
approaching 0.1°C at most wavelengths) lead to a wide range of important 
applications, including monitoring tropospheric trace gases, surface and cloud 
spectral properties, and surface temperatures, in addition to spectroscopy . 

2 . SPECTROSCOPY 

The long wavelength part of a HIS downwelling radiance spectrum is 
compared to a calculated spectrum in Fig. 1. The calculation uses the AFGL 
HITRAN/86 line file and FASC0D2 line -by-line program with atmospheric state 
data from in situ measurements. In general, agreement between HIS and 
FASC0D2 spectra is remarkably good, a tribute to the current state of spectral 
line files and line-by-line codes. 

The detailed radiance differences of observations from FASC0D2 
calculations are shown in Figs. 2 and 3 for a range of atmospheres. Note that 
the largest differences are reproducible over a wide range of conditions, 
including even uplooking and downlooking viewing conditions. 

The top two differences in Figs. 2 and 3 are from uplooking observations 
(note the reversed sign) . The in situ atmospheric temperature and humidity 
measurements for these observations are very good because NCAR Cross-chained 
Loran Atmospheric Sounding System (CLASS) balloons were launched from the same 
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site where the HIS was located* The last comparison labelled Gulf was flown 
at 400 mb over the Gulf of Mexico aboard the NOAA P3 aircraft during an 
integration test flight and is somewhat noisy* It is included because, again 
in this situation, the atmospheric measurements made by sensors onboard the 
plane and by NCAR dropsondes are very good and because it is an example of a 
reasonably wet atmosphere (the wetness of the atmospheres shown generally 
increases from the top of the figure to the bottom) . 

The largest residuals between observations and calculations have been 
identified with either spectroscopic oriatmospheric uncertainties. Two 
examples are labelled in Fig. 2* The effects of line mixing on the effective 
line shapes in the CO 2 Q-branches was not accounted for in FASC0D2. The 
latest version (FASC0D3), soon to be available from AFGL, does account for 
line mixing. The effects of absorption by chlorof luorocarbons (CFC's), not 
currently included in FASCODE, are also labelled. These features are shown 
in more detail in Fig. 4, along with absorption by carbon tetrachloride. 

The major new conclusions which have been drawn from these comparisons 
are illustrated in Figs. 4 and 5. Figure 4 shows that only a 10% adjustment 
of the H 2 O amount is needed to bring the. calculated and observed spectra into 
excellent agreement regarding the longwave H 2 O continuum. This adjustment is 
on the order of the uncertainty in the H 2 O profile. Over most of the spectral 
range, the residual differences are small after accounting for this 
adjustment, the larger effects discussed above, and small H 2 O line parameter 
uncertainties. However, a major difference still remains from about 700 to 
760 cm . The correlation of this region with the absorption of CO 2 and the 
symmetric effect on the other side of the CO 2 band (Fig. 2), strongly suggests 
that the line absorption of CO 2 needs to be increased in the model. 

Figure 5 shows that the positive discrepancy between 1250 and 1425 cm^ ; 
which is largest for the driest atmospheres, is related to the H 2 O continuum, 
because the continuum contribution (the difference between a FASC0D2 
calculation without and with a continuum) accounts for the shape of the 
positive discrepancy. The continuum in this region is dominated by the 
foreign-broadened continuum, unlike the longwave region where the self- 
broadened continuum dominates. Therefore, it is possible to reduce the 
continuum for this region by the large amount needed (60%) , without destroying 
the agreement already discussed In the longwave window. 

Finally, the consistent differences at the H 2 O lines from 1100 to 1230 
cm implies the need to increase the line strengths of water vapor in this 
region by about 30%. It is important to note that even large adjustments of 
the atmospheric water vapor amount cannot eliminate both this problem and the 
continuum problem simultaneously, because they require changes in the opposite 
sense . 

3. CONCLUSIONS 

Reproducible differences between HIS observations and FASC0D2 line-by- 
line calculations lead to the following conclusions: 

(1) The FASC0D2 water vapor continuum in the longwave window region 

from 10 to 13 microns (750 to 1000 cm’ ) gives reasonable agreement 
with radiance observations. 


The model H 2 O continuum from 7 to 8 microns (1250 to 1425 cm'^) 
needs adjustment to reduce its contribution by about 60%. 

CO 2 absorption in the region from 13.1 to 14.3 microns (700 to 760 
cm"'*') is too small in the model. 

Water vapor line strengths in the region from 8.1 to 9.1 microns 
(1100 to 1230 cm‘^) need to be increased about 30%. 
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VERTICAL RESOLUTION AND ACCURACY OF ATMOSPHERIC INFRARED 

SOUNDING SPECTROMETERS 

Hung-Lung Huang, William L. Smith and Harold M. Woolf 

Cooperative Institute for Meteorological Satellite Studies 
University of Wisconsin, Madison 

ABSTRACT 

A theoretical analysis is performed to evaluate the 
accuracy and vertical resolution of atmospheric profiles 
obtained with the HIRS/2, GOES I/M and HIS instruments. In 
addition, a linear simultaneous retrieval algorithm is used 
with aircraft observations to validate the theoretical 
predictions. Both theoretical and observational results 
clearly indicate that the accuracy and vertical resolution 
of the retrieval profile would be improved by high soectral 
resolution and broad spectral coverage of infrared radiance 
measurements . 

The HIS is found to possess the equivalent of 10 pieces 
eacn or temperature and water vapor independent precise 
measurements. The characteristics for temperature include a 
ver tical resolution of 1-6 km with an accuracy of IK, and 
for water vapor a vertical resolution of 0.5-3 km with an 
accuracy of 2K in dewpoint temperature. The HIS is a factor 
of almost five times better in vertical resolution and a 
factor of two times better in accuracy than the GOES I/M and 
HIRS/2 filter radiometers. 
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1. Introduction 

Radiance measurements by satellite borne multispectral 
infrared radiometers are used to sense atmospheric 
temperature and moisture profiles (Kaplan, 1959; Wark and 
Fleming, 1966; Smith and Woolf, 1976; Susskind et al., 
1984) . Because the radiances arise from very thick layers of 
the atmosphere, the poor vertical resolution and associated 
accuracy of the derived profiles limit their use in 
Numerical Weather Prediction (NWP) (Bengtsson 1979; Phillips 
et al., 1979; Smith 1991). The presence of clouds which 
absorb infrared radiation further complicates the derivation 
of the sounding profiles from the radiance observations 
(Chahine, 1974/77/82) and also limits their utility in NWP. 

The poor vertical resolution of current operational 
satellite sounding instruments is attributable, in part, to 
low spectral resolution. Due to this coarse spectral 
resolution, the radiometer does not sense radiation between 
individual absorption lines so that the emitted radiance 
reaching the satellite is a mixture of radiation from the 
high atmosphere, due to the strong absorption near line 
centers, and radiation from the lower atmosphere, from 
between the absorption lines. Thus, the relatively poor 
spectral resolution of current instruments causes a vertical 
smearing of the atmospheric structure. 

To minimize vertical smearing and thereby approach the 
ultimate vertical sounding resolution and accuracy 
achievable with passive measurements, a sounding instrument 
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must possess a spectral resolving power (defined by x/sx, 
where x is wavelength) of at least 1000. For example, a 
spectral resolution of 0.7 cm -1 is required in the 600-700 
cm -1 (15 jun) thermal emission band of CO2. This spectral 
resolution can be achieved using a Michelson interferometer 
(Smith et al., 1979/83/90) or a large detector array grating 
spectrometer (Chahine et. al, 1990). The High-resolution 
Interferometer Sounder (HIS) and Advanced InfraRed Sounder 
(AIRS) are experimental instruments designed to demonstrate 
improved sounding performance from high spectral resolution 
and continuous broad spectral coverage. An aircraft model of 
HIS is a Michelson interferometer with a spectral resolving 
power of approximately 2000 , covering a spectral range of 
3.7-16.7 fim (Smith et al . , 1983). The AIRS is to achieve 

measurement characteristics similar to HIS by incorporating 
new cooled fecal plane detector array technology into a 
grating spectrometer (Chahine et al.,1990). 

In this paper, the linear simultaneous retrieval 
algorithm (Smith et al., 1987/91) along with the 
theoretical vertical resolution (Backus and Gilbert, 1970, 
Conrath, 1972, Thompson, 1982, and Rodgers, 1987) and error 
analysis (Rodgers, 1987 and Smith et al., 1991) are used to 
quantify the vertical resolution and accuracy provided by 
three different sounding instruments: (1) the GOES I/M 

sounding radiometer; (2) the HIRS/2 aboard the NOAA 
sa ^®Hites; (3) the HIS sounding instrument proposed for 

future GOES implementation. 


4 


2 . Data 

HIS data have been successfully collected from NASA 
U2/ER2 aircraft flights (Smith et al., 1987; Revercomb et 
al., 1988). With observations at a six second sampling 
interval from the height of 65,000 ft, an instantaneous 
ground resolution of 2 km is achieved. Figure 1 shows a 
clear sky spectrum of brightness temperature observed during 
the cooperative Huntsville Meteorological Exp eriment 
(COHMEX) , with a spectral resolution of about 0.5 cm -1 from 
600-1100 cm ^ (9.1-16.7 pm) , and 1.0 cm ^ resolution from 

1100-2700 cm -1 (3.7-9 . 1 pm) (Huang, 1989). Major absorbing 
constituents are labeled accordingly, namely, CO 2 , H 2 O, O 3 , 
N 2 O, CH 4 and CO. The selected weighting functions for the 
three HIS infrared bands are presented in Figures 2 and 3 . 
These curves express the vertical derivative of atmospheric 
transmittance with respect to the natural logarithm of 
pressure and are calculated from a line-by-line model 
FASCODE II (Clough et al., 1936) using FASCODE standard 
atmosphere profile conditions. As for the GOES I/M and 
HIRS/2 data, they are simulated from the HIS radiance 
spectra by convolution with their spectral response 
functions. 

3. Retrieval Algorithm, Vertical Resolution and Error 

Analysis Methods 

a. Linear Simultaneous Retrieval Algorithm 
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A computationally efficient retrieval methodology has 
been developed by Smith et al., (1991), which is a unique 
linear algorithm that simultaneously retrieves temperature 
and absorbing constituent profiles (i.e. water vapor, ozone, 
methane et al.) from observations of spectral radiances. The 
Radiative Transfer Equation (RTE) linearization results from 
a definition for the deviation of the true gas concentration 
profiles from an initial state in terms of the deviation of 
their "effective" temperature profiles from the true 
atmospheric temperature profile. The "effective" temperature 
profile for any absorbing constituent is defined as that 
temperature profile which satisfies the observed radiance 
spectra under the assumption that the initial absorber 
concentration profile is correct. Differences between the 
effective temperature, derived for each absorbing 

constituent, and the true atmospheric temperature are 
proportional to the error in the initial state of the gas 
concentration profiles. The gas concentration profiles are 
thus determined after inversion of the linearized RTE from 
the retrieved effective temperature profiles assuming that 
one of the assumed concentration profiles is known (i.e. 
C0 2 ) . 

Following the algorithm of Smith et al., 1991, the 
linearized RTE is 

N p s 

5T B^ i/(Ps) r y(Ps) 5T S “ - f (P) <STj_ (p) r 0 ^ (p) dlnr 0 ^ (p) 



where 0° v (p) =[3B„ (T°)/3T]/ [3B„ (T° B )/3T] , 

N=the number of optically active atmospheric constituents, 
N=2 when only the uniformly mixed gases ( mainly CO2) 
and water vapor are considered, 
p s =pressure with subscript s denoting the surface, 
r j/(P) =t he total transmittance of the atmospheric above 
atmospheric pressure level p, 

5( ) =the difference between the true quantity and the 
initial state denoted by a superscript 0, 

E( ) i=the summation of i quantities, 
r l/ j_=the transmittance of the atmosphere for the ith 
absorbing constituent, 

Tg„=the brightness temperature at spectral frequency u 
T(p)=the true atmospheric temperature profile, 

T i(P) =t he effective temperature of the ith absorbing 
constituent, 

u i(P)=P ath length profile of the ith absorbing constituent 


Ui(p)=l/g 


P 

qi(p')dp', 

Jo 


g= acceleration due to gravity, 

<5i (P) ^mixing ratio of the ith absorbing constituent, 
5T i (p)=T i (p)-T°(p) . 

In matrix form, again, following Smith et al., 1991, the 
generalized statistical/physical solution of equation (1) is 
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t r =Ct b (2a) 

and the retrieval coefficient matrix C is 

C=(A*E~ 1 A+S~ 1 ) _i A*E _1 (2b) 

where subscript r denotes retrieved quantity. The elements 
of the weighting function matrix A are 
A y ,o=/3°^(Ps) r °^(Ps) 

K, j=-3 0 l /(Pij)’-°y(Pij)dlnr 0 I/i (p i . . j) , j=l,2, . . . ,N*M i; 


where N is the number of constituents and Mj_ is the number 
of quadrature pressure levels denoted to each constituent. S 
is the sample statistical covariance matrix and E is the 
covariance of the brightness temperature error ( ) * denotes 
the matrix transpose, and ( ) _1 denotes the matrix inverse. 

b. Analysis of Retrieval Vertical Resolution 

The vertical resolution of the various sounding 
instruments is analyzed using the concepts of Backus and 
Gilbert (1970), Conrath (1972), Thompson (1982) and Rodgers 
(1987). The retrieved temperature can be expressed as a 
smooth version of the true temperature profile, 

T r (Pi)= S R(Pi,?j)T(pj) (3) 

j 

where Tj-(pj_) is the retrieved temperature at pressure level 
i, T(pj) is the true temperature profile, and R(Pi,Pj) is 
the vertical resolution function or averaging kernel for 
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level pj^. (3) describes the retrieved temperature at any 
level Pi as a vertical average of the true temperature 
profile weighted by the vertical resolution function R for 
the level pj_. In terms of the sounding retrieval algorithm 
given by (2) t r =Ctfc, and the RTE (1), written in matrix 
notation as tfc,=At, then one can rewrite t r as 
t r =CAt 
=Rt. 

with 

R=(A*E" 1 A+S -1 ) _1 A*E" 1 A (4) 

It has been found that the concept of resolution 
is vague and has many possible definitions (Thompson, 1982 
and Rodgers, 1987). 3ackus and Gilbert (1970) defined the 
"spread" as a measure of the resolution and used to quantify 
the vertical resolution of a retrieval system using 
satellite radiance measurements. However, Newman (1979) and 
Thompson (1982) have pointed out that the definition of 
spread suffers from several mathematical deficiencies which 
are caused by the oscillatory sidelobes of resolution 
functions that results in an unreliable estimations of the 
resolution. The mathematical difficulties become increasing 
severe with increasing spectral resolution and measurement 
precision. The computations of spread of HIS, HIRS/2 and 
GOES I/M (not shown) also confirm these findings. 

It is these difficulties that motivate the development 
of an alternative technique to define an effective vertical 
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resolution, which would directly apply to any satellite 
retrieval system. The appendix gives an abbreviated 
derivation of this tec hn ique (Purser, personal 
communication) which will be discussed in complete detail in 
a subsequent article. If one lets J be the trace of matrix 
R, than J can be thought as the total effective degrees of 
constraint (independent precise measurement) imposed by the 
satellite data. The measure of the local effective data 
density is 

Pi= 2 F i, j R j , j 

3 

where 

F i,j= R2 j,i/( 2 R 2 j,k5Zk)- 
k 

5 Z]£ is the height increment at the level k. And the vertical 
resolution defined in terms of data density is simply, 

Wi-1/Pi (5) 

c. Analysis of Retrieval Error 

Within the framework of the linear simultaneous 
inversion theory, a retrieval error analysis algorithm is 
developed following the concepts of Rodgers (1987) . 

The covariance of the retrieval temperature error G, 
G=< ( t-t r ) (t-t r )*>, can be obtained from 

t r =ct b =C(t' b +e) ~ c ( AT+e) =CAt+Ce , 
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where angle bracket denotes expectation operator, and e is 
assumed to be the brightness temperature measurement error, 
as 

G= ( t-CAt-Ce) (t-CAt-Ce) * 

= [ ( I-CA) t-Ce ] [ ( I-CA) t-Ce ] * 

where I is the identity matrix. By assuming temperature, t, 
and random brightness temperature error, e, are 
uncorrelated, then 

G= ( I-CA) tt*( I-CA) *+Cee*C* 

= (I-CA) S (I-CA) *-C EC*=V+M. 


Here 


V= (I-CA) S (I-CA) * 

(6) 

M=CEC* 

(7) 


and S=tt* and ^=ee*. V and M can be defined as the vertical 
resolution component error and the measurement noise 
component error, respectively (Rodgers, 1987). The total 
root mean square (ms) error is the square root of the 
diagonal elements of matrix G, if ignoring inter-level 
correlations, and can be written as 

( 8 ) 

(8) is obtained by using the matrix identity 


G=(A*E- 1 A-S - 1 )- 1 
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I- (X+Y) -1 X= (X+Y) -1 Y . 

The vertical resolution error results from those 
components of the profile that cannot be recovered by the 
retrieval process. In other words, the accuracy of the 
retrieved profiles is limited by the vertical resolution of 
observations. This vertical resolution related error is the 
square root of the diagonal elements of the matrix V 
(Rodgers, 1987). The contribution to the retrieval error due 
to radiance measurement error is given by the square root of 
the diagonal elements of the matrix M (Rodgers, 1987) . It 
should be noted that the systematic retrieval error due to 
the imperfection of the linearized RTE and inverse model is 
not included in G. 

4. Results 

a. Results of Retrieval Vertical Resolution and Error 
Analysis 

Retrieval vertical resolution and error are determined 
for HIS, GOES I/M and HI RS /2 sounding instruments. A 
radiance error equivalent to a brightness temperature 
measurement error (or noise) of 0.25K at a scene temperature 
of 2 60K is assumed for all measurement situations. The 
sample covariance matrix S is derived from 400 mid-latitude 
global climatological atmospheric profiles (Smith et al., 
1974). The global sample mean profile and FASCODE II (Clough 
et al . , 1986) are used to calculate the spectrum of 

atmospheric transmittance profiles for all three instruments 
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considered and used for all vertical resolution and error 
analysis computations. The uniformly mixed constituents 
(CO 2 , N 2 O, CH 4 , SO 2 and CO) are combined to form a single 
transmittance profile so that together with water vapor the 
number of atmospheric transmittance components was two 
(ozone absorption regions are excluded in this analysis) . 
Forty pressure levels between 0.1 hpa to 1000 hpa are 
considered for each of the two constituents. 

Results of the retrieval vertical resolution functions 
R for temperature and water vapor are shown in figure 4 and 
5, respectively. For temperature vertical resolution, one 
can see little difference in these functions (except near 
the surface) between HIRS/2 and GOES I/M sounding 
conditions. In general, both HIRS/2 and GOES I/M sounders 
only are able to resolve very broad vertical layers of the 
temperature profile. However, examination of HIS temperature 
vertical resolution functions (figure 4c) reveals that HIS 
possesses many more fine layers of temperature profile 
information. Similar conclusions can be made for water vapor 
vertical resolution functions of HIRS/2, GOES I/M and HIS. 
Nevertheless, the GOES I/M sounder seems to be able to 
obtain a little more vertical water vapor information than 
the HIRS/ 2 , because of its additional tropospheric moisture 
channel (7.43 „m) . 

Table 1 presents the total effective degrees of 
constraint (precise independent measurements) for 
temperature and water vapor retrievals. It is found that 
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both HIRS/2 and GOES I/M are only able to obtain 
approximately two pieces each of temperature and water vapor 
retrievals. In particular, GOES I/M seems to possess less 
temperature information than HIRS/2 (1.8 vs. 2.2). In 
contrast to temperature, GOES I/M possesses more water vapor 
information than HIRS/2 (2.7 vs. 1.3). The total effective 
independent measurement possesses by HIS, however, has ten 
for each of temperature and water vapor. This indicates a 
factor of almost five improvement of the HIS over that of 
GOES I/M and HIRS/2. 

Vertical resolution functions of 300 hpa temperature 
and 700 hpa water vapor together with their vertical 
resolution W obtained from (5) are shown in figure 6. The 
lengths of the fine vertical bars are the vertical 
resolution W. The examples presented in figure 6 demonstrate 
the reliable estimation of the vertical resolution which 
their function has oscillatory sidelcbes. Figure 7 shows the 
vertical resolution computed from (5) which is defined in 
terms of data density and is in the unit of km. It can be 
seen that a major improvement in the vertical resolution of 
temperature and water vapor retrieval profiles is 
demonstrated by the HIS instrument approach. 

Figure 3 shows the total rms error of HIS and GOES I/M 
temperature and water vapor concentration profiles. This 
result suggests that, consistent with the vertical 
resolution analysis shown earlier, sounding accuracy is 
improved by a tactor of 2-3 over that of current radiometers 
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by using measurements with high spectral resolution and 
broad spectral coverage. 

Figure 9 demonstrates the resulting total random rms 


error 

of 

temperature 

and 

water 

vapor profiles 

for different 

spectral 

resolutions 

of 

HIS 

measurements. 

These 

three 

spectral 

resolutions 

are 

the 

full HIS, half 

HIS and one 

tenth 

of 

HIS spectral 

resolution, each with the same 

noise 

(0.25 

k) 

and spectral 

coverage (600-2700 

cm ■*■) . 

The 


sensitivity of the retrieval error to specrral resolution is 
large, with the higher spectral resolution providing the 
most accurate sounding product. 

Different "optimal" sets of HIS channels have been 
selected to investigate retrieval error sensitivity to the 
total number of spectral channels used. For each "optimal" 
case, those channels whose weighting functions have the 
greatest sensitivity to each of the forty sounding levels 
between 0.1 to 10C 2 hpa were selected for long wavelength 
(15 A.m) , mid wavelength (6.2 M m) and short wavelength (4.3 
^m) bands, respectively (Huang, 1989). Figure 10 presents 
rms error of temperature for different sets of "optimal" 
channels. Curve A in this figure is the result obtained by 
using 19 GOES I/M c.nannels and is presented as a reference. 
Curves B to G are cased on optimally selected HIS channel 
sets of 31, 104, 222, 604, 2181 and 6144 (complete set) . As 

shown in this figure, the addition of more "optimally 
chosen" channels toes provide more independent retrieval 
information and results in a smaller retrieval error. 
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b. Results Of Synthetic and Aircraft Observation Retrievals 
Figure 11 shows the rms error of sounding retrievals 
from HIS radiance spectra synthesized from an independent 
set of 38 special radiosondes during COHMEX. The temperature 
and dew point soundings were retrieved using the solution of 
(2) based on a dependent statistical set of 153 different 
COHMEX radiosondes. The guess profile used for retrieval is 
the mean of these 153 soundings. Figure 12 shows an estimate 
of the rms error of retrievals from actual HIS measurements 
from the NASA ER-2 aircraft for the two COHMEX days, 15 and 
19 June 198 6. The results in Figure 12 were obtained by 
Bradshaw and Fuelberg (1990) using statistical structure 
function analysis. As can be seen by comparing figures 11 
and 12 for the COHMEX region, observational results confirm 
the theoretically expected accuracies close to 1 K for 
temperature and 3 K for dewpoint temperature. 

3. Conclusions 

The theoretical retrieval vertical resolution and 
profile accuracy of GOES I/M, HIRS/2 and HIS are quantified 
and compared by using a linear simultaneous retrieval 
algorithm and vertical resolution and error analysis 
•techniques. These observed and simulated retrievals indicate 
that high spectral resolution and broad spectral coverage 
are indeed the key to both the improved vertical resolution 
and the profile accuracy that results. 
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In conclusion, the passive remote atmospheric sounding 
of HIS is found to possess the equivalent of 10 pieces each 
of temperature and water vapor independent precise 
measurements. The characteristics for temperature include a 
vertical resolution cf 1-6 km with an accuracy of IK, and 
for water vapor a vertical resolution of 0.5-3 km with an 
accuracy of 3K in dewpoint temperature. The HIS is a factor 
of almost five times better in vertical resolution and a 
factor of two times better in accuracy than the GOES I/M and 
HIRS/2 filter radiometers. 
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APPENDIX 

A Theory of Effective Resolution 

Let vertical resolution functions matrix R be defined 
as in (4) , and from the fact that the precision matrices 
(A*E~~A) and S - * are nonnegative symmetric matrices it 
follows that each eigenvalue \ of R is bounded within the 
interval [0,1] and since the sum of nonvanishing eigenvalues 
is simply the trace of the matrix, 

J-Tr (R) (Al) 

it can be found that 


0 < J < M 


(A2) 


17 


where M is the number of independent measurements. The case 
J=0 occurs in the singular limit of vanishing measurement 
precision, E -1 -> 0, and implies no effect of the satellite 
measurements on the retrieval, while the case J=M occurs 
only in the case of infinite measurement precision, E -> 0, 
and implies a retrieval that always fits all measurements 
exactly. Thus it seems at least intuitively reasonable to 
consider J to be a measure of the total effective degrees of 
constraint imposed by the satellite data used in the 
retrieval . 

It would be useful to be able to extrapolate from the 
measure J of the effective total quantity of satellite data 
to a measure of the local effective retrieval density. The 

simplest choice that integrates to J would be, 

/ ^ 

p i =R i, i/ > z i (A3) 
where 6Zi is the height increment of the level i. In theory 
there is no guarantee that seme comoonents Ri ,• are not 
negative and in practice the distribution p is invariably 
noisy. However, i - is observed (Purser, personal 
commun ica t ion ) that, averaged locally over a scale 
comparable with that implied by the structure of the 
corresponding rows and columns of R, the typical magnitudes 
of p genuinely reflect the local density of distinct 
measurements when these measurements a re prec i se and are 
themselves localized. A serviceable measure of local 
effective data density p is therefore given by applying a 
unimodular nonnegative smoothing operator F to /, 


that is, 
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Pi~ s . F i,j R j,j 

1 



(A4) 

with 




Z Fi,jSZi=l 
i 



(A5) 

For the univariate data one need look no 

further 

than R 

itself to find the material from which 

to 

construct a 

smoother with the desired qualities. Of 

the 

two 

natural 

definitions , 




i, j=R 2 i, j/ I R^k,j 5Z k 
k 



(A6) 

F (b) i,j=R 2 j,i/ 2 R 2 j, k *Z k 
k 



(A7 ) 


the first is tore universal in_ the sense that it continues 
to apply when the retrieval formalism is generalized to an 
"optimal interpolation" (Daley, 1991) that incorporate 
point-like measurements (such as radiosonde temperature 
observations) , but the second method has been found to yield 
superior estimates when, as in the present study, attention 
is confined exclusively to remotely sensed data whose 
weighting functions (the rows of A) are mainly much smoother 
-han the covariance profiles (rows of S) . It is the second 
form, ( A7 ) , that is therefore adopted here. 
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Table Legends 


Table 1: The total effective degrees of constraint 

(precise independent measurements) for temperature 
and water vapor retrievals obtained from GOES I/M, 
HIRS/2 and HIS. 


TEMPERATURE 


WATER VAPOR 


GOES I/M HIRS/2 HIS 


1.3 2.2 10.1 


2.7 2.1 10.3 


iable 1, the total effective degrees of constraint (precise 
independent measurements ) for temperature and water vapor 
retrievals obtained from GOES I/M, HIRS/2 and HIS. 


Figure Legends 


Figure 1: Brightness temperature spectra measured by HIS 
from U2 aircraft over Huntsille, Alabama on 15 
JUNE 1986 Of COHMEX. 

Figure 2: Temperature component Planck radiance weighting 
functions for a small set of HIS band 1 and 
band 3 channels. 


Figure 3: Water vapor component Planck radiance weighting 
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(a) Water vapor vertical 
of GOES I/M sounder. 
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(b) Water vapor 
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(c) Water vapor 
of HIS. 
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Vertical resolut 

ion functi 
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resolution (the fine vertical bars) of HIS 
300 hpa temperature and 700 hpa water vapor 
retrievals . 


Figure 7: (a) Temperature vertical resolution (km) of HIS 

and GOES I/M sounder. 

Figure 7: (b) Water vapor vertical resolution (km) of HIS 

and GOES I/M sounder. 

Figure 3 : RMS temperature and water vapor concentration 
retrieval errors of KIS and GOES I/M sounder. 

Figure 9: Temperature and water vapor concentration 
RMS errors of HIS in full resolution, half of 
the HIS resolution, and one tenth of the HIS 
resolution. 

Figure 10: Temperature RMS errors for the use of 
different optional selections of HIS channels 
and GOES I/M sounding channels. 


Figure 11: Simulated temperature and dew point 

temperature RMS errors and standard deviation 
of statistical profile set from the mean 
profile. 

Figure 12 : Temperature and dew point temperature RMS 
error of ER-2 HIS retrievals during COHMEX (15 
& 19, JUNE 1986) . 
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ABSTRACT 

A high spectral resolution interferometer sounder (GHIS) has been designed for flight on future geosuuo nary 
meteorological satellites. It incorporates the measurement principles of an aircraft prototype instrument, which 
has demonstrated the capability to observe the eanh-emuted radiance spectrum with high accuracy. The aircraft 
results indicate that the theoretical expectation of 1*C temperature and 2°-3 8 C dewpoint retrieval accuracy 
will be achieved. The vertical resolution of the water vapor profile appears good enough to enable moisture 
tracking in numerous vertical layers thereby providing wind profile information as well as thermodynamic 
profiles of temperature and water vapor. 


I. Introduction 

Because of the temperature and water vapor profiling 
capabilities needed to meet the 1 990s initiative for im- 
proved mcsoscaie weather prediction, the remarkable 
performance history of Michclson interferometry in 
space (Hand et al. 1967, 1971, 1972, 19S0; Schafer et 
ah 1989), and the recent successful aircraft demon- 
stration of the High-rcsolution Interferometer Sounder 
(HIS) technology (Revercomb et ah 1988a. 1988b, 
1988c; Smith et ah 1988), the time is right to proceed 
with plans to place interferometer sounders on geo- 
synchronous satellites. The objective of this paper is 
to show the performance improvement of the GHIS 
(GOES High-resolution Interferometer Sounder ) over 
the current satellite sounding instruments, the char- 
acteristics of the instrument that emerged from a Phase 
A study, and the results achieved from tests of the air- 
craft prototype interferometer sounding instrument. 

Higher vertical resolution temperature and moisture 
soundings from geosynchronous orbit are needed to 
make full use of the improved wind observations 
planned for the 1990s. For mesoscale weather predic- 
tion, it is important to have observations of tempera- 
ture, humidity, and wind distributions that together 
provide information on atmospheric stability, and 
horizontal and vertical motion fields. Improved wind 
measurements will soon be available to enhance 
weather prediction capability in the United States, since 
NOAA is installing a demonstration network of 30 
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Doppier radar wind profilers, spaced 500 km apart, 
across the midwestem United States. The Doppler ra- 
dars provide wind profiles on an hourly basis with an 
accuracy of I m s' 1 and a vertical resolution of I km. 
However, the profiler will not have the spatial coverage 
needed for mesoscale or regional models. GOES ob- 
servations arc needed to provide coverage on a finer 
spatial grid, '•'•’-i- 

To allow atmospheric stability to be determined 
from GOES-acquired temperature and water vapor 
measurements with an accuracy consistent with profiler 
wind accuracy, ft can be shown ( using the thermal wind 
relationship) that the temperature and water vapor 
distribution must be observed with a vertical resolution 
of 1-2 km and an accuracy of 1°C for temperature 
and 3°C for dewpoint. The GOES I/M filter wheel 
instrument will fall far short of meeting these require- 
ments. as will be demonstrated in the following section. 
Experience with GOES-VAS data confirms the fact 
that geostationary satellite temperature profile data 
with the poor spectral (vertical) resolution provided 
by filter wheel radiometers offer little information on 
mesoscale phenomena. The primary uses of the VAS 
instrument for mesoscale applications have been in the 
definition of the gross water vapor structure from the 
sounding retrieval process and the diagnosis of the 
horizontal wind field by tracking cloud and upper-level 
water vapor motion. 

Aside from improving thermodynamic observations 
from satellites, improved vertical resolution water va- 
por profiles are needed to supplement wind profile ob- 
servations over coastal oceanic regions. From geosta- 
tionary satellites, the GHIS improved water vapor pro- 
files should enable the vertical profile of horizontal 
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wind to be achieved by tracking the water vapor motion 
observed within 1-2 km vertical layers. 

In summary, greatly improved temperature and wa- 
ter vapor profile resolution and accuracy from the 
GOES sounding instrument are needed to meet the 
mesoscale meteorological data requirements of this 
decade and beyond. As will be shown, the GHIS can 
provide the improved performance. 

2. Performance improvement 

To obtain the vertical temperature and moisture 
profile resolution required, the sounding instrument 
must achieve near continuous spectral coverage 
throughout the 600-2600 cm’ 1 region with a spectral 
resolution of 0.1% (AX/X) (Smith et al. 1979, 1983a). 
The high spectral resolution is needed to avoid smear- 
ing the upw'elling radiance contributions from relatively 
opaque absorption lines with the contributions irom 
more transparent regions in between the absorption 
lines. For example, in the thermal emission bands of 
C0 2 , a spectral resolution of 0.7 cm" 1 is needed in the 
600-700 cm’ 1 ( 15 ) region and 2cm" 1 in the 2300— 

2400 cm’ 1 (4.3 jxm) region. Although the required 
spectral coverage and spectral resolution are beyond 
the capabilities of filter wheel radiometers, they can be 
achieved using an interferometer. 

The HIS aircraft instrument flown aboard the NASA 
U2 and ER2 aircraft has proven the performance ca- 
pabilities of a high spectral resolution interferometer. 
Ground calibration tests and airborne science missions 
have demonstrated the ability of the HIS to measure 
radiometrically accurate emission spectra with a res- 
olution far exceeding the performance of contemporary 
filter radiometers (Revercomb et al. I988a,b.c). Ra- 
diance spectra have been obtained during the more 
than 40 aircraft missions with a resolution of 0.35 cm’ 1 


in the 600 to 1 1 00 /im region, as compared to the typ- 
ical 15,0 cm -1 resolution of filter radiometers, such as 
the one to be flown on GOES I/M. 

Figure I shows a typical spectrum of infrared radia- 
tion brightness temperature sensed t the HIS. The 
ban dwidth s of the filters planned for me GOES I/M 
sounding radiometer are superimposed. As can be seen, 
the filter wheel radiometer severely smears the fine- 
scale spectral radiance structure of the atmosphere. The 
spectral smearing causes unwanted absorption con- 
tamination in atmospheric “windows” used for sensing 
the earth's surface temperature, and it greatly limits 
the vertical resolution of temperature and w-atcr vapor 
profiles because it broadens the atmospheric weighting 
functions ( Fig. 2). 

Figure 3a shows more clearly the effect of spectral 
resolution on radiance smearing. It can be seen from 
these brightness temperature spectra observed with the 
HIS and portrayed at three different spectral resolutions 
that the 15 cm’ 1 resolution of the GOES I/M filter 
wheel instrument causes an extreme smearing ol energy 
emission from the lower atmospheric levels (brightness 
temperatures in excess of 240 K) with energy emission 
from upper atmospheric levels (brightness tempera- 
tures below 240 K). 

The vertical resolution of atmospheric temperature 
profiles as sensed by the GOES I/M filter wheel radi- 
ometer is compared to that of the interferometer 
sounder in Fig. 3b. (The vertical resolution was deter- 
mined using the procedure described by Rodgers 1987.) 
As can be seen, a factor of 2-3 improvement is achieved 
by the GHIS. Figure 4 shows results of a retrieval error 
analysis conducted using the “linear simultaneous re- 
trieval method” and associated error diagnosis scheme 
described by Smith et al. (1990), and the expected 
sounding accuracy of the GOES I/M filter instrument 
compared to the GHIS. Simulations also demonstrating 
these substantial improvements in accuracy have been 


FiG. 2. Comparison of a small selection of GOES Interferometer (GHIS) spectral channel Planck radiance prohle weighting 
functions with those for the GOES I/M filter radiometer. The 15 mhi. 6.3 ym, and 4.3 designations correspond to the 
600-1 100 cm' 1 , ! 1 00-2000 cm “ 1 , and 2000-2300 cm"* regions, rcspectively. < H: o and 'dry refer to the atmospheric trans- 
mittance for water vapor and the uniformly mixed (e g., COj, N ; 0) optically active gases, respectively. 
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Fro. 3. (a) Spectra showing how emissions from different atmospheric levels arc smeared by low-resolution 
measurements, which cannot resolve the COj lines, (b) Improvement in vertical temperature profile resolution 
from modifying GOES I/M to use an interferometer. 


performed by Fleming and Barnes (1987). (Experi- 
mental confirmation of these theoretical expectations 
is shown in section 4.) As can be seen from Figs. 3 and 
4, dramatic improvements in vertical resolution and 
accuracy will result from the interferometer modifi- 
cation. 

Another very important application of the GOES 
interferometer will be the tracing of fine scale features 
of the atmosphere’s water substance. As will be shown 
in section 4, the interferometer sounder is capable of 
sensing water vapor emission with very high vertical 
resolution. Observations from a geostationary satellite 
of the dynamic displacement of fine vertical scale water 
vapor features might be used to construct vertical pro- 
files of horizontal wind velocity as well as to provide 
a direct measure of the moisture convergence respon- 
sible for weather development. This wind profiling ca- 
pability could be extremely useful over coastal oceanic 
regions, which are devoid of ground-based wind ob- 
servations. The moisture structure sensing capability 
of the interferometer, together with its associated im- 
proved temperature profiling performance, are consid- 
ered crucial for the use of geostationary satellite data 
in the mesoscale numerical weather prediction models 
planned for the next decade and beyond. 

In addition to greatly improving atmospheric 
sounding resolution and accuracy, the modified GOES 
Sounder will provide more accurate measurements of 
cloud altitudes ( Smith and Frey 1 990 ) and atmospheric 
soundings beneath thin and/or broken clouds. More- 
over, detailed observation of the infrared spectrum 
yield radiances highly sensitive to the concentrations 


of atmospheric gases such as ozone, nitrous oxide, car- 
bon monoxide, methane, and Freon 1 1 and 12 (Rev- 
ercomb et al. 1988c; Clough 1988), as well as the ther- 
mal emission of land and sea surfaces. The interfer- 
ometer sounder, once implemented aboard the 
geostationary satellite, will heip satisfy observational 
requirements for climatology, oceanography, and 
geology as well as operational meteorology. 

3. GHIS instrument characteristics 

One practical approach for implementing the GOES 
HIS is to modify the existing GOES I/M filter wheel 
sounder (Koenig 1989). This approach would permit 
its implementation as early as possible, in order to 
achieve the urgently needed improved sounding per- 
formance. Feasibility studies confirm that the replace- 
ment o t the GOES sounder filter wheel subassembly 
with an interferometer module is relatively straight- 
forward and that it is a very attractive way to obtain a 
major enhancement in vertical sounding resolution and 
accuracy. 

The primary operating characteristics of the GOES 
sounder as described by Koenig ( 1989) will be retained. 
The only functional change will be to expand the flex- 
ibility of the scan system to permit an optional longer 
dwell time at each sounding location, and another op- 
tional skip-step mode to permit faster soundings over 
large areas. As a consequence of the richer dataset col- 
lected by the interferometer, the output data rate will 
increase from 40 000 bits per second to as much as 
500 000 bits per second. 
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Fig. 4. Improvements in vertical temperature profile (a) and water 
vapor profile ( b ) rms retrieval errors from modifying GOES L/M 
to use an interferometer. 


Key features of the GHIS are given in Table 1 . Items 
modified or added for the new system are starred. It 
may be noted that these changes are a small part of 
the total system, indicating the compatible nature of 
the interferometer. The interferometer module will fit 
in place of the GOES I/M filter wheel assembly, leaving 
the telescope and detector-cooler assemblies essentially 
unchanged. 

a. GHIS operating modes 

The current filter wheel design for GOES I/M con- 
tains two scene scanning modes. Figures 5a, b show the 
scan patterns for the existing single step, normal line 
and single step, skip line modes. Two new scene scan- 
ning patterns have been introduced to fully utilize the 
capabilities of the interferometer while maintaining 
geographical coverage. The new scan patterns (skip 


step, normal line and skip step, skip line) are shown 
in Figs. 5c,d. 

In any one of the four scene scanning patterns, the 
sounder operator has the choice of three vertical (spec- 
tral) resolution modes: low, medium, or high. The three 
vertical (spectral) resolution modes are defined in Table 
2 in terms of the unapodized resolution attained in 
each band. The normal step, low-resolution mode is 
designed to provide the same area coverage currently 
achievable with the filter wheel radiometer, but with 
improved spectral resolution and dramatically lower 
noise performance in the long wavelength band (Table 
3b). The higher resolution modes require longer dwell 
times to allow the interferometer moving mirror to 
traverse the longer optical path range. 

The definition of these operating modes is based on 
a desire to maintain most of the operational charac- 
teristics of the GOES I/M sounder. 

b. Interface to existing system 

The modularity of the GOES sounder makes im- 
plementation of the interferometer modification quite 
feasible. The existing scene-mirror optics and its control 
circuitry will remain unchanged. Some reprogramming 


Table I. Key features of sounder. 


Simultaneous FOVs 

4 

FOV denning element 

Field stop, 242 ^rad diam (8.6 km, 
nadir) 

Telescope aperture 

31.1 cm (12.25 in) diameter 

Channel separation. LW- 

Dichrotc 

SW-MW 

IR spectral definition 

Interferometer* 

Visible sounding 

Fixed detectors 

Radiometric calibration 

Space and 290 K IR blackbody 

Frequency of space clamp 

2 min 

Frequency of IR calibration 

20 min 

IR detector operational 

92*, 97*, 102, 107 

temperature (K) 

Field sampling 

Four arras N-S on 10-km centers 

Scan step angle 

280 Mfad (10 km nadir); Skip-step* 
and skip-line optional 

Step and dwell time 

0.1 Seconds; 0.2 s, 0.4 s, 0.8* 
optional 

Scan capability 

Full earth and space 

Sounding arras 

10 km X 40 km to 60°N-S and 
60*E-W 

Optical location 

Star sensing 

Output data quantizing 

12 bit ail channels* 

Output data rate 

<500 kbits per second without 
data compression (lower by 2 to 
5 times with compression)* 

System power average 

112 W 

System weight: 

Sensor assembly 

202 lb (92.0 kg) 

Electronics module 

61 lb (27.9 kg) 

Power supply 

14 lb (6.2 kg) 

Interferometer processor* 

16 lb (7.3 kg) 

Total 

293 lb (133.4 kg) 


• Added for interferometer. 
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Figure 5a. Single step, normal line 


Figure 5b. Single siep, skip line 
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Figure 5c. Skip siep, normal line 
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Figure 5d. Skip siep, skip line 


FiG. 5. Sounding in different scan modes (dwell limes of 0.1. 0.4, 0.8 s are selectable at each 
location). The numbers indicate the order by which the four simultaneous fields of view are 
sampled. 


of the scene mirror will be required to accommodate 
the new scan modes; however, the GOES design is well 
suited to accepting such changes. The existing passive 
cooler, used to cool the aft optics, detectors, and focal 
plane will be maintained. Also, the physical dimensions 
of the sensor assembly will remain unchanged. The 
filter wheel passive cooler will be used to cool the in- 
terferometer module. The major change is to remove 
the niter wheel and relay optics from the input beam 
and to replace them with an interferometer and its relay 
optics. Another change is to reolace the filter wheel 
electronics circuitry with circuitry appropriate for con- 
trolling the dynamically aligned interferometer and for 
on-board processing of the data. The last change is to 
replace the current set of detectors with ones optimized 


for the interferometer’s long-wavelength (LW), mid- 
waveiength (MW), and short-wavelength (SW) band 
passes. 


The mass of the interferometer and the mass of the 
electronic control circuitry being added ( 10 lbs) nearly 
equals the mass of analogous material being removed 
(11 lbs). The main additional mass (16 lbs) being 
added in the interferometer design is for the on-board 
processing system. 

The total system mass is expected to be only about 
5% greater than the mass of the current filter wheel 
sounder. 


- — Cr-Mass^power. and volume 


Table 2. Vertical resolution modes. 


Spectral 

resolution 

Dwell ume 
(s) 

Time/step 

(s) 


Resolution (cm -1 ) 


Band I 

(620 to 1150 cm' f ) 

Band 2 

( 1 2 1 0 to 1740 cm’ 1 ) 

Band 3 

(2150 to 2721 cm' 1 ) 

Low 

0.075 

0.1 

3.33 

4.63 

4.63 

Medium 

0.375 

0.4 

0.67 

1.15 

1.15 

High 

0.775 

o.s 

0.32 

L15 

1.15 


ORJGjJV ;: % l p.4Q£ |5 
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Table 3a. Expected GHIS NEN performance. 


Spectral 

resolution 

Dwell time 
(s) 

Time/step 

(s) 


NEN (mW/m : srcnT 1 ) 


Band i 

(620 to 1150 cm _i ) 

Band 2 

(1210 to 1740 cm -1 ) 

Band 3 

(2150 to 2721 cm") 

Low 

0.075 

0.1 

0.14 

0.06 

0.007 

Medium 

0.375 

0.4 

0.32 

0.13 

0.014 

High 

0.775 

0.8 

0.46 

0.13 

0.014 


Table 3b. GHIS low resolution mode NEN performance reduced to GOES I/M resolution. 

Channel 

number 

Center 

wavenumber 

Resolution 

dv 

GOES I/M 

NEN 

specified 

NEN 

desired 

GHIS 

NEN 

expected 

1 

680 

13 

0.66 

0.25 

0.083 

2 

696 

13 

0.58 

0.25 

0.084 

3 

711 

13 

0.54 

0.25 

0.085 

4 

733 

16 

0.45 

0.25 

0.080 

5 

748 

30 

0.44 

0.2 

0.060 

6 

790 

50 

0,25 

0.15 

0.048 

7 

832 

50 

0.16 

0.06 

0.051 

8 

907 

25 

0.16 

0.06 

0.079 

9 

1035 

25 

0.35 

0.12 

0.090 

10 

1345 

55 

0.16 

0.05 

0.021 

II 

1425 

80 

0.12 

0.04 

0.019 

12 

1535 

60 

0.15 

0.04 

0.023 

13 

2188 

23 

0.013 

0.004 

0.0037 

14 

2210 

23 

0.013 

0.004 

0.0038 

15 

2245 

23 

0.013 

0.004 

0.0038 

16 

2420 

40 

0.008 

0.003 

0.0032 

17 

2513 

40 

0.0082 

0.003 

0.0032 

18 

2671 

100 

0.0036 

0.002 

0.0022 


The modification requires an additional 10.7 W of 
average power, including 5.2 W for the data processing 
on board. The interferometer adds no additional peak 
power requirements, and the additional average power 
is within the current margin for power available to the 
sounder. 

d. Detector performance 

Table 3 summarizes the expected noise performance 
for each spectral resolution mode specified in Table 2. 
This noise performance has been shown to meet the 
needs of high resolution sounding (Smith et al. 1983a), 
and in the low resolution mode will offer a substantial 
improvement over GOES I/M for the long wavelength 
band (Table 3b). 

Note that the NENs for soundings will be better than 
those shown in Table 3 by a factor of the square root 
of the number of field of views (FOVs) used per 
sounding retrieval. 

e. Data rate and processing 

The expected dataware is based on the use of Digital 
Signal Processor ( DSP) modules to do on-board Fou- 


rier transforms of the interferograrn data in real time. 
Table 4 gives the data rates for the full spectra (all 
bands included). For each item in the table, the higher 
number holds if complex spectra are needed for phase 
corrections of the data (Revercomb et al. 1988a), the 
lower number if only magnitude spectra are used. In 
all cases, the data rates are less than 500 kbps. Further 
reductions in data rate (by a factor of 2 to 5) could be 
realized with data compression. Another factor of 5 to 
10 could be achieved if the radiance spectra were con- 
voluted with a “standard atmosphere'’ weighting func- 
tion matrix, as required for sounding retrieval using 
the linear simultaneous approach (Smith et al. 1990) 
prior to telemetry to the ground (see section 4). In 
summary, on-board processing will enable the sound- 
ing information to be compressed to fit almost any 
prescribed data rate constraint. 

4. Aircraft prototype instrument 

In order to prove the theoretical expectation of im- 
proved sounding performance (Smith et al. 1979, 

1 983a^ Smith 1984), NOAA and NASA jointly funded 
the development of an aircraft prototype High-reso- 
lution Interferometer Sounder (HIS). The instrument 
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Table 4. Data rates for the interferometer sounder 


by operating mode. 


Spectral 

resolution 

Total data rate (kbps) 

Single step, normal line/ 
skip step, skip line 

Skip step, normal line/ 
skip step, skip line 

Low 

452/226 

226/113 

Medium 

‘452/226 

362/181 

High 

336/168 

299/149 


was completed and tested in 1985 and since then has 
flown on more than 40 airborne missions of the NASA 
U2/ER2 high altitude (65 000 feet) aircraft. All of its 
research flights to date have been flown on the NASA 
U2/ER2, but the instrument was recently checked out 
on the NOAA P-3 research aircraft, which has excellent 
in situ capabilities for temperature, trace gas and aero- 
sol measurements. Here we summarize the character- 
* istics and the radiometric performance of the aircraft 
HIS instrument, and present examples of the sounding 
results. 

a . Description 

The specific hardware configuration of the HIS air- 
craft instrument bears little resemblance to the design 
for GOES, but the principles are the same. The aircraft 
HIS views directly downward from inside a pod (about 
3 m long and 0.5 m in diameter) mounted either under 
the wing or under the center line of the fuselage. Many 
of the important design parameters of the HIS are 
summarizes in Table 5. 

Calibration is accomplished by viewing two high 
emissivity blackbodies. servo controlled at altitude to 
300 K and about 240 K After 12 scans of the earth 
over the full range of optical path, a 45 degree scene 
switching mirror rotates the fieid-of-view from the open 
earth viewing port to give four scans ol the hot and 
four scans of the cold blackbodies. The blackbodies, 
built and tested by Epplev Labs, are blackened cavities 
with thermoelectric cooler/ heaters for temperature 
control and platinum resistance thermometers for 
monitoring. The temperatures of the interferometer 
optics are nor actively controlled. 

The BOMEM Michelson interferometer (Buijs 
1979) as modified for this application provides double- 
sided interferograms from both scan directions. Its 
auto-alignment system makes it possible to operate in 
the ambient thermal environment of the pod and in 
very close proximity to the aircraft jet engine. The op- 
tical bench is shock-mounted to dampen high fre- 
quency vibration and the interferometer is evacuated 
to protect the beamsplitter during descent. 

The three spectral bands, covering most of the region 
from 3.8 to 16.6 microns (Table 5), are split inside a 
single liquid helium dewar, which contains three sets 
of bandpass cold filters, focusing optics, and arsenic- 


doped silicon detectors. The preamplifiers are external 
and operate near the ambient pod temperature of about 
260 K. The gain of each channel is fixed and the signals 
are digitized with a 16 bit A/D. On-board numerical 
filtering is used to reduce the sample rate from the 
HeNe laser rate by factors of 14, 8 and 8 in bands I, II 
and III. 

The data system is controlled with a 6809 micro- 
processor-based system built at the University of Den- 
ver. The three bands of interferometer data and house- 
keeping parameters are combined and recorded on 
formatted cassette tapes. Two drives with a capacity of 
67 megabytes each are used to provide 9 hours of con- 
tinuous recording time. 

Processing of selected data in the field is performed 
on IBM-compatible personal computers. Data is trans- 
feree! to hard disk and is processed with custom soft- 
ware, wftich displays the measured interferograms and 
corresponding spectra, and performs calibration to 
yield radiance or brightness temperature spectra. The 
calibration procedure uses full complex spectra to avoid 
errors that can arise from radiance emitted by the warm 
interferometer (Revercomb et al. 1988a,b). 


Table 5, Characteristics of the HIS aircraft instrument. 


Spectral range (cm* 1 ) 

Band I 

590-1080 

Band I! 

1080-1850 

Band ill 

2000-2700 

Field of view diameter (mr) 

Telescope 

100 

Interferometer 

30 

Blackbody reference sources 

Emtssmty 

>0.998 

Aperture diameter (cm) 

■ 1.5 

Temperature stability (K) 

— 0T"“ 

Auto-aligned interferometer 

modified Bo 


~ BBDA2.I 

Beamsplitter 


Substrate 

KG 

Coatings { 1/4 g at 3.3 mm) 

Ge + Sb:Sj 

Maximum delay (double-sided) current 
configuration (cm): 

Band i (hardware limit is — 2.0) 

-1.8 

Bands II and III (limited by data system) 

+ 1.2, -0.8 

Micheison mirror optical scan rate (cm $"*): 

0.6 

Aperaturc stop (at interferometer exit window): 

Diameter (cm) 

4.1 

Central obscuration area fraction 

0.17 

Area (cm*) 

10.8 

Area-solid angle product (enr sr): 

0.Q076 

Detectors: 

Type 

Ar doped Si 

Diameter (cm) 

0.16 

Temperature (K) 

6 
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b . Radiometric performance 

The rms detector noise for a single interferogram 
scan determined from in-flight calibration data is 
shown in Fig. 6. The noise is expressed in brightness 
temperature units for the calibration blackbody scene 
temperature of 260 K, For bands I and II the noise is 
background limited, with radiation from the instru- 
ment providing most of the photons. Band III noise is 
considerably higher, with an electrical interference 
contribution between 2150-2350 cm l , which is not 
detector noise. 

The other type of noise encountered in flight is sam- 
ple-position-error noise caused by the effect of aircralt 
vibrations on the velocity of the scanning Michelson 
mirror. Since it originates from very near zero path 
difference, this noise is highly correlated with wave- 
number, causing a_small rocking of the spectra. It is 
absent on the ground and would not be present in a 
spacecraft application. Even in the hostile aircralt en- 
vironment, this noise can be made small. Itsampiitude 
varies with instrument configuration, but is generally 
comparable to the detector noise level. 

Accurate radiometric calibration over the full spec- 
tral range was demonstrated on the ground using a 
third blackbody. A liquid nitrogen blackbody was used 
as the cold reference, one on-board blackbody served 
as the 300 K hot reference, and the temperature of the 
other on-board blackbody (set between 260 and 280 
K) was determined. The unknown temperature can be 
determined routinely to within about 0.2 to 0.4 K. This 
procedure relies on a careful optical design and align- 
ment to prevent optical-path-diffcrence dependence of 
the rcsponsiviiy, which would degrade the high reso- 
lution integrity of the spectra. Self-apodization is kept 
small by the relatively small ficid-of-view of the inter- 
ferometer. and can be accounted for accurately. Ad- 
ditional errors could be present in flight, but compar- 
sons with aircraft altitude temperatures and water sur- 
face temperatures are generally within i K. An accuracy 
of better than 0.5°C is possible, except for regions of 
low brightness temperature in the 4.3-^m CO: band. 

Figure 7 shows a portion of the HIS spectra recorded 
15 April 1986 over the desert Southwest compared to 


a spectrum calculated using radiosonde data with the 
AFGL FASCOD2 line-by-line program (Clough et al. 
1986). The finite FOV of the interferometer is ac- 
counted for in the calculated spectra. The agreement 
is remarkably good. The major differences are believed 
to be due to incomplete information on the tempera- 
ture, water vapor, and ozone (the large discrepancy 
between 1010-1070 cm' 1 ) of the observed atmosphere, 
deficiencies in the line-by-line calculations, and the 
presence of constituents not included in the calcula- 
tions ( e.g., CFC- 1 1 and CFC- 1 2 ) that contribute to the 
discrepancies seen between 840-860 and 910-925 
cm' 1 in Fig. 7. Noise and calibration errors are not 
believed to be significant. 

Another comparison of HIS observations with cal- 
culations is shown in Fig. 8. The observations were 
taken from the NOAA P-3 aircraft over the Gulf of 
Mexico on 18 November 1988. This flight provided 
the best supporting measurements of temperature and 
humidity of any of the HIS observations, via both air- 
craft sensors and dropsondes. The spectral agreement 
is generally within two degrees with the exception of 
the 600-680 cm' 1 region, which is influenced by 400 
Hz noise from the aircraft power supply, and the 2260- 
2380 cm' 1 region, which is due to the large influence 
of the electrical and detector noise noted earlier for this 
spectral region on brightness temperature when viewing 
low scene temperatures. Note especially the excellent 
agreement over the 6.3 micron water vapor band oc- 
cupying most of the HIS Band II. 

c. Example experimental sounding results 
1) Retrieval methodology 

For the retrieval of atmospheric profiles from the 
HIS, it is desirable to utilize all spectral observations 
in a simultaneous solution for temperature, water vapor 
and other desired absorbing constituents (e.g., ozone, 
methane, etc). Such a general solution has been for- . 
mulated for treating the HIS spectra and is described 
in detail elsewhere (Smith et al. 1988, 1990). - 

For the results shown here, the following perturba- 
tion form of the “physical simultaneous' 1 retrieval al- 
gorithm is used 



Fjg. 6. Detector noise observed with the HIS aircraft instrument. 
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Fig. 7. Blackbody equivalent temperature of upwelling radiance from the nadir viewing HIS instrument aboard the 
NASA U2 at an altitude of 20 km over the desert southwest on 1 5 April 1 986. Also shown is a line-by-line calculation 
with AFGL FASCOD2 using a nearby radiosonce. The spectra are apodized to yield an approximate resolution of 0,7 
cm - 


t = DAo\ = Dt* (1) 

where t is a vector of the deviations of the desired at- 
mospneric variables (e.g., temperature and dewpoint) 
from an ensemble mean condition, t* is a vector of 
observed brightness temperature deviations from the 
mean condition, Ao r is the transpose of the profile 
weighting function matrix, partitioned into "dry' 1 and 
“wet” components as required for the simultaneous 
solution, and D is a “deconvolution” matrix. In its 
application, we define 

t* 3 Ao T (t 4 + 0) (2) 

where d is the systematic error in calculated brightness 
temperature due to uncertainties in the spectroscopy 
and errors in the numerical methods used to define 
the atmospheric transmittance spectra. The elements 
of Aq are 

Af 

f- 1 

with 

w <i ~ f<) r ( ‘'i . Pj)d In r k {v h Pj'f/ff.i,,) 
and 

S, = dB(t, h Tj°)/dT \ 


D is Planck radiance, T is temperature, and M is the 
number of spectral radiance observations. The zero su- 
perscript refers to the mean temperature profile con- 
dition. The scripts ij, k denote spectral wavenumber, 
atmospheric pressure level, and atmospheric absorbing 
constituent, respectively. The transmittance of the at- 
mosphere between the instalment and pressure level, 
P)> is t( p } ) whereas that due to an individual ab- 
sorbing constituent (e.g., water vapor) is rd v n p ; ). The 
derivative is with respect to the vertical coordinate. 
Here we treat the spectrum as composed of uniformly 
mixed constituents (k = 1 ) and water vapor (k = 2). 
Strictly speaking, the first half of the vector t is the 
atmospheric temperature profile which satisfies that 
portion of the observed radiance spectrum influenced 
by the uniformly mixed gases whereas the second half 
is the temperature profile which satisfies that portion 
of the observed radiance spectrum influenced by the 
water vapor, assuming that the mean water vapor pro- 
file used to define W tJ is correct. Thus, if the actual 
water vapor profile is equal to the mean profile, then 
_the_twQ temperature profiles, uniformly mixed gas and 
water vapor, are identical. For all other conditions, the 
actual moisture profile is a function of the discrepancy 
between the two temperature profile components 
(Smith et al. 1990). The dewpoint temperature can be 
approximated from a linear combination of the two 
temperature profiles. The quantity is the standard 
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Fig 8 Spectra for ail three bands of the HIS observed from the NOAA P3 aircraft at 7.J 
km altitude ( 400 mb ) on 1 8 November 1 988. The difference spectra show the excellent agree- 
ment between HIS observations and spectra calculated using FASC0D2. 
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deviation of the radiance noise due to measurement 
noise and error in theoretical radiative transfer com- 
putations. It is included to weight those spectral regions 
that have the lowest noise most heavily in the profile 
solutions. 

The vector t*, as defined by (2), can be thought of 
as the vertically smoothed temperature profile as de- 
picted by the convolution of the spectral weighting 
functions with the radiant brightness temperatures. 
Numerically, 

M M 

t*(Pj)= 2 w n% ( 3 ) 

i- 1 /■ i 

that is, the convoluted temperature profile deviation 
from the mean is given by the spectrally weighted-av- 
erage brightness temperature spectrum deviation from 
the mean spectrum. This form for t* is identical to the 
generalized iterative retrieval solution (Smith 1970) 
used for processing grating spectrometer and filter 
wheel radiometer data. However, because of the large 
number of spectral radiance observations (~4000) 
with the HIS, iteration of (3) is prohibitive due to the 
very large amount of radiative transfer computation 
involved to update and the mean brightness tem- 
perature spectrum with the spectrum computed from 
a prior retrieval. Instead, t* serves here as the smooth 
first approximation of the true profile whose enhanced 
vertical structure is retrieved through the application 
of the deconvolution matrix, D. The D matrix is de- 
termined by matrix inversion of ( 1 ) 

D = (TT* r )(T*T* T )-' (4) 

where T is a matrix of dry-bulb temperature and dew- 
point temperature profile deviations from the mean 
for a statistical ensemble of soundings and T* is a ma- 
trix of the smoothed temperature (uniformly mixed 
gas and water vapor component) profiles from their 
ensemble mean values. T* is specified by radiative 


transfer computation of and /j ; . Random errors of 
measurement of the same magnitude observed from 
calibration spectra with the HIS aircraft instrument 
are simulated using a random number generator and 
added to Iti in order to stabilize the solution against 
actual instrument noise. For the generation of T* used 
in (4), £ is assumed to be equal to zero since here T* 
is based solely on radiative transfer theory. However, 
when D is applied to real spectral brightness temper- 
ature observations, 0 must be applied to t* [Eq. (3)J 
in order to account for systematic errors in the radiative 
transfer model. 

During the COmbined Huntsville Meteorological 
Experiment (COHMEX), the HIS flew aboard the 
NASA U2/ER2 aircraft. A few sample results achieved 
during COHMEX are presented beiow to illustrate the 
sounding capabilities of the HIS. 

2) Simulation results 

Soundings have been derived using HIS radiance 
spectra simulated from the special radiosondes 
launched during the COHMEX. The purpose of first 
investigating retrieval accuracy using synthesized data, 
rather than actual spectral observations, is that the 
“truth” is known since the radiosonde profile was used 
to compute the “observed” spectrum. ( Random errors 
of the s: me magnitude associated with actual HIS ob- 
servations are added to the computed radiances to 
simulate “observed" radiances for the retrieval process.) 
As shown later, when comparing profile retrievals from 
actual spectral radiance observations with radiosondes, 
it is uncertain as to whether the difference is due to 
errors in the retrieval, the radiosonde, or space and 
time discrepancies between the two independent ob- 
servations. Thus, the use of synthetic spectra serves to 
illustrate the errors expected to result from the sounding 
retrieval process. 




Fig. 9. Rms retrieval errors for (a) temperature and (b) dewpoint compared to standard 
derivation of the 193 COHMEX radiosondes used in the retrieval error analysis. 
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Fig. 10. Comparison of rms errors from simulated temperature 
retrievals using all GHIS spectral bands (a), omitting the 4.3-um 
COj band (b), and omitting the 15-um CO; band (c). The long 
wavelength band is important at ail altitudes, while the short wave- 
length band only contributes significantly near the surface. 


Figure 9 shows rms errors of the sounding retrievals 
from HIS radiance spectra synthesized from a set of 
193 COHMEX radiosonde profiles. The temperature 
and dewpoint soundings were retrieved using decon- 
volution matrix. D, based on the same COHMEX sta- 
tistical set of 193 radiosonde soundings, and as a result. 
Fig. 9 represents the lower limit of retrieval error using 
the linear simultaneous retrieval algorithm. As can be 
seen for the COHMEX observation region, HIS re- 


trieval accuracies close to 1 °C could be achieved for 
temperature, and 2°-3°C for dewpoint temperature. 
Figure 10 shows rms error curves derived from a global 
set of soundings (Smith et al. 1990) for three different 
measurement conditions: l ) all spectral bands ( 600- 
2700 cm -1 ), 2) bands 1 and 2 (600-2000 cm' 1 ), and 
3) bands 2 and 3 ( 1 100-2700 era -1 ). Best results are 
achieved when all spectral measurements are used in 
the retrieval. 

Figure 1 1 shows two example comparisons of syn- 
thetic radiance retrievals with individual COHMEX 
radiosondes. As can be seen, much of the fine scale 
detail is retrieved from the full HIS radiance spectrum. 
The vertical structure of water vapor dewpoint is par- 
ticularly noteworthy. 

3) Experimental results 

Figure 12 shows two examples of comparisons of 
radiosondes with retrievals of temperature and dew- 
point profiles from actual HIS spectral radiances ob- 
served on 19 June 1986. It is noted that the D matrix 
and initial (mean) profile used for these retrievals are 
the s am e as those used for the simulated radiance re- 
trievals shown earlier. We can see from these compar- 
isons that the agreement of the radiance derived profiles 
with the radiosonde is close to the theoretical expec- 
tation (Fig. 9). Once again, the fine-scale water vapor 
structure retrieved is particularly noteworthy. 

Figure 13 is presented to show the improvement in 
vertical resolution and accuracy of HIS profile retrievals 
over those achieved with the current VAS geostationary 
satellite sounder. The GOES I/M filter wheel sounder 
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FiC. 11. Examples of two COHMEX radiosondes, (a) and (b) above, compared to the retrievals of temperature and dewpoint 
obtained from synthetic radiances calculated from them. Notice that most of the fine scale structure is recovered, particularly for 
water vapor. 
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Fig. 12. Examples of two COHMEX radiosondes, (a) and (b) above, compared to retrievals 
of temperature and dewpoint from actual HIS spectral radiances observed on 19 June 1986. 


will produce soundings of vertical resolution similar to 
VAS. As can be seen, whereas VAS depicts the broad 
scale profiles of temperature and moisture, HIS is ca- 
pable of resolving finer scale features of vertical struc- 
ture. Once again the most noticeable improvement ap- 
pears in the moisture profile retrieval. The moisture 
discrepancy between the HIS and radiosonde moisture 
above the 600-mb level is believed to be due to space 
and time discrepancies of the measurements rather 
than due to errors in either sounding. Evidence for this 
arising from strong space and time gradients of at- 


mospheric moisture in this layer is given in Fig. 14 
described below. 

In order to illustrate the water vapor sensing capa- 
bility of the HIS, Figs. 14 and 15 show vertical cross 
sections along the flight tracks of the NASA U2 and 
ER2 on 15 and 19 June 1989, respectively. Both flight 
tracks are over northern Alabama and Tennessee, but 
as can be seen from the HIS data, the meteorological 
characteristics are quite different for the two days. On 
15 June the HIS cross section reveals an intense 
boundary layer whose top was near the 850-mb level. 
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Fig. 13. The improvement in vertical resolution of the HIS (a) is apparent when compared 
to a coincident satellite retrieval using the VAS filter wheel sounder (b). 
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Indeed, cloud observations revealed that the top of a 
low-level scattered cumulus deck was near the 850-mb 
level. The rapid drying and stabilizing temperature 
structure just above the boundary layer ( 850-700 mb ) 
is clearly evident in the HIS cross-section analysis. Also, 
the largest horizontal variability in moisture occurs 
within the 200-400 mb layer where the aircraft tra- 
verses a moisture front located across central Tennes- 
see. On 19 June the meteorological condition is quite 
different. The HIS cross section does not reveal the 
existence of an intense boundary layer and this was 
supported by the fact that no boundary layer cumulus 
clouds or haze were observed. Contrary to the 15 June 
observation, the largest horizontal variability of at- 
mospheric moisture occurs lower in the atmosphere 
within the 400-700 mb layer. In both cases, the vertical 
resolution of the retrieved water vapor structure is 
striking. In this regard, the HIS sounding capabilities 
are unique; that is, no other sensor is capable of pro- 



ms TEMPERATURE ANO DEWPOINT DEPRESSION 05 JUNE 1986) 


Ftc. 14. (a) Flight track of the NASA/U2 over northern Alabama 
and Tennessee on 15 June 1986. The urnes of various aircraft are 
denoted in hours, minutes, and seconds of central daylight time, (b) 
Vertical/Ume cross section of atmospheric temperature (dashed) and 
dewpoint depression (solid) beneath the U2 flight track as retrieved 
from HIS spectra. 




HIS TEMPERATURE ANO DEWPOINT DEPRESSION (19 JUNE 1986) 
Fig. 15. As in Fig. 14. but for 19 June 1986. 

viding the vertical and horizontal sounding detail as 
achieved with the airborne HIS. From the geostauonarv 
satellite, the HIS will provide an unprecedented four- 
dimensional sounding capability. 

5. Conclusions 

A design exists to enable implementation of a greatly 
improved sounding capability from geostationary sat- 
ellite orbit. Experimental results achieved with an air- 
craft prototype instrument indicate that the theoretical 
expectations of 1 °C temperature and 2°-3°C dewpoint 
accuracy can be achieved. Most exciting is the great 
improvement in sensing the water vapor profile struc- 
ture. It appears that the vertical resolution is good 
enough to enable water vapor tracking in numerous 
vertical layers thereby providing water vapor tracer 
wind profile data as well as thermodynamic profiles. 
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1. Introduction 

The purpose of this paper is to demonstrate the trace 
gas profiling capabilities of future passive high spectral 
resolution (1 cm 1 or better) infrared (600 to 2700 cm -1 ) 
satellite tropospheric sounders. These sounders, such as the 
grating spectrometer, Atmospheric InfRared Sounders (AIRS) 
(Chahine et al., 1990) and the interferometer, GOES High 
resolution Interferometer Sounder (GHIS) , (Smith et al., 

1991) can provide these unique infrared spectra which enable 
us to conduct this analysis. 

In this calculation only the total random retrieval 
error component is presented. The systematic error 
components contributed by the forward and inverse model 
error are not considered (subject of further studies) . The 
total random errors, which are composed of null space error 
(vertical resolution component error) and measurement error 
(instrument noise component error), are computed by assuming 
one wavenumber spectral resolution with wavenumber span from 
1100 cm -1 to 2300 cm -1 (the band 600 cm -1 to 1100 cm' 1 is 
not used since there is no major absorption of our three 
gases here) and measurement noise of 0.25 degree at 
reference temperature of 260 degree K. Temperature, water 
vapor ozone and mixing ratio profiles of nitrous oxide, 
carbon monoxide and methane are taken from 1976 US Standard 
Atmosphere conditions (a FASCCDE model) . Covariance matrices 
of the gases are "subjectively" generated by assuming 50% 
standard deviation of gaussian perturbation with respect to 
their US Standard model profiles. 

Minimum information and maximum likelihood retrieval 
solution s a re u sed. All terms and definitions except 
weighting functions outlined in Sec. 2, are identical to 
those defined in Rodgers (1989), and are not repeated here. 
Retrieval random error is presented in terms of "level 
explained variance" for each individual trace gas. However, 
these preliminary results are by no means complete and are 
only to be considered as "illustrative". 


2. Outline of the Computations 


The transmittance is computed from a fast regression 
model (developed at CIMSS) which is based on line-by-line 
FASCODE 2. Component transmittance profile spectra are 
created for each radiatively absorbing gas such as H 2 O, O 3 , 
N 2 O, CO, CH 4 , and CO 2 plus all other absorbing gases. The 
effective temperature weighting function is then defined as 
described in Smith et al. (1991), which has the form of 

Ki(p) = S(p)r(p)dlnrj.(p) (1) 

where i represents ith absorbing gas, p denotes pressure, r 
is the total transmittance, rj_ is the component 
transmittance of ith gas and 3 denotes scaling factor among 
radiance, brightness temperature and air temperature. The 
wavenumber dependence of ( 1 ) is to be understood. 

Following the development of Rodgers (1989) , total 
random error is the sum of null space error and measurement 
error. The null space error covariance matrix is 

S N = (A-I)S X (A-I) T (2) 


and the measurement error covariance matrix is 

S M = DS s D T ( 3 ) 

where the averaging kernel A=DK and D is the inverse model 
dependent contribution function that will be defined later. 
The total random error covariance matrix is 


S T = S N “ S M (4) 

The "level explained variance" can be expressed as one minus 
the diagonal element of Srp divided by the diagonal element 
of S x for each gas x at each level. Note the neglect of the 
off-diagonal elements of Sp. 

For the minimum information solution, the contribution 
function is 

D = (K T K - 7 I) “ 1 X T ( 5 ) 

where 7 is arbitrarily set to io~ 3 . For the maximum 
likelihood solution, 

D = (K T S e _1 K + S x _ 1 ) - 1 K T S e -1 ( 6 ) 

3. Results and Discussions 

Figure 1 represents the retrieval explained variance 
profile for nitrous oxide, carbon monoxide and methane. Two 
curves in each panel represents the minimum information and 
maximum likelihood solutions, respectively. In general, the 
mid— atmosphere (400 to 850 hpa) has high explained variance 



for all three gases. The subtle differences in details 
between two solutions require further studies. 

Random error contributions from the null space and 
measurement components will be presented in the conference. 
The averaging kernels will also be discussed. The 
uncertainties of the a priori statistics and the effects of 
the smoothing parameter 7 on the component errors will be 
explored. The results presented here are preliminary. 
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Figure 1, Retrieval explained variance profile of nitrous 
:xide, carbon monoxide and methane for minimum 
information and maximum likelihood solutions. 
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ABSTRACT 

High resolution multispec tral aircraft data (MAMS) and 
geostationary sounding data (VAS) have been combined in a physical 
retrieval method to produce 100 meter resolution derived images of 
precipitable water, lifted index, and skin temperature for the COHMEX 
day of June 19, 1986. The method benefits by combining the 
multispectral sounding capability of VAS with the 100 meter spatial 
resolution of MAMS. The derived images - show significant small scale 
features including a moist tongue intruding into the COHMEX region and a 
heat island at Nashville, Tn. Structure function calculations indicate 
that significant mesoscale variability is not being captured by current 
geostationary sounding instruments. Geostationary sounding data with at 
least 1 km resolution or better is needed to monitor this additional 


mesoscale variation with time. 


1. 


INTRODUCTION 


Monitoring atmospheric moisture and stability at small 
scales is an important step to understanding and forecasting 
mesoscale processes. Few ground-based observing systems can 
provide the high temporal and spatial coverage necessary to study 
mesoscale variation; satellite observing systems have shown 
promise in this endeavor. Moisture and stability products from 
the geostationary platform VISSR Atmospheric Sounder (VAS) exhibit 
Interesting detail at the meso (20-200 km) scale (Chesters et 
al. 1986; Mostek et al . 1986; Smith et al . 1985; Hayden 1988). 
However, observations of variation at the meso 7 (2-20 km) scale 
remain elusive. Atmospheric processes involved in meso 7 scale 
weather events (such as airmass thunderstorms, cloud development, 
inversions and land/water interface circulations) are not well 
observed by VAS data because of its 8 km resolution and 
radiometric noise. Data averaging reduces the noise but also 
removes mesoscale variation. Therefore, information leading to a 
better understanding of such mesoscale phenomena must come from a 
higher resolution observing system. 

The Multispectral Atmospheric Mapping Sensor (MAMS) is well 
suited to provide mesoscale information about the atmosphere and 
surface of the earth. This instrument provides high resolution 
(-100 meter) data over a significant swath (-37 km) and has been 
used to depict mid- tropospheric water vapor (Menzel et al. 1986); 
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Jedlovec et al . 1986), total precipitable water (Jedlovec 1990) 
and sea surface temperature (Moeller et al. 1989) through the use 
of its 6pm, 11pm, and 12pm channels. Similar VAS channels have 
contributed significantly to precipitable water and lifted index 
products (Chesters et al. 1986; Robinson et al . 1986). However, 
Chesters et al., and Robinson e L al . also show that other VAS 
channels are important contributors to these products. This 
suggests that high resolution derived products be produced by 
combining multispectral VAS data and high resolution MAMS data. A 
physical retrieval algorithm developed by Smith et al. (1985) 
provides the framework to combine VAS and MAMS data by replacing 
VAS 6pm, 11pm, and 12pm spectral data with similar MAMS spectral 
data to generate the derived product images. 

This paper presents the procedure by which MAMS and VAS data 
are used to produce high resolution derived products of 
precipitable water (PW) , lifted index (LI), and skin temperature 
(TS) . Section 2 addresses the characteristics of the VAS and MAMS 
instruments. Section 3 describes the procedure used to make the 
derived products, while section 4 presents the application of the 
procedure to real data, with comparisons to radiosondes collected 
during the Cooperative Huntsville Meteorological Experiment 
(COHMEX) (Williams et al. 1987). The spatial and temporal 
variability of the derived products is evaluated in a structure 
function analysis (Hillger and Vonder Haar 1979). Data of 


June 19, 1986, a day of low level cumulus development over the 
COHMEX field experiment, is presented. 

2. INSTRUMENT DESCRIPTIONS AND COMPARISONS 

2.1. Both the MAMS and VAS instruments have been described 
extensively in the open literature so that only a brief discussion 
of each instrument is presented here. 

2.1.1. MAMS 

The MAMS instrument is a 12 channel scanning radiometer that 
flies onboard NASA's high altitude ER2 and U2 aircraft at an 
altitude of about 20 km. Figure 1 shows a depiction of the MAMS 
scanning geometry. The instrument scans across an 86 degree field 
of view for a swath width of about 37 km. At 20 km altitude, the 
data has 100 meter resolution when the 5.0 mrad aperture is used. 
The MAMS spectral bands are summarized in Table 1. Channels 1-8 
are visible/near infrared channels. Channels 9-12 are infrared 
channels. The infrared channels include surface viewing channels 
at 3.7pm and 11.2pm as well as water vapor channels for the mid 
(6.6pm) and lower troposphere (12.5pm). As noted in the footnotes 
of Table 1, only three infrared bands can be used on a given 
flight (the 3.7pm and the 6.6pm channels are readily exchanged 
between flights). Generally, the 3.7pm channel is used for 
surface applications (e.g., sea surface temperature); the 6.6pm 
channel is employed for atmospheric studies (e.g., atmospheric 
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water vapor). For this study, the 6,6^m detector was used in 
channels 9 and 10. Channels 11 and 12 form the well-known split 
window combination. 

MAMS visible/near IR data is calibrated by laboratory 
measurement, usually right before and after aircraft deployments. 
Infrared data is calibrated on each scan line by viewing two 
onboard blackbodies of known temperature, A linear detector 
response is assumed to define the calibration slope and offset. 
Blackbody counts may be averaged over several lines to reduce 
noise in the line to line calibration. 

Additional discussion of the MAMS instrument and its 
performance is found in Jedlovec et al. (1989). 

2.1.2. VAS 

VAS is a 12 channel spin scan radiometer that was deployed 
on the GOES series of satellites for the purpose of atmospheric 
sounding and multispectral imaging (Menzel 1980; Smith et al. 
1981). The VAS channels are shown in Table 2. VAS consists of 
temperature sounding channels in the 15^m CO 2 band (channels 1-5) , 
a split window combination (7 and 8) , mid- tropospheric water vapor 
channels (9 and 10), and shortwave channels (6, 11, and 12). In 
general, VAS channels 1-10 are used to make retrievals with 
channel 12 used for a skin temperature guess and cloud checks 
(Hayden 1988). VAS channels 4, 5, 7, 8, 9, and 10 are available 
at 8 km resolution (others at 16 km resolution) . Noise is 


4 



generally reduced to less than 1,5 K for all channels by multiple 
sampling (the dwell sounding mode of operation). Further 
information on the VAS instrument may be found in Montgomery and 
Uccellini (1985) . 

2.2, MAMS/VAS Instrument Comparisons 

The VAS instrument has several atmospheric sounding channels 
that the MAMS instrument does not. This additional spectral 
coverage provides VAS with a better vertical resolution sounding 
capability than MAMS has. The MAMS data has high horizontal 
resolution. This resolution advantage over VAS, while perhaps not 
so important for mid- tropospheric data, is well suited for 
evaluating the variability of surface and lower tropospheric 
fields. The spatial variability of MAMS ll^um and 12/irn data has 
been demonstrated by Jedlovec (1990) for June 19, 1986 COHMEX 
data. In that work thermal variance was shown to level off 
significantly beyond a field of view size of about 3 km. Thus, 
the MAMS shows a capability to identify small scale thermal 
features that are below the resolution of VAS data, which may be 
important in identifying preferred areas of local convection or 
low level cloud development. 

In replacing VAS spectral data with MAMS spectral data, it 
is important, because of absorption constituent considerations , 
that the MAMS and VAS channels have close spectral agreement. 
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Figure 2 illustrates the spectral characteristics of the similar 
channels. Both the 6pm and 11pm channels of VAS and MAMS show 
excellent agreement while the 12pm channels show some deviation 
which is more noticeable because of their narrow bandwidths , 

Water vapor absorption in the MAMS 12pm channel will likely be 
greater than in the VAS 12pm channel. However, this can be 
accounted for in a MAMS/VAS instrument bias adjustment. 

An instrument bias was calculated by comparing collocated 
VAS and MAMS observed data of June 19, 1986 over the entire MAMS 
data domain. This bias remov es variance between the instruments 
arising from spectral response, absolute calibration, viewing 
geometry, and time of observation differences, A set of 69 MAMS 
and VAS collocated data were considered for the bias calculation. 
VAS data observed at nominal times of 1600 and 1730 UTC were used 
to compare with the MAMS data collected between 1530 and 1800 UTC. 
In order to obtain spatial agreement in footprint size, the VAS 
spatial weighting function (which includes the effects of impulse 
response and diffraction) was applied to the MAMS data to simulate 
the spatial response of a VAS field of view. The results of this 
comparison showed biases of 2.21, 5.69, and 3 . 24 K for the 6pm, 
11pm, and 12pm comparisons, respectively. The positive values 
indicate that the MAMS data is warmer than the VAS data. These 
bias values are applied to the MAMS data when it is substituted 
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for the VAS data in the derived imagery production process which 
is described in the next section. 

3. ALGORITHM 

Previous efforts to make derived imagery products from 8 km 
VAS data have used retrieval statistics generated from a 
simultaneous retrieval technique (Smith and Woolf 1984) for an 
area average (e.g. , 100 km by 100 km) of VAS radiances. The 
retrieval statistics are then applied to individual VAS fields of 
view (fov) to produce 8 km resolution derived product Imagery 
(Smith et al . 1985). In the formulation, a set of coefficients, 
ct ^ , are produced and are used with pressure dependent basis 
functions, 0(p) , to produce departures of skin temperature T s , 
precipitable water vapor U, and atmospheric temperature T, from 
some first estimate: 

<5T S — Oq<£q 

N p 

5U(p) = 2 a i f q(p) 4>i (p) dp 

i-1 o 

M 

ST(p) - 2 aj: (p) 

i=N+l 

where precipitable water vapor U and mixing ratio q are related by 

N 

S q(p) - g S q (p) 6i (p) 
i-1 
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1 p 

6 (P) f 6 q(p) dp 

g o 

and g is gravity. See Smith and Woolf for the formulation of the 
at£ coefficients. Since is dependent on the channel brightness 
temperatures, a ^ can be recalculated at each fov to capture pixel 
to pixel variation in the skin temperature, precipitable water 
vapor, and atmospheric temperature fields. ct£ dependence on the 
atmospheric transmittance (including the transmittance dependence 
on water vapor), the slope of the Planck function (for converting 
radiance to temperature) , and the expected errors of the radiance 
observations are not recalculated at each fov (Smith and Woolf 
1984). Channel weighting functions are used as basis functions. 

In the MAMS/VAS effort, a three channel (6//m, ll^m, 12/im) 
retrieval is produced for a 5 x 5 VAS fov area, with 
coefficients then computed at each HAMS fov, using the bias 
corrected HAMS 6.6/xm, 11.2^um, and I2.5/im channel radiances in 
place of the similar VAS channels. A 5 x 5 retrieval size was 
chosen because it roughly matches the HAMS 37 km swath width. 

Those VAS channels not replaced by MAMS data influence the 
MAMS/VAS retrieval through its first guess, a VAS -only 5x5 fov 
10 channel retrieval (which uses the NMC Aviation model forecast 
for its first guess) . As with previous retrieval efforts (Hayden 
1988) , surface data from National Weather Service SVCA 
observations are used in the retrieval. This data is updated to 
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the resolution of the MAMS data by adjusting a grid of the SVCA 
data for MAMS 11/im temperature variation. Cloud checks are 
performed on the VAS and MAMS data. Cloudy VAS fovs are removed 
from the 5x5 retrieval. No soundings are made at cloudy MAMS 
pixels . 

4. EVALUATION OF TECHNIQUE 
4.1. J une 19, 1986 

The surface weather of June 19 over the COHMEX region is 
shown in Fig. 3, Dry Canadian high pressure which had moved over 
the COHMEX region on the 18th began to slide off to the east, 
reducing the subsidence over the region as well as increasing the 
low level moisture, especially in the southern part of the COHMEX 
region. Surface temperatures reached into the upper 80s while 
surface dew points ranged from the upper 50s to the upper 60s from 
north to south over the COHMEX region. Synoptic scale atmospheric 
forcings were weak as surface and upper level winds were light. 
Although no significant convection occurred, sufficient moisture 
and surface heating were present to produce low level cumulus over 
eastern Tennessee beginning about 1700 UTC. During the next 
several hours this cumulus field expanded to cover much of east 
central Tennessee and northern Alabama. During the afternoon a 
few of these cumulus did grow sufficiently to produce brief light 
showers over the region. 
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On June 19, the MAMS was flown over the COHMEX region from 
about 1530 UTC to 1800 UTC. Figure 4 shows the COHMEX region and 
the MAMS straight line flight tracks that were flown on that day. 
Table 3 lists information about each flight track. Note that 
flight tracks F-G and H-I were repeated. During the flight, 100 
meter resolution radiometric data with 8-bit precision was 
collected for the MAMS 6.6/im, 11.2/im, and 12 . 5/im channels. On 
this same day, VAS dwell sounding pairs were collected at 1600 UTC 
and 1730 UTC, These dwell sounding data were combined with the 
MAMS data as indicated in Table 3 to produce the derived imagery 
products . 

Both the MAMS and VAS navigation were inspected to assure 
good collocation between the two instruments. MAMS 6 . 6/xm data was 
corrected for zenith angle effects under an assumption of 
homogeneity across the MAMS swath width; no correction was applied 
to the MAMS 11.2/im or 12.5/im data because of a relatively weaker 
effect in these channels and uncertainty in the correction. 

MAMS/VAS retrievals were made using the procedure described 
in section III for five of the six straight-line flight tracks 
flown on June 19 by MAMS (the B-C Flight track was not processed 
because it fell in a VAS data void). Figure 5 shows a 
representative profile as it evolves during the retrieval process. 
The top panel of Fig. 5 shows the impact of the VAS-only retrieval 
(which is used as the first guess for the MAMS/VAS retrieval). A 
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slight modification of the NMC guess moisture profile has occurred 
while the temperature profile is affected at mid and upper 
tropospheric levels only. This temperature profile impact is at 
least partly a result of the VAS CO 2 sounding channels which are 
not common to the MAMS instrument. Without including these VAS 
channels, this modification to the guess profile probably would 
not occur. The bottom of Fig. 5 demonstrates the resulting 
profile after the MAMS/VAS retrieval. The guess profile shown in 
the bottom of Fig. 5 is the profile that was produced in the VAS- 
only retrieval. Changes to the guess profile are confined to 
lower levels of the troposphere. This reflects the impact of the 
MAMS llfim and 12/im data on the MAMS/VAS retrieval. Note that the 
mid and upper tropospheric profile has not been modified in the 
bottom panel . 

4.2. Analysis of Derived Products 

Fields of precipitable water (PW) , lifted index (LI) and 
skin temperature (TS) from the MAMS/VAS retrieval process are 
shown for June 19 in Figs. 6-8. The data has been remapped into a 
Mercator projection with a degradation of resolution down to . 6 km 
so that adjacent flight tracks can be shown in a single image 
format. Each figure consists of a composite of three initial 
flight tracks (D-E, F-G, H-I) and a composite of two repeated 
flight tracks (F-G, H-I). The northernmost two flight tracks of 


the three flight track image are collocated with the image of the 
repeated flight tracks . Note that there is about 30 to 40 minutes 
time difference between adjacent flight tracks which reduces 
continuity of clouds and features between flight tracks. There is 
also an indication of a limb effect across each flight track in 
the figures. 

The precipitable water ( PW) analysis shows generally higher 
moisture to the east and south in the imagery. Undetected thin 
cirrus exists in the NW corner of both the three track and two 
track composites and is likely responsible for the moisture 
maximum in that portion of the PW imagery. A moisture maximum of 
about 34 mm can be seen in east central Tennessee southeast of 
Nashville, TN (BNA) . The gradient to the northwest shows a 
decrease of about 8-10 mm over a distance of about 100 km with a 
more concentrated gradient of 6-8 mm over about 20 km to the 
northwest of the moisture maximum. This moisture pattern was also 
found by Jedlovec (1990) using a split window variance ratio 
technique on MAMS June 19, 1986 data. This general moisture 
pattern shows up in the COHMEX Beta Network radiosondes for 1800 
UTC in Fig. 9. It is evident however that the maximum and minimum 
locations in Fig. 9 are strongly influenced by the location of the 
COHMEX radiosondes, as would be expected; the high resolution 
MAMS/VAS product is much more caoable of accurately identifying 
the location of the moisture maximum than the -100 km spaced 
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COHMEX Beta Network radiosondes. Radiosonde PW values show fairly 
good relative agreement with the MAMS/VAS PW field, while absolute 
comparisons show a small (-2 mm) bias in the MAMS/VAS PW product 
(for this comparison, the repeated flight track data is used since 
it is the closest collocation in time to the launch of the 
radiosondes) . 

Higher moisture to the southeast is also identifiable in the 
repeated flight tracks which provide some insight to the time 
evolution of the field. The moisture has extended somewhat 
northward and expanded in the east-west direction in the repeat 
overpass composite. The analysis shows PW values have increased 
by 1-2 mm in east central Tennessee and that the gradient to the 
northwest has smoothed somewhat. This time continuity lends 
credence to the MAMS/VAS product. Low level cumulus began to form 
at the time of the repeat overpass (-1700 UTC) in east central 
Tennessee . In the next couple hours this field of cumulus 
expanded westward to about the center of flight track F-G, in good 
correlation with the increasing moisture. Figure 10 shows the 
cumulus field as it looked at 1900 UTC, about one hour after the 
MAMS data collection ended. 2100 UTC COHMEX radiosondes (not 
shown) indicated a continuing increase of moisture over east 
central Tennessee, in agreement with the trend depicted in the 
MAMS/VAS PW imagery. While the influence of other factors such as 
topography, surface conditions, etc., should not be overlooked, 
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atmospheric moisture seems to have played an important role in the 
development of the cumulus . 

The lifted index (LI) analysis (Fig. 7) of the initial 
flight tracks shows a broad unstable maximum (-2) in east central 
Tennessee near the location of the moisture maximum. This 
generally agrees with the relative variation of the COHMEX 1800 
UTC radiosonde LI, although a 3-4 K bias towards unstable values 
is apparent. The general location of the low level cumulus field 
that developed agrees with the relatively unstable location in 
east central Tennessee. The repeated flight tracks show that much 
of the area is destabilizing over the one hour interval. This is 
likely the result of surface heating. Interestingly, the most 
rapid destabilization is taking place in the western portion of 
the flight tracks. This might lead one to expect cloud 
development further west than actually occurred. However, the 
cumulus development remained in and around the moisture (Fig. 6) 
concentrated in east central Tennessee. 

Around the cumulus cloud development in east central 
Tennessee (repeat overpass of the F-G track) the LI values appear 
to remain constant or actually increase (become more stable) with 
time. This seems contrary to cloud development and contrasts to a 
general destabilization over the rest of the COHMEX region. 

Evidence is insufficient to say whether this is a true atmospheric 
trend or a result of cloud contamination in the VAS radiances, 
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which have a larger impact on the LI product than on the PW or TS 
products . However, this does point out the difficulty of using 
relatively low resolution data in a highly variant field (in this 
case impacted by clouds) . The need for high resolution data in 
atmospheric sounding for cloud clearing purposes has been pointed 
out by Menzel et al. (1984) in a study of VAS retrievals using 
large detector versus small detector data. 

The skin temperature (TS) product of Fig. 8 shows 
considerable high resolution variability. Land slope and the 
various land uses in the COHMEX region (e.g., forest, agriculture, 
urban) contribute strongly towards this variability. In 
particular, a definite heat island is observed at the location of 
Nashville, TN (BNA) with skin temperatures 5 to 10 K warmer than 
the surrounding region. Surface emissivity variation may also 
contribute to the variability in the skin temperature product 
since emissivity is assumed uniform at .96 for all the retrievals 
(Hayden 1988) . Note that individual pixel values of TS range by 
more than 20 °C. 

Time change of TS shows good continuity of thermal features 
in Fig. 8. Samples show that individual pixel skin temperatures 
increase by as little as < 1 K to as much as 5 K or more over the 
one hour time interval. Over a selected larger clear region 
(-50000 pixels; -20 X 25 km) the average skin temperature increase 
was 2.75 K. Many land types are included in this average, which 
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is thus not representative of the temperature changes that might 
occur at individual 100 meter resolution pixels. For comparison, 
Tarpley (1988) found 30-day average three hour (0900-1200 LST) 
heating rates of 4-8 K across Kansas in June using 40 X 40 km GOES 
target areas. Such skin temperature monitoring is useful in 
radiation budget studies of the surface of the earth. Significant 
small scale variability in the skin temperature field makes it 
difficult to identify relationships between cloud development and 
skin temperature at the meso 7 scale. In a study of Oklahoma land 
types and their relationship to convective cloud formation, Rabin 
et al . (1990) found convective clouds tended to form earliest over 
warm, dry land surfaces. In a broad sense, the cumulus on June 19 
formed over a region of relatively warm skin temperatures, 
although there was no cloud development over many other regions of 
comparably warm skin temperatures in the TS imagery. Given a lack 
of rainfall in Tennessee over the previous eight davs (only 
isolated convection along the Tennessee/Alabama and 
Tennessee/Georgia borders on June 17), it is unlikely that surface 
residual moisture was a significant factor in the location of the 
cumulus development. This indicates that other factors such as 
the atmospheric moisture pattern and topography played an 
important role in the cloud development on June 19. Topography 
data show that significant relief (500-1000 feet) exists in SE 
Tennessee and gradually diminishes towards the MAMS flight area to 
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the west (<500 feet) . The initial cloud formation did occur in SE 
Tennessee over the region of greatest topography; however, the 
cloud formation continued over much weaker relief as the day wore 
on. This suggests that topography and atmospheric moisture 
combined to produce the earliest clouds, but that sufficient 
instability and moisture was present to produce clouds over 
relatively weak relief later in the day. 

Other noted relationships included a positive correlation 
(-.75) linking MAMS channel 5 (similar to VISSR visible channel) 
visible brightness to warmer skin temperatures. This is 
presumably a relationship between surface vegetation type and skin 
temperature. A similar relationship was noted by Rabin et al. 
(1990) over the Oklahoma region. Also, in the cloud development 
region the skin temperature actually showed a small decrease with 
time (-.35 K) over a large area, possibly indicating cloud 
shadowing of the surface. 

4.3. Effect of MAMS Noise 

The effect of MAMS single sample noise on the derived 
products was evaluated by using single sample noise estimates for 
each MAMS channel in the <*]_ coefficient calculation. Since the 
coefficients represent a departure from a mean condition, a 
direct impact of noise can be determined by simply allowing the 
departure to be solely due to noise. An evaluation using several 


MAMS/VAS retrievals indicated that the MAMS noise contribution to 
the PW, TS , and LI products was about ,75 mm, .5 K, and .5 K, 
respectively. This represents about 15%, 5%, and 5% of the 
variation in each product respectively. Thus, real variation of 
PW, TS , and LI is easily observed using the MAMS/VAS retrieval 
methodology. 

4.4. Variability of Derived Products 

The spatial variability of the PW, LI, and TS fields was 
evaluated in a structure function calculation. The structure 
function has been described by several authors (e.g., Hillger and 
Vonder Haar 1979; Fuelberg and Meyer 1986). Under the assumptions 
of homogeneity and isotropy, structure becomes a function of only 
the scalar separation between data points. A graph of structure 
versus scalar separation depicts spatial gradients. Where the 
slope of the plot is greater, mean horizontal gradients are 
stronger over separation distance. 

The random noise of the derived products was evaluated by 
extrapolating the structure function to zero separation (Hillger 
ana Vonder Haar 1979) . With an assumption of no spatially 
correlated error, the structure function at zero separation 
becomes 

STR(O) =2 a 2 
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where a is the variance of the random noise. Structure function 
calculations were made at three locations using a 5 km by 5 km 
domain with a 100 meter separation interval (i.e., 50 X 50 pixel 
domain) . The structure between 100 m and 1 km separation 
distances was fit with a second order curve and extrapolated to 
zero separation distance. A second order curve fit was used 
because it best fit the shape of the structure data. Random noise 
estimates for the three locations of each product are shown in 
Table 4 along with the average of the three estimates for each 
product. The three estimates for each product show some variation 
making the random noise estimates somewhat uncertain. However, 
the averaged noise results for PW and LI are consistent with those 
obtained by using single sample MAMS noise estimates in the 
retrieval (Section 4.3), while the TS structure noise value is 
somewhat larger. The explanation for the TS discrepancy may be 
that the rapidly increasing TS structure at small separations 
creates uncertainty in the extrapolation to the noise value. 

PW, LI, and TS isotropic structure functions for a repeat 
overflight (1630 and 1730 UTC) are shown In Figs. 11-13, 
respectively. Flight track F-G was selected for the computation 
because it encompassed the moisture maximum identified in section 
4.2 as well as good visual variability in the LI and TS products. 
The structure was calculated out to 180 km separations with a 
separation interval of 300 meters (calculations have indicated 
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that using the data at higher resolution did not appreciably alter 
the results). Data pairs per separation interval exceeded 100,000 
for all separations. All three products demonstrate significant 
structure at meso 7 separations. The rapid structure increase at 
small separations reduces the possibility of noise dominating the 
structure. Hillger and Vender Haar (1988) noted that simulations 
showed gradient information was detectable at a level only 10 % 
above the extrapolated structure minimum. Calculations showed the 
ratio of the structure at the smallest separation (300 m) to that 
of the average extrapolated minimum is - 1 . 75 , 2 . 10 , and 4.60 
respectively for the PW, LI, and IS products. Thus real gradients 
are being observed at even the smallest separations. The moisture 
maximum in east central Tennessee is likely associated with the 
structure maximum at about 30 km separation distance. The LI 
structure shows a relatively weaker maximum at this separation, 
probably because the instability maximum in east central Tennessee 
is somewhat more spread out. shifting the variability to larger 
separations At large separations PW and LI structure both 
increase, indicating that atmospheric gradients are present at 
larger scales. This is consistent with the fact that the 
atmosphere varies significantly at larger scales (e.g. , across air 
mass boundaries), which has been demonstrated with structure 
functions bv Fueiberg and Meyer (1986) using VAS derived 
precipi table water products for separations out to 700 km. The 
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dominance of small scale gradients in the TS structure is not 
surprising given the highly variant nature of the thermal field in 
the COHMEX region. 

The one hour time evolution of the structure is also shown 
in Figs. 11-13. The atmospheric quantities, PW and LI, show 
changes in the shape of the structure curve with time, while the 
surface temperature, TS , shows a fairly uniform increase of 
structure with time. While it is difficult to explain changes at 
particular separations, the overall trends may be informative 
about basic atmospheric and surface responses to forcings. PW and 
LI structure slope decrease with time at small scale separations 
(up to 25 km), Indicating a smoothing of gradients with time. The 
structure maxima at about 30 km separation associated with the 
moisture and instability features in east central Tennessee are 
now weaker. Since PW and LI are affected by changes throughout 
much of the atmospheric column, there is uncertainty concerning 
the explanation of the reduced gradients with time. However, it 
can be reasonably speculated that mesoscale moisture expanded 
sDatiaily with time and smoothed the smaller scale gradients (as 
is indicated by the PW analysis of Fig. 6). Also, boundary layer 
vertical and horizontal mixing may be affecting the structure. 
Surface inhomogeneities cause spatial variations of sensible and 
latent heat fluxes which lead to horizontal variability of 
boundarv layer temperature and moisture. Mixing processes may 
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reduce this variability during the day as air is exchanged 
horizontally between thermals. Physical evidence of boundary 
layer mixing on June 19 is given by the onset of low level 
cumulus. This boundary layer mixing would be superimposed on 
other atmospheric processes such as mesoscale moisture advection. 
The net effect has been to reduce the smaller scale horizontal 
temperature and moisture gradients in this case. At larger 
separations (> 50 km) PW gradients become nearly constant 
(constant slope) . This also indicates gradient smoothing with 
time, however caution is advised as the separation distance 
approaches half of the domain size (-200 km). Features with 
wavelengths approaching 100 km and larger may be aliased in the 
structure data. Larger separation LI structure shows gradients 
weakening with time. This reflects the rapid destabilization that 
took place in western Tennessee over the one hour interval Fig. 
7). The result of this was a smoothing of larger scale gradients. 
This rapid destabilization may have been caused by larger sensible 
heat t luxes m western Tennessee. The LI structure change may 
also reflect temperature cnanges In the mid- troposphere with time 
although given the synoptic conditions, such changes would be 
expected to be small over a one hour Interval. The fairly uniform 
increase of TS structure can be attributed to the inabilitv of the 
surface to quickly mix thermal gradients ; warm surfaces heat up at 
a faster rate than cool surfaces (which are presumably cooler 
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because of considerations involving vegetation type, soil moisture 
content, etc.). The increased surface thermal gradients are seen 
as increased structure. It is expected that these thermal 
gradients decrease when the solar input is removed. A slight 
indication that differential surface heating occurred over larger 
scale surface inhomogeneities is shown by small maximums of 
structure at about 70 km and 90 km separation distance. 

In the context of these physical explanations, the evolution 
of the structure functions demonstrate reasonable atmospheric and 
surface trends. While the PW and LI structure functions might be 
expected to show more continuity of features with time, these 
results indicate the atmosphere's ability to evolve rapidly at 
small scales. Such rapid mesoscale evolution has been pointed out 
by Fuelberg and Meyer (1986; in both satellite as well as 
radiosonde PW. However, while most of the structure trend is 
attributable to physical processes, influences such as limb 
effects, undetected cloud, noise charac teris tics , and aliasing at 
large separations are also affecting the structure. 

A comparison of MAMS/VAS and VAS-oniy structure is shown In 
Fig. 14 for the PW product. The VAS PW data is taken from special 
data set soundings produced for the COHMEX scientific community 
(Veldon et al. 1988). The primary difference between the MAMS/VAS 
and VAS COHMEX special data set is the spatial resolution (i.e., 

100 meter vs . 8 km). It is immediately obvious that the 
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instrument resolution impact on the PW product is quite large. 

The MAMS/VAS product demonstrates larger gradients at relatively 
small separations that the VAS COHMEX data set does not capture. 

At separations less than 50 km the MAMS/VAS PW structure is 
roughl" twice that of the VAS structure. In addition, larger 
sample sizes per separation interval should make the MAMS/VAS PW 
structure less susceptible to noise. Thus, high resolution data 
provides both a capability to identify small scale atmospheric and 
surface variability, and reduces the effect of noise on the data 
through larger sample sizes for analysis at given spatial scales. 

A similar comparison between MAMS /VAS LI structure and VAS -only LI 
structure (not shown) showed much more similar structure functions 
at small separations, probably because the VAS data contributes 
more to the MAMS/VAS LI product than to the PW product; had a 
comparison been mace tor TS , it is expected that results would 
nave een similar to those of the PW product structure comparison. 
*his comparison, while representing only a single case, clearly 
suggests that graaient information at small scales is not being 
captured by relatively coarser resolution instruments such as VAS. 

f . CONCLUSIONS 

A technique for combining 100 meter resolution MAMS and 8 km 
resolution VAS radiances to produce high resolution maps of 
precipitable water, lifted index and skin temperature has been 


24 


demonstrated. Each quantity exhibits real spatial and temporal 
variation az the meso 7 and larger scales. In particular, a 
moisture maximum in east central Tennessee during the COHMEX day 
of June 19, 1986 has been identified, as well as a heat island at 
Nashville, Tennessee. The June 19, 1986 case study of these 
quantities and their evolution in the context of a developing low 
level cumulus field Indicated that returning moisture on the 
backside of high pressure had an apparent influence on the 
development and existence of low level cumulus in the region. 

This suggests that atmospheric variability is influential in cloud 
development. A lifted index quantity also demonstrated some 
ability to identify the location of the cumulus development. 
Topography is also indicated as an important contributor to the 
initial cloud development. Skin temperature data was more 
-unbiguous, with clouds forming over a relatively warm region while 
not forming over other comparably warm regions. A possible effect 
of cloud shadowing was also noted in the skin temperature data. 
Structure function analyses were used to present a 
quantitative measure of the spatial variability of each product. 
Real atmospheric gradients of moisture and stability have been 
identified at scales down to 300 ra. Surface variability 
(dependent on local land use) was found to be dominant at small 
scales. A comparison between VAS 3 km resolution and the MAMS/VAS 
100 meter resolution precipitable water structure demonstrated the 
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additional variability identified by using high resolution data. 
These results indicate chat a geostationary^ sounding capability 
with better than 1 km resolution is needed if most of the raeso 7 
scale variability in the atmosphere (and on the surface) is to be 
monitored over short temporal intervals (< 1 hr). 
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Table 1. 
channel 

MAMS spectral 
bandwidth ($ 

channels . 

50% resDonse 

l a 

0.42 - 

0 . 

.45 

2 

0.45 - 

0 . 

52 

3 

0.52 - 

0 . 

60 

4 

0.57 - 

0 . 

67 

5 

0.60 - 

0 . 

73 

6 

0.65 - 

0 . 

83 

7 

0.72 - 

0 . 

99 

8 

0.83 - 

1 . 

05 

9 b 

3.47 - 

3. 

86 

10 

6.28 - 

6. 

98 

11 

10.55 - 

12 

.24 

12 

12.32 - 

12 

.71 


a Channel I is not available when 10 bit data is 
collected . 

k Only three infrared bands are available at once. The 3/xra 
and 6/im bands are interchangeable. 


Table 2. VAS spectral channels. 


channel 


bandwidth (& 50% response) 


1 

14.63 - 

14.84 

2 

14,31 - 

14.65 

3 

14.09 - 

14.42 

4 

13.82 - 

14.21 

5 

13.16 - 

13.51 

6 

4.48 - 

4.57 

7 

12.30 - 

12.82 

8 

10.36 - 

12.12 

9 

7.15 - 

7.37 

10 

6.40 - 

7.08 

11 

4.41 - 

4.48 

12 

3.34 - 

4.06 



Table 3. MAMS flight track data for June 19, 1986. 


Time Coincident 


Track 

(GMTl 

VAS time 

Remarks 

B-C 

1521 - 

1545 

— 

no VAS overlap 

D-E 

1553 - 

1610 

1600 

clear skies 

F-G 

1621 - 

1643 

1600 

clear skies 

H-I 

1650 - 

1713 

1600 

some cirrus 

F-G 

1720 - 

1743 

1730 

repeat overpass 





cumulus in 





eastern quarter 

H-I 

1750 - 

1814 

1730 

repeat overpass 


some cirrus 


Table 4. Precipitable water (PW) , lifted index (LI), and skin 
temperature (TS) random noise estimates using structure functions. 
Estimates are presented for each of three locations and an average 


of the three locations for each product. 
Product Extranolated noise 

Average 

PW! 

0 . 

z 2 

;n 



pw 2 

0 . 

65 

mm 

0.74 

mm 

PW 3 

0 . 


mm 



LIi 

0 . 

- o 

K 



LI 2 

0 . 

:6 

K 

0.23 

K 

LI 3 

0 . 

20 

K 



TSi 

2 . 

— 

:< 



ts 2 

0 . 

£ 5 

K 

1.55 

K 

TSt 

1 _ . 

- 5 

K 
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FIGURE CAPTIONS 


Figure 1 - MAMS instrument scanning geometry. At 20 km altitude, a 
MAMS instataneous field of view (ifov) is about 100 meters on a 
side. Insert demonstrates overlap of ifovs from scan line to scan 
line. 

Figure 2 - Normalized spectral response functions of MAMS and VAS 
6pm (top), 11pm (middle) and 12pm (bottom) channels. 

Figure 3 - Surface pressure and temperature (top) analysis and 500 
mb height and temperature analysis (bottom) for June 19, 1986 at 
1200 UTC. Temperatures in 3 C. pressure in millibars, heights in 
decameters (from Jedlovec 1990), 

Figure 4 - GOES 1600 GMT visible imagery of the COHMEX region on 
June 19, 1986 with the COHMEX Beta Network radiosonde locations 
(numbered) and the MAMS flight tracks (lettered). Note the 
generally clear conditions in the MAMS flight region. 

Figure 5 - Guess (solid) and departure (dashed) for retrieved 
atmospheric profiles of temperature (T) and moisture (T^) for VAS 
(top) and MAMS/VAS (bottom) physical retrievals. The guess 
profile used in the MAMS/VAS retrieval is the retrieved profile in 
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the top panel, thus deviations from the guess in the MAMS/VAS 
retrieval are due to the MAMS data. 

Figure 6 - MAMS/VAS precipitable water product (mm) with analysis 
for flight tracks D-E, F-G, and H-I (top) and for repeat overpass 
of flight tracks F-G and H-I (bottom) one hour later. Vhite areas 
are cloud cover. See Fig. 4 for Location of flight tracks. 
Location of Nashville. 7N (BNA) is shown. 

Figure 7 - MAMS/VAS lifted index^. stability product (deg) with 
analysis for sa a region as Fig. 6. Positive values are most 
s tab 1 e . 

Figure 3 - MAMS/VAS skin temperature product (deg) with analysis 
for same region as Fig. 6. 

Figure 9 - COHMEX Beta Network radiosonde precipitable water (top) 
and lifted index (bottom) values with analysis for 1300 GMT on 
June 19, 1986. Radiosonde data in large numerals. 

Figure 10 - GOES 1900 GMT visible imagery for June 19, 1986 over 
the COHMEX region with MAMS flight tracks and COHMEX Beta Network 
radiosonde locations overlain. Note the cumulus cloud development 
that has taken place from 1600 GMT (Fig. 4) to 1900 GMT. 


n c 
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Figure LI - MAMS/VAS precipitabLe water isotropic structure 
function for flight track F-G initial (1630 UTC) and repeated 
(1730 UTC) overpasses of June 19, 1986. 

Figure 12 - MAMS/VAS lifted index isotropic structure function for 
flight track F-G initial (1630 UTC) and repeated (1730 UTC) 
overpasses of June 19, 1986. 

Figure 13 - MAMS/VAS skin temperature isotropic structure function 
for flight track F-G initial (1630 UTC) and repeated (1730 UTC) 
overpasses of June 19, 1986. 

Figure 14 - MAMS/VAS and VAS-oniy precipitable water isotropic 
structure functions over the CGKMEX region on June 19, 1986. Note 
the MAMS/VAS data shows much greater structure at small 
separations . 




ORIGINAL PAGE !S 
OF POOR QUALITY 


ORIGINAL PAGE IS 

OF POOR QUALITY 






PRESSURE (mb) 










25 50 75 

100 125 150 

175 200 


SEPRRRTION 

DISTRNCE 

(KM) 

— 


Fig. 11 









SEPRRfiTION DISTANCE ( KM) 




= n 

w sj 


/3 


100 125 150 


j | — i — i — i — i | 

175 200 


SEPRRRT lOM DISTANCE (KM) 


13 






HIGH RESOLUTION ATMOSPHERIC AND SURFACE VARIABILITY 
FROM COMBINED MAMS AND VAS RADIANCES 



Christopher C. Moeller 1 
W. Paul Menzel 2 
Kathleen I. Strabala 1 

Cooperative Institute for Meteorological Satellite Studies 
2 NOAA/NE3DIS Advanced Satellite Products Project 
1225 West Dayton Street 
Madison, Wisconsin 53706 



I. INTRODUCTION 

Measuring of the variation of atmospheric and surface quantities such as moisture and stability at small scales is 
an important step to monitoring, understanding and forecasting mesoscaie processes. ""'While few ground based observing 
systems provide the high resolution coverage with sufficient areal extent necessary to study mesoscaie variation, satellite 
observing systems have proven to be extremely useful in this endeavor. HFor instance, moisture and stability products from 
the geostationary platform VISSR Atmospheric Sounder (VAS) exhibit interesting detail at the meso B (20-200 km) level 
(Chesters et al., 1986; Smith et al., 1985). However, variation at the meso C scale (2-20 km) and at the lower end of the 
meso B scale ( about 20-100 km) suffer from the effects of noise. Noise can be greatly reduced by using a high resolution 
data collection system. The Multispectral Atmospheric Mapping Sensor (MAMS), with its 100m resolution over a 36km 
swath, is well suited to provide high resolution mesoscaie information about the atmosphere and surface of the earth. A 
more complete depiction of the atmospheric state at the smaller scales is possible through the combination of the high 
horizontal resolution MAMS data with the vertical sounding data of the VAS. 

This paper demonstrates a procedure for combining MAMS and VAS data in a physical retrieval to produce high 
resolution derived products of precipitable water (PW), fiftetTindex (LI), and skin temperature (TS). The variability of 
these products is evaluated in a structure function analysis similar to that of Hillger and Vonder Haar (1979). Results for 
19 June 1986 from the Cooperative Huntsville Meteorological Experiment (COHMEX) are presented. 


II. ALGORITHM 

The MAMS/VAS combined retrieval is possible because the MAMS and VAS have spectrally similar 6, 11, and 12 
micron spectral bands (Menzel et al., 19861. The MAMS high resolution data can be incorporated into the VAS retrieval 
process by replacing the appropriate VAS data with the MAMS data. The resulting retrieval benefits from the MAMS high 
resolution information as well as from the multispectral sounding information from VAS. More complete descriptions of the 
MAMS can be found in Jedlovec et al. (1989) and of the VAS in Montgomery and Ucceilini (1985). 

A physical retrieval process for producing derived imagery (Smith et al., 19851 i« used, which entails a one step 
simultaneous solution (Smith and Woolf, 1984) to produce a retrieval for VAS data covering an area of 5 x 5 fields of view 
;FOV). The initial full VAS retrieval at coarse resolution is then used as the first guess for combined MAMS/VAS 
retrievals (where the 6, 11, and 12 micron MAMS data has replaced the equivalent VAS data) at 100m resolution. Biases 
with respect to the VAS radiances are removed from the MAMS radiances. Cloud checks with the infrared window are 
performed at each MAMS FOV; no MAMS - VAS retrieval is produced in a cloudy FOV 


III. EVALUATION OF TECHNIQUE 

A. Processing of 19 June 1986 Data 

On June 19, dry Canadian high pressure which had moved over the COHMEX region on the 18th began to slide 
off to the east, reducing the subsidence over the region as well as increasing the low level moisture, especially in the 
southern part of the COHMEX region. Although no significant convection occurred, sufficient moisture and surface heating 
were present to produce low level cumulus at the top of the boundary layer beginning about 1700 GMT. During the day 
this cumulus field expanded horizontally to cover much of east central Tennessee southward into northern Alabama. 

During the afternoon, a few of these cumulus clouds did grow suffTciently to produce brief light showers over the region. 

The MAMS was flown over the COHMEX region from about 1530 GMT to 1800 GMT on this day. Figure 1 shows 
the COHMEX region with the straight line flight track* that were flown. VAS dwell sounding pain were collected at 1600 
GMT and 1730 GMT. 
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B. _ “ Analysis of D^ved Pfodu<J* f ^ 

The PW, LI, and TS derived product! are shown in Fig. 3. Factors common to all 3 products include a varying 
degree of zenith angle effect, a weak left hand side bias (as the aircraft flew) and thin cirrus contamination in the NW 
sector. 

The derived PW image shows generally higher values to the east and south and generally drier conditions to the 
north and west. A surge of moisture from the south can be seen penetrating into east central Tennessee. The gradient to 
the west and north shows a decrease of about 6 mm over a distance of about 100 km with a more concentrated gradient of 
about 3-4 mm over about 20 km to the northwest of the moisture maximum. This moisture surge also shows up in the more 
coarsely spaced COHMEX Beta Network radiosondes (not shown), though without as much certainty about its location. 

The radiosonde PW values show fairly good relative agreement, while absolute comparisons show a 2 to 4 mm high bias in 
the MAMS/VAS PW. The moist tongue is also identifiable in the repeated flight tracks one hour later, which reveal the 
time evolution of the feature. Interestingly, an onset of low level cumulus is seen in the repeat overpass in the moist tongue. 
In the next couple hours this field of cumulus expanded westward to the center of flight track F-G, in good visual 
correlation with the moist tongue feature. 

The LI image shows that the most unstable portion of the flight exists over east-central Tennessee. This is in 
good agreement with the relative variation of the COHMEX 1800 GMT radiosondes, but a 3 to 4 unit bias towards unstable 
values is present in the MAMS/VAS LI. The general location of the low level cumulus field that developed agrees with the 
relatively unstable area in east central Tennessee. The repeated flight tracks indicate that much of the area is becoming less 
stable with time. This would be expected as solar heating warms the surface over the one hour interval. However, it 
appears the most rapid destabilisation is taking place in the western portion. This suggests cloud development further west 
than actually occurred. However, the most unstable conditions (disregarding the thin cirrus region) prevail in the 
southeastern portion of the repeated flight tracks, in the location of the moist tongue. While sufficient low level instability 
(from surface heating) is apparently present as indicated by the expansive region of unstable LI values, cloud development 
is restricted to the moisture concentrated in east central Tennessee. 

The skin temperature (TS) image shows considerable high resolution variability. The various land uses in the 
COHMEX region (e.g., forest, agriculture, urban) contribute strongly towards this variability. In particular, a definite heat 
island is observed at the location of Nashville, TN (BNA) with skin temperatures 5 to 10 K warmer than the surrounding 
region. Time change of TS indicates good continuity of the thermal features. Indications are that individual FOV skin 
temperatures increase by as little as 0 K to as much as 5 K over the one hour time interval. This is a result of the thermal 
properties of the various land types in the scene. Over a selected larger region (about 20 X 25 km) the average skin 
temperature increase was 2.75K. Such skin temperature monitoring is useful in radiation budget studies of the surface of 
the earth. Very warm TS values exist throughout the imagery and don't seem to be correlated to the low level cloud 
development. This once again points out the apparent importance of the moisture on this day to cloud development. 

C. Effect of MAMS Noise 

The effect of sample noise on the derived product* was evaluated by applying uniform single sample noise 
estimates to the retrieval coefficients. Since the retrieval coefficients were generated to represent a departure from a mean 
condition, the direct impact of noise can be observed by simply allowing the departure to be solely due to noise. For two 
retrievals in flight track F-G, the noise contribution to the PW, TS, and LI products was found to be about .75mm, .5K, 
and ,5K respectively. This represents about 15%, 5%, and 5% of the deviation in each product respectively. Thus, real 
gradients of PW, TS, and LI are easily observed by the MAMS/VAS retrieval methodology. 

D. Variability of Derived Products 

The spatial variability of the PW, LI, and TS fields was evaluated in a structure function calculation. This 
provides insight into the scales of natural variability that each derived product exhibits (Hillger and Vonder Haar, 1979). 

In a graph of structure versus separation distances, the slope of the plot is greater where mean horizontal gradients are 
greater. 

The results of structure function calculations for PW, LI, and TS are shown in Fig. 2. Flight track F-G(l) was 
selected for the computation because it encompassed the moisture tongue identified in section III.B as well as good visual 
variability in the LI and TS products. Also, there were no clouds in this flight track. To reduce computation time of the 
structure calculation, the domain was subsampled at every 5th line and element to produce an effective resolution of 500 m; 
this left a domain of about 500 lines by 72 element* (250 km by 36 km). A separation interval of 1500m was used. Because 
the dimensions of the domain are greater in the along track direction, isotropic structures beyond 36 km will become 
increasingly biased in the along track direction. 


The PW structure function display* a generally smooth increase from 1 km up to about 30 km separation. In this 
interval, PW increaaes from about 1.5 mm up to about 3 4 mm (including noite contributions). Small amplitude 
fluctuations of the structure in this interval are less than the noise level and cannot be regarded as significant. From SO km 
to about 120 km separations the structure function flattens out, indicating that the gradient of PW is not increasing in this 
range. Beyond 120 km the PW structure again begins to increase, as PW varies significantly at larger scales (e g., across 
frontal boundaries, etc.). This also indicates that the flat structure function between 30 and 120 km is not a systematic bias 
of the calculation, but rather represents a true lack of atmospheric variation at that scale. The indication here is that 
significant atmospheric variation is taking place at small scales. While this makes no statement about the meteorological 
importance of this variability, it does demonstrate its existence. 

The LI structure shows a step-like pattern with strong increases in structure from 1-25 km, 60-75 km, and 100- 
120 km. This behavior is significantly different from that oTthe PW structure. Instead of leveling off as the PW structure 
does, the LI structure continues to increase as separations increase. Because the LI product is primarily dependent on the 
coarser resolution VAS data, the LI structure may be somewhat limited in its ability to depict the true high resolution 
atmospheric variability of LI. However, it is still indicated that small scale variability is certainly present and measurable. 

The TS structure function exhibits behavior both similar to and different from the PW and LI structure functions. 
Like the PW and LI structure, the TS structure increases rapidly at small scales with most of the increase contained in the 
1-15 km separation range. This rapid increase in TS structure at small separations is not surprising given the different land 
uses in the region. The structure flattens out over the remainder of the separation range. This is different from both PW 
and LI structure, and indicates a strong preference for high small scale variability of skin temperature in the COHMEX 
region. It is noteworthy that the TS variation of 5.6K at 1500m separation is significantly greater than the noise 
contribution. The implication here is that small scale variation is dominant and measurable. 

IV. - ^CONCLUSIONS 

This pa^er has presented the results of merging the 6 km resolution sounding data from the VAS with the 100m 
resolution moisture data from the MAMS in the physical retrieval of precipitable water, atmospheric stability (lifted index), 
and skin temperature. The VAS provides the vertical information of the atmospheric state and the MAMS indicates the 
horizontal variability. Structure function analysis of the resulting derived product fields indicates that both atmospheric 
PW and LI products have significant variability at less than 30 km spacing. Inspection of the surface skin temperature 
reveals structure at less than 15 km due to differences in local land use. The implications for mesoscale monitoring of the 
atmospheric state are that definition of small scale variations require horizontal resolutions of better than 5 km. 
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SEPRRRTION DISTANCE (KM) Figure 1 - MAMS flight path on June 

19, 1986 over the COHMEX region. 
Letters identify beginning and end of 

Figure 2 - Isotropic structure straight line flight tracks. The two 

functions of PW, LI and TS for June northernmost flight tracks (F-G,H-I) 
19, 1986 flight track F-G(l). were repeated. Dashed line indicates 

NW boundary of cumulus cloud 
development at about 19 GMT. 




Figure 3 - Fields of precipitable water (PW), lifted index (LI) and skin temperature 
(TS) for the entire flight of June 19, 1986. The data has been remapped into a 
Mercator projection with a degradation of resolution down to .6 km so that 
adjacent flight tracks could be shown in a single image format. Each product 
consists of a composite of 3 flight tracks (D-E, F-G, H-I) and a composite of the 2 
repeated flight tracks (F-G, H-I). Thus the northernmost two flight tracks of the 3 
flight track composite are collocated wi th the composite of the repeated flight 
tracks. White pixels arc cloud cover. Thin white vertical dashes in the imagery 
arc a result of the remapping process. 
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1 . INTRODUCTION 

The Multispectral Atmospheric Mapping Sensor 
(MAMS) is a new instrument designed to produce high 
resolution Imagery in eight visible and three 
Infrared spectral band s f rom an aircraft platform. 
Thermal emission from the earth's surface, clouds, 
and atmospheric water vapor is measured at 12.3, 

11.2, and 6.5 microns at up to 50 meter horizontal 
resolution. Similar infrared spectral bands are part 
of the VISSR Atmospheric Sounder (VAS) which is 
capable of 6.9 km horizontal resolution from its 
geostationary orbit. Sequences of 6.7 micron VAS 
Images have been used to track water vapor structures 
and to infer atmospheric motion vectors (Stewart et 
al., 1985). They have also been used to investigate 
smail scale atmospheric moisture variability 
(Jedlovec, 1984). The high resolution MAMS data give 
credence to the small scale variations in the VAS 
water vapor imagery. 

The purpose of this paper is two-fold. We 
shall (1) describe the MAMS instrument and the 
radiance characteristics, and (2) intercompare some 
VAS and MAMS radiances and images. 


2. MAMS INSTRUMENT CHARACTERISTICS 

The MAMS is a modification of NASA's Airborne 
Thematic Mapper (Daedalus, 1982). The MAMS has two 
major components, the scan head and the spectrometer 

(Fig. 1 ) . The scan head consists of tha primary_ 

collecting telescope, a rotating scan mirror, a motor 
encoder assembly, and two controlled thermal source**. 
The field of view is defined bv the aperture which is 
available in two different sizes, 2.5 mrad and 5.0 
mrad. The spectrometer consists of optical elements 
which spectrallv separate the polychromatic input 
energy, lenses which focus the separated energy onto 
sensors, detectors which convert optical energy to an 
electrical signal, and pre— amplifiers which condition 
the signal. The energy collected by the spectrometer 
is separated into four distinct optical paths using 
dichroic filters; three paths tor the infrared data 
(one for each spectral band) and one for the visible 


and near infrared data. The spectrometer 
detectors are silicon for the visible radiation 
and mercury cadmium telluride for the thermal 
infrared radiation. The visible detector is an 
eight element arrav which converts oprical 
energy to electrical energy. The spectral 
response of each band is determined by the 
dispersion characteristics of the prism placed 
in the optical path and the location of the 
arrav in the dispersed beam. The thermal 
infrared detectors are each housed in a vacuum 
dewar which contains a cooled long wave filter 
which defines the particular spectral band. 

Each infrared detector has an associated pre- 
amplifier to process the electrical signal. The 
spectral bands are summarized in Table 1. 

The spectral response functions for the three 
infrared charnels of the MAMS (Fig. 2) were 
specified to be similar to the 6.7, 11.2, and 
12.7 micron channels on VAS (also shown in Fig. 
2 ). 

The MAMS produces high resolution imagery 
in these visible and infrared spectral bands 
from a high altitude aircraft at 20 km as shown 
in Fig. 3. The horizontal ground resolution, 

GR, of each individual field of view is 100 

Table 1 


Rangle of MAMS Spectral Bands 






8 visible 

3 Infrared 


(microns) 

(microns) 


.42- 

.45 

6.20- 6.85 (upper level 

water vapor) 

.45- 

.52 

10.32-12.02 (window > 


.52- 

.60 

12.20-12.56 (lower level 

water vapor) 

.60- 

.62 



.63- 

.69 



.69- 

.75 



.76- 

.90 



.90-1 

.05 
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SPECTRAL RESPONSE (PERCENT) 




WAVELENGTH (MICRONS) 


Fig, r, Spectral response curves indicating the 
<*erMtlvity or the MAMS and VAS instruments In 
the 12 micron, 11 micron, and b micron regions. 


r-eters (50 meters if the smaller 2.5 tnrad optics 
used). The total field of view for one scan 
*s 8 6° which covers the width, w, of the entire 
cress path of roughly 40 km. The scan rate Is 
6.25 revolutions per second (12,50 when the 2.5 
mrad optics are used; and produces 67Z overlap, 
I-P, of each scan line with the previous one for 
romiral aircraft speea (740 km/hr). On each 
scan line 716 samples are taken so that there is 
52T overlap of neighboring pixels along a scan 
line (4% when the 2.5 mrad optics are used). 
These overlaps can be used to good advantage to 
reduce the noise in the scene by spatial 
averaging. 

Radiometric calibration of the infrared 
data is accomplished bv viewing two blackbodv 
sources and by assuming the detector response is 
T inear with respect to radiance (a fractional 
nonlinearity of less than 5x10 14 is observed). 
Calibration occurs everv scan line by viewing a 
vara and then a cold blackbodv of known 
temperature. For a given spectral band, the 
raaiance is determined from the temperature 
through the convolution of the normalized 
speerrai response ard the Planck function. The 
calibrated radiances are transformed back into 
temperatures using the inverse Planck relation 
with the nominal spectra) hand center wave- 
1 ength . 

Prior co flight, the gain and offset for 
each spectral band are adlusted electronically 
in order to minimize rhe likelihood of satura- 
tion over hoc scenes and tc maximize the effec- 
tive operating range. The 6.5 micron band is 
adjusted to cover 210 to 280*K, while the 11.2 
and 12.3 micron bands operate from 220 to 340°K. 
An eight bit digitizer is used so chat for the 
6.5 micror band the gain is approximately . 3°K 
per count, and for the remaining surface viewing 
bands it is approximately . 5 # K per count. 
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Fig. 3a. A simplified schematic of the inflight 
configuration of the optical elements of the 
HAMS . 



Fig. 3b. Scanning geometry for the MAMS instru- 
ment onboard a NASA U2/ER2 aircraft. Lower 
insert shows position of the center of the nadir 
ifov for two adlacenr scans. 

Truncation errors are assured to be roughiv half 
these values. 

The infrared detectors on HAMS exhibit 
some low frequency noise so that samples taken 
at 52 microsecond intervals for the soin rate of 
6.25 revolutions per secona (26 microsecond 
sampling occurs at 12.5 rpfl) experience noise 
contributions with significant correlation. An 
analysis of covariance (Menzei, I960) or the 
noise with the 5 mrad aperture reveals that 
after averaging ten consecutive samples noise is 
reduced by a factor of .54, . :2, ana .37, 
respectively, in the three infrared spectral 
bands in order of increasmz wavelength __ 

(uncorreiated noise would have a factor of 1/*10 
or . 32) . 

Inflight single sample noise vaJues were 
determined from multiple samples over a uniform 
portion of the ocean. For a sample sire cf 100, 
the single sample noise values with the 5 mrad 


aperture were found to be .2, . !, and .4°K, 
respectively , for the spectral bands in order of 
increasing wavelength (for the 2.5 mrad aperture 
the values were .5, .2, and .9°K). 


3. DATA I NTEP COMPARISONS 

Radiances for the MAMS and VAS spectral 
responses were simulated for each band and 
sensor in a radiative transfer calculation. A 
representative atmospheric transmittance was 
computed using a bard model and radiosonde 
temperature and moisture profiles. Table 2 
shows that the MAMS detects radiances within one 
to two mW/ster/nr/cm” 1 of the VAS detected 
radiances. The MAMS spectral responses are 
shifted to shorter wavelengths (as shown in Fig. 
2), hence less radiance is detected. This 
translates into brightness temperature 
differences of Jess than a degree Kelvin. 

Table 2 

Intercomparison of Simulated Data 



6 micron 

11 micron 

12 micron 

R(MAMS) 

7.86 

96.5 

109.5 

T(MAMS) 

254.8 

287.3 

286.0 

R(VAS) 

8.73 

97.5 

Ill.fi 

T(VAS) 

253.9 

287.3 

285.4 — 

^"'J4MC_U4C 

0. 9°K 

0.0 

0.6 


radiances in mW/ster/m 2 /cm~ 1 , temperatures in °K 


Radiances from the MAMS and VAS were 
collected simultaneously over a relatively 
isothermal area of ocean off the California 
coast on 18 May 1985. Table 3 presents the 
comparison. After compensating for the inherent 
spectral differences of the two instruments, the 
MAMS data is within .5, 2.5, and 1.5°K of the 
VAS data, respectively, for the spectra] bands 
in order of increasing wavelength. The observed 
split window channel differences (11 micron less 
12 micron) are comparable, 2.5°K for MAMS and 
2.I 3 K for VAS and thus showing similar low level 
moisture concentrations. The observed water 
vapor channel data (6 micron) show MAMS 
detectine upper level moisture concentrations 
verv much like VAS, with MAMS showing somewhat 
less attenuation than VAS. These comparisons 

Table 3 

Intercomparison of Observed Data 



6 micron 

11 micron 

12 micron 

= 

R(MAMS) 

6.01 

99.9 

111. I 

- _ 

T(MAMS) 

247. 1 

oo 

c/1 

287.0 


R(VAS) 

fi.80 

97.0 

111.1 


T (VAS) 

246.5 

287.1 

285.0 

— 

at mams-vas 

. 6"K 

2.4 

2.0 

■ 

radiances 

in mW/ster/m 2 

/cm 1 , temperatures in "K 
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are quite reasonable, since the spatial resolu- 
tion of each sensor is very different (VAS at 7 
km and MAMS at .05 km at nadir ) , the navigation 
of the MAMS data is from visual siting of 
landmarks only, and non-nadir viewing of the 
MAMS has not yet been taken into account. 

Imagery from the MAMS and VAS were 
compared from flights over Oklahoma and Kansas 
on 12 May 1985. The goal of the intercoraparison 
is to verify small scale features observed in 
the VAS water vapor images with the higher 
resolution MAMS data. These water vapor inhomo- 
gereities in the VAS data have been used to 
obtain motion vectors when several images are 
considered in sequence. While the imagerv is 
often too amorphous to track winds reliably with 
a correlation technique, it has produced good 
wind fields by single pixel tracking (Stewart et 
al. , 1985). The MAMS data tends to verify that 
snail scale gradients (down to several single 
pixels) in the VAS data are real. Fig. 4a shows 
the VAS 6.7 micron images with the MAMS flight 
track superimposed while Fig. 4b shows the 6.3 
micron MAMS data (at 400m resolution for dis- 
play) gathered in flight from C to H. These 
images were enhanced to delineate the moisture 
features. Brighter (colder) regions represent 
clouds or relatively high water vapor content 
and darker (warmer) areas represent drier air ih~ 
the middle layer of the atmosphere. The MAMS 

images are not limb corrected and, therefore, 

Images appear brighter (colder) toward the 
edges. The moisture features compare very 
nicely; starting at G, MAMS shows a dry tongue 
giving way to p moist ridge which becomes less 
moist at H. Fig. 5a shows single pixel resolu- 
tion of the VAS image covering the MAMS crack 
from E to F of Fig. 4a; Fig. 5b displays the 
MAMS data covering the boxed region of Fig. 5a. 
The temperatures from brightest to darkest 
pixels within the boxed-in area differ by 5®K 
indicating that these fluctuations are not 
instrument noise, but real atmospheric features. 
Six of the VAS features are corroborated nicely 
by the MAMS image. This lends credence to the 
assertion that the "'AS distinguishes moisture 
features which are small enough to represent a 
wind vector, yet large enough to exist in 
several successive Images. The features we are 
examining were discernible for at least two 
hours in the VAS images. 


The MAMS has been developed to combine 
high resolution visible and Infrared measure- 
ments to study small scale atmospheric moisture 
variability, to monitor and classify clouds, and 
to investigate the role of surface characteris- 
tics in the production of clouds, precipitation, 
and severe storms. As part of the Cooperative 
Huntsville Meteorological Experiment (COHMEX) , 
the MAMS will be flying this summer to gather 
data for these studies. 
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4. CONCLUSION 

The MAMS measures thermal emission from 
the earth atmosphere systen in e ight visible and 
three infrared soeccral banas (12.3, 11.2, and 
6.5 microns j at up to 50m horizontal resolution. 
These infrared bands are very similar to three 
of the VAS infrarea spectral bands. The instru- 
ment was flown aboard a NASA U-2 in Mav-June 
1985 ar.d corresponding VAS data was collected. 
MAMS radiometric performance is comparable to 
that of VAS. The M-AMS exhibits somewhat Jess 
attenuation from water vapor than VAS because 
its spectral bands are shifted to shorter wave- 
lengths away from the absorption band center. 
Intercomparisons of these images give credence 
to small scale variations in the VAS water vapor 
imagery. Attempts are underway to identify the 
scales of which these features must be tracked 
to depict the preconveccive environment. 


Ill 



U 40 
<0 O' 


00 >. 

O' 1C 


JK 40 
* — 

1 c 

ao o 


u 

o o 


TJ L3 
01 

to e 
c o 
a u 
6 
T* 

V XL 
CL DC 
3 

« *”H 


£ 4J 
00 £ 


3E E 


OJ 

-3 T3 


-c c 

*-> 1C 
f-» £ 

3 c 

V 

01 

00 c 

<g o 


c 

a. en 


I 2 

oc 
on c 
< 

> "C 

c 

■o o 
4) a. 

O CO 
C 4> 

flj Ini 
XL U 

e o 

UJ o 


0^ OG 

u. u. 


ORIGINAL PAGE IS— 
OF POOR QUALITY = 








OF POOR QUALITY 



. S7A113 97 


DETECTION OF MOUNTAIN- INDUCED MESOSCALE 


WITH ETCH RESOLUTION >«ISTURZ IMAGERY 


Gary J. Jedlovec 
CXMSS, Univ. of Vise, 
Had Ison, Vise* 


V. Paul Henzel 
NOAA/ASPP 
Madison, Wise. 


Gregory S. Wilson 
NASA/MS FC 
Huntsville, Al. 


Robert J* Atkinson 
General Electric Co. 
Huntsville, Al. 


1 . INTRODUCTION 

Radiosonde observations near the surface 
and aloft often detect pressure, temperature, and 
wind patterns associated with large scale atmospner' 
ic wave features which govern our daily weather 
conditions. The detection of small scale waves 
such as gravity and bouyancy oscillations however; 
usually require special observations to measure the 
perturbations directly. These waves need not al-— ^ 
ways be measured or sensed directly and their ob- 
servable signatures (i.e., clouds) can be detected 
in visible photography or satellite imagery. 

The detection of mountain and lee wave 
structures is often done by observing cloud sig- 
natures in regions where vertical oscillations are 
forced by terrain. Satellite observations or these 
small scale clouas are possible with high resolution 
visible sensors. Direct observation of wind ana 
thermodynamic parameters is achieved by instrument- 
ed aircraft, Ildar, and radar measurements. From 
these ooservations, quantitative information is 
available aDout mountain and lee waves, their form- 
ation and three dimensional structure. However, 
because a moist environment 13'fidt hieessirf , aanv 
wave features may go unaetected by these convention- 
al measurements because of the absence of clouds. 
During one of the flights of the Multispec- 
tral Atmospheric Happing Sensor (HAMS) on a NASA 
U2 high altitude aircraft in 1985, wave features 
were found in water vapor channel imagery in close 
proximity to the Sierra Nevada mountains of Cali- 
fornia (Fig. 1). Visible channel imagery from the 
flight was examined for cloud features but none 
were found. This paper presents evidence that the 
waves in the Imagery reflect changes in the vert- 
ical tesoerature and moisture distribution pro- 
duced by the vertical oscillations. These features 
on the downwind or lee side on the mountains are 
called lee waves. 

2. BACKGROUND 

2 . I Instrument Characteristics 

The HAMS is a multlspeccral scanner capable 
of producing high resolution imagery In visible 
and infrared bands with detailed spectral resolu- 
tion. Table i summaries the characteristics of the 
instrument and its spectral banns. The MAMS is 
flown on a NASA U2/ER2 at an altitude of 20km pro- 
viding a resolution cell of 50-10(te on the ground. 


Table r~~ 


MAMS: The Hultispectral Atmospheric Mapping Seni or 

Provides high resolution visible and 
IR scanning capabilities from a high 
altitude aircraft platform. 


Instantaneous fisld-of-view 
Crouad resolution (nadir ,19. 8km) 
Total fleld-of-viev 
Roll correction ~ 

Calibre cion sources 


or u.^u rps 
2.5 or 5.0 arid 
50m or 100m 
86 ° 


Pixels per scan line 


n s* 

IK, 2 controllable bbs 
716 


Spectral Bands 


Channel Wavelength Spectral Region 



(microns) 

_ ' ' 

2.5/5. 0(36. 25rpa 

1 

.42 - .45 

blue 

1.0 

2 

.45 - .52 

blue/ green 

.5 

3 

.52 - .60 

green 

.4 

4 

.60 - ,62 

yellow/ orange 

.5 

5 

.63 - .69 

red 

.5 

6 

.69 - .75 

red/near IR 

.4 

7 

.76 - .90 

near IR 

. 4 

8 

.90 -1.05 

near IR 

.4 

9 

6.20 - 6.90 

water vapor 

.40/. 20 

10 

10.32 - 12.02 

window (cold) 

.15/. 10 

11 

10.32 - 12.02 

window fwara; 

.15/. 10 

12 

12.20 - 12.56 

water vapor 

.75/. 40 

[1) NEAT (nr/s ter/mVcm 1 ) for 1-8, 

NEAT CC) for 9-L2 


The width of the entire cross crack f ield-of-view 
is 40ka. providing detailed resolution over e rel- 
atively large area* 

Eight channels of the instrument sre used to 
aspirate the visible end near- infrared radiation • 
into narrow bands which sense reflected solar energy 
^ rom the earth's surface sod clouds. The three 
infrared channels (one of which is redundant) cover 
selected regions where infrared emission from the 
esrth and atnoaphsre are of interest. Thus bands 
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Figure 1. 


HAMS water vaoor and visible channel conooaice 1 mi terr from fmir j , , 

22 ,99, 1942-2109 «T. th. JZ . “ 

Nevaaa mountains over which the imagery was collected. tne blerr * 


include the mid-eropoaoneric (channel 9, 6.5 microns; 
ana low-tropospheric (channel 12, 12.3 microns) water 
vaoor region, and a window region (channel 10/11, 

11 microns) where the atmosohere is relatively trans- 
parent. The 11 and 12 micron bands lorn what has 
become known as the spilt window channels ( Chesters 
» 1983). Further details about the inacruaent 
and spectral bands may be round in Jedlovec et al. 
(1986). 

2 • 2 Structure of Lee Waves 

Stationary waves such ss mountain or lee wavee 
are a local phenomena vnich may have significant 
et feet on local weather conditions because of their 
associated cloud patterns and wind distributions. 

They may also feedback to the larger scalee of mo- 
tion. The presence of either mountain or lee waves 


is due to vertical oscillations of sir parcels as 
air is forced over the mountain. Mountain waves 
usually consist of one vertical oscillation of the 
flow pattern between two mountain peaks and are 
mainly governed in size by the spacing of the ridges. 
Lee waves, although forced by mountains, are not 
restricted in eize to one or two vertical oacil^ . 

1st ions as with ths mountain waves. In fact, the 
overshooting of the downward flow (ths restoring 
force is gravity) on the lee tide of the obstacle 
produces vertical oscillations which may extend for 
many wavelengths downstream. 

Tha water vapor image in Pig. 1 shows wave 
structures In ths lee of the Sierra Nevada moun- 
tains. Thrae or four wavea are dlacernahle, each 
having an estimated wavelength of 13-20km. The 
wavee extend away from tha fsca of tha Muncaln 
range for about 70km and along ths face for at Isaac 



120ka. The waves appear aomevhtt Irregular, poa- 
albly due co variations in the mountain peaks, tao- 
gencal distortion from the scanner geometry 
(Jedlovec et al . t 1986), limb brightening (due co 
increased water vapor path length at the edges of 
the flight track), as well as other unknown causes. 

Th® occurs of the vertical oscillation pro- 
ducing lee waves was studied in detail by Brunt (1927) 
who derived a set of equations which can be used to 
describe this flow pattern. From Brunt's work, a 
simple vertical oscillation in the atmosphere has 
a period given by 


consist of layers of varying stability and an ab- 
icnce of horizontal wind shear. In tha lowest 
layer, the acaosphsre would ba neutral with a 
•trong layer of itabillty aloft. Thia atable lay- 
er la ncceeeary for tha formation of tha wavea and 
■hould occur oaar tha level of the nountain peak. 

A layar of weafcar atabillty uaually exlata above 
thia level. The flow ahould ba perpendicular to 
tha nountain at all levela and Increaae in apeed 
with height. Significant horizontal ahear reduces 
the llklihood for tha deveiopaant of lea waves. 


0(c- 

where 

K is the environmental lapse rate, 

5* the dr y adiabatic lapse rate, 
g the gravitational acceleration, and 
T the environmental temperature. 

An air parcel moving horizontally with a velocity 
of V and under this vertical oscillation will 
complete one cycle over a distance of 




[-Hc-Op' 

Upon examining thia equation, it becomes apparent 
that the horizontal wavelength QJ is governed by 
the wind speed and atmospheric stability. Strong 
winds and low stability favor long wavelengths, 
while weak winds and strong stability favor shorter 
ones. 

A number of investigators have studied lee 
waves by observing their visible signatures and 
with direct measurements of temperature and wind 
velocity. Atkinson (1981) presents a fairly ex- 
tensive review of this research. TTie character- 
istics of these waves vary from occurance to oc- 
cur ance but representative values are presented 
below. The vavelengchs of lee waves range from 
2-70ka but most occurances are in the range of 
5-20km. The wavelengths tend to increase with 
height and seem to be related to the mean wind 
speed. Corby (1957) performed a linear regression 
of wind speed on wavelength and suggested that a 
high enough correlation existed for the wavelength 
to be determined from the wind speed alone. 

Another feature which is of interest is the 
wave amplitude, defined as the vertical distance 
between the peak and trough of a streamline. This 
Is more difficult to determine from the cloud it- 
self and requires a mincer of special measurements 
of tha wind field. From the available empirical 
data, these wave amplitudes arc typically about 
1km but amplitudes of 7km have been measured. The 
^ ar *? er amplitudes may occur when the wavelength 
closely matches the size of the obstacle forcing 
the wave. 

From these studies, it appears chat tvo 
P* ra *®ters control the appearance and size of lee 
waves: namely, stability and wind. Of course the 
obstacle (mountain) i3 necessary to initiate the 
vertical oscillation and mav have some control over 
the wave amniltude. Stability and wind are also 
important parameters in (1), however, the import- 
ance of stability in determining wavelength seems 
less important than indicated by (2). A typical 
environmental temperature end wind profile would 



■ 2.3 Atmospheric Condi tlona 

i 

In the previous section, the atmospheric 
conditions required for the development of lee waves 
were summarised. Although they are local conditions, 
they ere governed by the large scale atmospheric 
patterns. Thus, a discussion of tha large scale 
features on 22 January 1985 corresponding to the 
MAMS observations la necessary. The upper air 
data (Fig. 2) exhibited a pronounced ridge of high 
pressure over the eastern Pacific ocean and over 
the northwestern portion of the U. S. A trough was 
positioned over the central part of the country 
with a portion of it hanging back into Utah and the 
southern portions of Nevada and California. Thia 
* mt U P ■ prevailing flow pattern from the northeast 
over such of the Sierra Nevada mountains. 




Figure 2. a) SOOmb rawinaonde station plots and « 
height analysis, and b) 700nb plots and streamline* 
for 23 Jan 1985 at 0000 CHT. Plots and analyses 
•re la their standard units. 




The region of Interest (the Sierras) Is In 
a data void region of the r* win sonde observing net- 
work. However, since only 3h difference exists 
between the HA MS and radiosonde observations, and 
given che etaclonarlty of the weather systems, a 
reasonable estimation of the temperature, moisture, 
and wind profile may be available from the upwind 
sites. For this reason, data from two Nevada loca- 
tions la presented in Fig. 3. A moderate to strong 
inversion exists around 700sb In both soundings. 

The atflK>sphere Is quite dry above and below this 
level. The nearly saturated layer at the level of 
the inversion Indicates a shallow layer of clouds. 
Clouds are noc observed over the Sierras so the 
actual dewpoint will be somewhat colder. The winds 
are weaker at 72486 but both profiles Indicate 
flow from che mortheaac without much directional 
shear with height. From these diagrams, it seems 
reasonable to assume that conditions would be 
favorable for lee waves over che Sierras. 


3. DISCUSSION 


Although Fig. 1 was briefly discussed In 
previous sections, further discussion seems appro- 



Flgure 3. Skevr-LogP diagrams for two Nevada 
rawinsonde sites a) 72486, and b) 72383 at 0000 
CMT 23 Jan 1963. Winds arc plotted in knots. 


prlate at this time. The MAMS 6.5 micron and vis- 
ible /near- Infrared channel imagery la displayed 
from four adjacent flight tracks on 22 January 1985. 
The data were collected from 1942-2109 CMT over the 
plateau region of Nevada, the Sierra Nevada moun- 
tains and the central valley of California as in- 
dicated by the inset map. The peaks of the Sierras 
climb to about 4000m and drop sharply to the south- 
vest. Each flight crack la about 40km vide making 
the entire ground distance of the vmst-eesc extent 
160km. The north-south extent la roughly 400km. 

For display purposes, che resolution has been de- 
graded to about 400m. 

Many very small scale features ara apparent 
in the water vapor imagery. In the center of this 
imagery, alternating bright (cold/molat) and dark 
(warm/ dry) features exist in vave-llka patterns. < 
They ere orientated In a northwest -southeast dir- 
ection parallel to the mountains. It la not pos- 
ible to determine the position of the waves with 
respect to the position of the mountains from this 
Image alone. The accompanying flight tracks of 
visible data capture the cloud/land features at ! 
the same Instant as the water vapor imagery. 

Therefore thla data can be overlaid and used to 
determine the exact position of the waves. Several 
landmarks are visible in the scene, the beet being 
the quasi-circular Lake Mono In the right hand 
flight track. The bright regions in the lower por- 
tion of the vlaible image are fog and low clouds 
over the central valley. The bright areas In the 
center represent snow In the mountains above the 
2400m level. Thus the snow region can be used to 
outline the major mountain features. From this 
and the terrain map. It la clear that moat of the 
water vapor wave features reside on the downwind 
aide of the mountains where surface elevations 
are below 2000m. It should also be noted that no 
clouds are present in the visible imagery over 
this region. 

The estimated wavelengths of these vsvea 
(described in section 2.2) can be compared to that 
obtained from (2) or estimated from the mean tropo- 
spheric wind speed. The environmental lapse rate, 
temperature, and wind speed used in (2) is somevnat 
arbitrary. A value of Y of O a K/km (T*250*K), ana 
V of 10 m/a, yields 3.2 km, a rather short wave- 
length. Since the wave structures occupy more 
than a single level in the vertical, it may be more 
appropriate to usee valuea representative of a 
layer of the atmosphere. Using values from the 
800-40(kri> layer (T»25(f, t = 6.5, and V=20) t a wave- 
length of 10.7 km is obtained. This is somevnat 
less than che estimated length of 15-20km. As pre- 
viously mentioned, Corby (1957) used an empirical, 
relationship between wind speed and wavelength 
with a good correlation. The linear fit of his 
data (wavelength versue wind speed) produced a slope 
of roughly 0.45. For a mean environmental wind 
speed of 2CWs, a wavelength of 9km is predicted. 

In either case, a stronger wind speed would be need- 
ed to match the one estimated from che water vapor 
imagery. Considering chat the actual wind profile 
over the Sierras la somevnat uncertain, these 
calculations are relatively close to the estimates 
from che observed data. 

It la much more difficult to estimate the am- 
plitude of these lee wave features without direct 
in situ measurements • However, the next section 
describes an approach to estimate the amplitude 
by slmulatm* che HAMS water vapor channel radiances. 
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RADIATIVE TRANSFER SIMULATIONS 

The water vapor iatgi presented in Fig. I 
represents energy evicted by atmospheric water va- 
por somewhere in the aid- troposphere. The pro- 
cedure to determine this layer or level has been 
the subject of Investigation by researchers for 
many years. In the cose of cloud height determina- 
tions using window channel radiances, the altitude 
of the cloud can be estimated by comparing the 
measured blackbody temperature to a local ravinsonde. 
The water vapor signal does not come from a single 
level but actually a layer of the atmosphere 
(Jedlovec, 1984), and depends on the amount, ver- 
tical distribution, and temperature of the vapor 
Itself. It is unclear whether blackbody temper- 
atures in this channel can be used directly to in- 
fer a mean height from which the signal originates. 

In an attempt to estimate the amplitude of 
the lee waves observed in Fig. 1, several approaches 
were considered. Table 2 presents the observed 
blackbody temperatures measure in the peeks and 
troughs of the lee wave features. Typical values 
indicate a 6* temperature difference from the peak 
to the trough of the wave. It is difficult to use 
one of the soundings in Fig. 3 to directly Infer 
the levels of the peak and trough of the lee waves. 

It may be reasonable to assume that the level of 
the neutral point of the wave is 246* Jt (-27*C) and 
the peak and trough occurlng at the level of their 
measured temperature (Table 2), however. If this 
Is done, the estimated wove amplitude is about 1km. 
While this is consistent with the amplitudes of a 
typical lee wave, it depends heavily on the lapse 
rate of the founding used. Because of the lifting 
of the lover level Inversion by the mountains, a 
smaller lapse rate of temperature probably existed 
over the Sierras, producing a greater vertical 
amplitude for the lee wave. 

Another approach to estimate the lee wave 
amplitude used a composite of the two soundings in 
Fig. 3. From the composite flouncing, MAMS water 
vapor channel radiances were simulated (McMillan 
and Fleming, 1976; Fleming and McMillan, 1977; and 
Weinreb and Neuendorf fer, 1973). In order to ob- 
tain a representative profile In both the peak and 
trough of the waves, the inversion layer near 700mb 
and the accompanying moisture distribution were 
raised and lowered adlabatically in increments of 
1 km. Blackbody temperatures were computed from 


Table 2 

Observed and simulated blackboov temperatures 
for the lee waves. Units are in degrees'* K. 


Observed Blackbody Temperatures 

mean extreme 

Wave peak 243.5 242.0 

Wave crougn 249.5 251.0 

Simulated Blackbody Temperatures (4ka) 

243.9 
249.5 


the channel radiances. Temperature differences 
of 6*K were observed (similar to the measure values 
from Fig. 1) as veil as the appropriate blackbody 
temperatures vtimn the inversion was at 2 and 6km 
level for the trough and peak, respectively. This 
4km vertical amplitude, while larger than most lee 
wave amplitudes. Is mors conslstsnt with the level 
of the mounts in peaks which produce them. 

5. CONCLUSIONS 

i 

The sbovs work has presented evidence for i 
the existence of lee waves in the absence of clouds 
using high resolution MAMS moisture imagery. Al- 
, though verifying measurements do not exist, on eval- 
luatlon of wind profiles and stability in the sur- 
rounding region support the existence of these waves. 
Estimations of the horlxontal wavelengths of the 
, observed lee wave features are not inconsistent 
with those from previous studies. The estimated 
; vertical amplitude, although less certain, seems • 
realistic for the setting. Further work Is neces- 
sary to mors accurately determine this vertical 
amplitude. 
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1 . INTRODUCTION 

Accurate absolute calibration of radiometric data is an important 
component of useful remotely sensed data. Lack of absolute accuracy in 
the past has hindered the effectiveness of satellite radiances in budget 
studies, data assimilation efforts, and climatology studies, which can 
suffer from long term temporal instability of the calibration. Because 
of this problem, many efforts confine conclusions to relative variations 
only (e.g. Chesters et al. 1983). It is desirable that the absolute 
calibration of new and future satellite platform instruments be 
evaluated. Evaluation may be made by directly comparing collocated 
fields of view between two (or more) instruments. While a favorable 
comparison of this sort does not necessarily imply good absolute 
calibration, agreement between independent instruments tends to indicate 
that they are fairly well calibrated. 

The High resolution Interferometer Sounder (HIS) and the 
Multispectral Atmospheric Mapping Sensor (MAMS) can be directly compared 
with each other for absolute calibration accuracy. The MAMS contains 3 
infrared channels that are simulated by convoluting HIS high spectral 
resolution radiances with the MAMS spectral response functions. HIS and 
MAMS flew together. during the cooperative Huntsville Meteorological 
Experiment (COHMEX) during the summer of 1986. This note presents the 
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results of the absolute calibration study between MAMS and HIS 6jxm, 11/im 
and 12 pm data of June 15, 1986. 

2. INSTRUMENT DESCRIPTIONS 

A. MAMS 

The MAMS instrument is an aircraft borne 12 channel spectrometer 
which collects 100m resolution data across a swath of about 37 km. 

During COHMEX, the MAMS configuration collected 8 bit precision data in 
8 visible and 4 infrared channels which include 6pm, 11/im, and 12pm 
bands similar to those of the VISSR Atmospheric Sounder (VAS) (Menzel et 
al. 1986). The infrared data is absolutely calibrated on each scan by 
viewing two onboard blackbodies of known temperature. Calibration noise 
is reduced through multiple viewings of each blackbody during a given 
scan. A linear relationship between blackbody radiance and digital 
count values is assumed. MAMS has been extensively discussed in 
Jedlovec et al. (1989). 

B. HIS 

The airborne High-resolution Interferometer Sounder (HIS) is a 
Michelson interferometer which measures thermal emission from the 
^^rth’s surface and clouds and atmospheric constituents in longwave 
(151m), midwave (6.71m) and shortwave (4.31m) bands (Fig. 1). Table 1 
(Smith et al. 1983) summarizes the characteristics of the HIS aircraft 
instrument. With its high spectral resolution (< 1 cm' 1 ), HIS has 
several thousand spectral channels; many more than the ten or twenty 
spectral channels on current filter radiometers such as VAS and HIRS/2. 


Detection in each of three spectral bands is performed using arsenic- 
doped silicon detectors contained in a single detector/dewar assembly. 
Three detector/filter systems are used to optimize the signal- to-noise 
ratio of the observations throughout the broad spectral region measured. 
High emissivity, temperature controlled reference blackbodies are used 
for an absolute calibration. Every two minutes of HIS observations 
consist of four cold blackbody (CBB) views, four hot blackbody (HBB) 
views, and twelve earth views. The four blackbody views are averaged 
for each blackbody and used for the calibration of the earth views. The 
basic approach for determining absolute radiances from the HIS nadir- 
viewing Michelson interferometer is the same as that used for filter 
radiometers and has been used successfully for other interferometric 
applications (Revercomb et al. 1988). The HIS is designed to yield an 
output which is linear in the incident radiance for all channels in the 
optical path difference of the instrument. The averaged HBB and CBB 
temperatures and counts determine the slope and offset which define the 
linear instrument response at each wave number (Huang 1989). The slope, 
S t , is given by 



where V t E is the count value in the earth scene of channel t, and V t H 
C 

and are the count values In the HBB and CBB scene respectively of 
channel t. The earth scene radiance at channel t, R t ^, is 

R t E - S t * (B t H - B t C ) + B t C 
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where B t is the Planck function of channel t. More detailed and 
complete calibration procedures concerning the effects of phase 
ambiguity, finite field of view, and other instrument effects on 
radiance calibration may be found in Revercomb et al. (1988). Figure 2 
shows the noise spectra for the June 15, 1986 of COHMEX. In general, 
good signal-to-noise ratios have been realized In COHMEX HIS flights, 
with an average noise level of 0.2 - 0.3 K for all three bands. 

3. HIS AND MAMS COMPARISONS 

A. COHMEX data set 

The MAMS and HIS instruments were flown repeatedly during COHMEX, 
which was conducted during the summer of 1986 over a portion of the mid- 
eastern United States. On June 15, MAMS and HIS were flown together 
from 1730 GMT to 2100 GMT on the NASA U2 aircraft over the central 
Tennessee / northern Alabama region. The flight path is shown in Figure 
3. The flight region, a region of highly variant land uses and types, 
was partially obscured by a scattered field of low level cumulus. 

For further information on COHMEX and on the weather of June 15, the 
reader is referred to Williams et al. (1987). 

B. Comparison procedure 

In order to properly compare the MAMS and HIS radiances the data 
sets were equated in factors of footprint registration and size, time of 
observation, viewing angle, and spectral response of channels. Since 
MAMS and HIS flew together on June 15, footprint registration, 
observation time, and viewing angle considerations were easily accounted 
for. However, because of different data recording rates, it is unlikely 
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that the two data sets are exactly collocated in time and space. It is 
estimated though that collocation errors are less than .5 km, which is 
less than a quarter of a single HIS field of view. Instrument 
resolution differences were accounted for by averaging the observed MAMS 
data over the size of a HIS 4 x 2 km averaged footprint (2 fovs). The 
HIS and MAMS instrument impulse functions were both assumed to be 
symmetrical. To equate the spectral channels, MAMS spectral response 
functions were convoluted with two averaged HIS fov (4x2 km) radiances 
to produce "MAMS equivalent” channels from the HIS radiances. MAMS 
spectral response function data is given at 1 wavenumber resolution, 
comparable to the HIS spectral resolution. The comparison was then made 
between the spatially averaged MAMS radiances and the HIS radiances that 
had been used to produce "MAMS equivalent” channels. 

C. Results of calibration comparison 

The results of comparing the MAMS and HIS equated radiances are 
shown in Figure 4 for the 6^m, ll^m and 12^m MAMS channels. Statistics 
of the comparisons are shown in Table 2. 


Table 2. - MAMS/HIS comparison statistics for June 15, 1986 (354 
collocations) . 


channel 

6fim 

11/im 

12/im 


bias(k) rms(k) correlation 

-0.03 0.88 0.98 

0.93 3.30 0.94 

-2.41 3.48 0.95 


High correlations in Table 2 indicate excellent agreement in relative 
variation between the two instruments. This lends credibility to the 
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collocation accuracy between the two instruments. This agreement is 
easily seen in the 6/m comparison of figure 4 but less so in the ll/m 
and 12/m comparisons because of the high thermal variability of the 
scenes. The agreement is more readily apparent in the ll/m and 12/m 
data when the U2 aircraft flew over a cumulonimbus cloud, causing both 
the MAMS and HIS temperatures to drop sharply. In the absolute 
calibration comparison, Che 6/m data shows excellent agreement between 
the two instruments with a total flight bias of only .03 K. The ll/m 
data also shows good results with a bias of less than 1 K. The absolute 
calibration of the 12/m data shows less favorable agreement. It Is 
noteworthy that the ll/m and 12/m biases and RMS may be affected by some 
scene misregistration error, especially when the aircraft went into and 
came out of turns. In order to test the effect of possible 
misregistration errors, a filtered data set of comparisons was generated 
by removing comparisons that exceeded a 5 K threshold difference between 
the two instruments. This should exclude most misregistered 
comparisons. The statistics for 269 samples are shown in Table 3. 


Table 3. - Misregistration filtered MAMS/HIS comparisons for June 15, 
1986. (269 samples). 

channel bias(k) rms(k) correlation 

6/m -0.03 0.78 0.98 

ll/m 1.10 2.05 0.98 

12/m -2.27 2.62 0.98 


The biases of the comparisons are somewhat changed, though not 
substantially. This indicates that collocation errors are not affecting 
the bias very strongly, suggesting that collocation errors cancel out. 
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That the 12pm bias remains significantly larger than the 11pm bias is 
probably due to the smaller bandwidth and thus noisier data in that 
channel. The reduction in RMS reflects the removal of the misregistered 
comparisons. Note that the smallest effect was on the 6pm RMS. This is 
not surprising since the 6pm spatial variability is relatively weak at 
small scales, and there was no cloud contamination of the 6pm data. In 
other words, small scale misregistration errors should not have a 
significant effect on the 6pm MAMS/HIS comparison. Also note that the 
12pm bias now accounts for almost all of the RMS while considerable 
unexplained RMS remains in the 6pm and 11pm statistics of Table 3. An 
investigation of MAMS data showed that the calibration "floated" with 
time in the 6pm and 11pm channels. This probably accounts for the 
unexplained RMS and is likely also present in the 12pm MAMS data though 
it may be masked by the noisier nature of the 12pm data. This 
calibration float is an undesirable characteristic of the MAMS June 15, 
1986 data, however improvements to the instrument since 1986 have 
probably improved the performance significantly (Jedlovec et al . 1989). 
Also, since the calibration float introduces an uncertainty to the bias, 
the RMS (which includes both calibration float and true bias) serves as 
an upper limit to the true bias. Thus, in the worst case (and assuming 
that none of the remaining RMS of Table 3 is due to misregistration) , 
the bias of the 11pm comparison might approach 2 K (comparable to 12pm 
bias) while the 6pm bias is still less than 1 K. 

This analysis suggests that the biases of Table 1 are the result 
of absolute calibration disagreement between the two instruments. This 
calibration disagreement may be caused by the choice of blackbody 
temperatures for data collection. MAMS blackbody temperatures were 
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approximately 243 K and 273 K (HIS were approximately 240 K and 300 K) , 
thus forcing a significant linear extrapolation to the thermal scene 
temperatures (-300 K) of the MAMS 11pm and 12pm data on June 15. Linear 
extrapolation errors for AVHRR 11pm data exceed 1 K when scene 
temperatures exceed the warm blackbody temperature by 10 K (Lauritson et 
al. 1979 and subsequent revisions to that document). While an 
extrapolation correction has not been determined for the MAMS channels, 
it is suggested that the MAMS/HIS bias is at least partly a function of 
the choice of blackbody temperatures. On the other hand, the 6pm MAMS 
data, which Is essentially unaffected by cloud cover or substantial 
small scale thermal variance, requires almost no extrapolation (scene 
temperatures - 245 K) , and thus would be expected to be well calibrated. 
The HIS cold blackbody being set at 240 K makes the HIS calibration also 
very reliable in the 6pm region. Thus, of the three compared channels, 
the 6pm comparison should be the best calibrated, and indeed shows the 
best results. This result indicates that when blackbody temperatures 
are set to encompass the expected scene temperature, that good absolute 
accuracy of calibration is achieved. 



4. CONCLUSIONS - 

An absolute comparison of HIS and MAMS radiances has been 
performed in an effort to evaluate the absolute calibration of each 
instrument. While not conclusive, good agreement between the two 
independent instruments indicates that both instruments are well 
calibrated. Results show excellent agreement between the 6pm and 11pm 
data (bias < IK) with lessor agreement in the 12pm data (bias - 2K) . 
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Bias results suggest that selecting blackbody temperatures close to the 


expected scene temperatures improves absolute calibration. 
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CHARACTERISTICS OF THE HIS AIRCRAFT INSTRUMENT 


Spectral range (cm -1 ) : 

Band I 
Band II 
Band III 

Field of view diameter (mr) : 
Telescope 
Interferometer 

BlackBody reference sources : 
Emissivity 

Aperture diameter (cm) 
Temperature stability (K) 


590 - 1070 
1040 - 1930 
2070 - 2750 


100 

30 


>0.998 

1.5 

± 0.1 


Auto-aligned Interferometer : modified BOMEM BBDA2.1 


Beamsplitter : 

Substrate 

Coatings 

Maximum Delay (double-sided) (cm) 
Band I 

Band II & ill 


KCL 

Ge+Sb2S ^ 


± 1.8 

+ 1 . 2 , - 0.8 


Michelson mirror optical scan rate (cm/s) : 0. 6-1.0 


Aperture stop (at exit window) : 

Diameter (cm) 4.1 

Central obscuration area fraction 0.17 
Area (cm 2 ) 10 .s 


Area-solid angle product (cm 2 -sr) 


0.0076 


Detectors : 

TyP e Ar doped Si 

Diameter (cm) 0.16 

Temperature (K) 6 


Table 1, characteristics of the HIS aircraft instrument 







Figure 1, Brightness temperature spectra measured by HIS 

from NASA U2 aircraft over Huntsville, Alabama on June 
15 , ] 
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enl Ire flight of June 15, 1986. 






